The Role of Symmetry and Separation in Surface Evolution and Curve Shortening
With few exceptions, known explicit solutions of the curve shortening flow (CSE) of a plane curve, can be constructed by classical Lie point symmetry reductions or by functional separation of variables. One of the functionally separated solutions is the exact curve shortening flow of a closed, conve...
Збережено в:
| Дата: | 2011 |
|---|---|
| Автори: | , |
| Формат: | Стаття |
| Мова: | English |
| Опубліковано: |
Інститут математики НАН України
2011
|
| Назва видання: | Symmetry, Integrability and Geometry: Methods and Applications |
| Онлайн доступ: | http://dspace.nbuv.gov.ua/handle/123456789/147169 |
| Теги: |
Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
|
| Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| Цитувати: | The Role of Symmetry and Separation in Surface Evolution and Curve Shortening / P. Broadbridge, P. Vassiliou // Symmetry, Integrability and Geometry: Methods and Applications. — 2011. — Т. 7. — Бібліогр.: 28 назв. — англ. |
Репозитарії
Digital Library of Periodicals of National Academy of Sciences of Ukraine| id |
irk-123456789-147169 |
|---|---|
| record_format |
dspace |
| spelling |
irk-123456789-1471692019-02-14T01:23:27Z The Role of Symmetry and Separation in Surface Evolution and Curve Shortening Broadbridge, P. Vassiliou, P. With few exceptions, known explicit solutions of the curve shortening flow (CSE) of a plane curve, can be constructed by classical Lie point symmetry reductions or by functional separation of variables. One of the functionally separated solutions is the exact curve shortening flow of a closed, convex ''oval''-shaped curve and another is the smoothing of an initial periodic curve that is close to a square wave. The types of anisotropic evaporation coefficient are found for which the evaporation-condensation evolution does or does not have solutions that are analogous to the basic solutions of the CSE, namely the grim reaper travelling wave, the homothetic shrinking closed curve and the homothetically expanding grain boundary groove. Using equivalence classes of anisotropic diffusion equations, it is shown that physical models of evaporation-condensation must have a diffusivity function that decreases as the inverse square of large slope. Some exact separated solutions are constructed for physically consistent anisotropic diffusion equations. 2011 Article The Role of Symmetry and Separation in Surface Evolution and Curve Shortening / P. Broadbridge, P. Vassiliou // Symmetry, Integrability and Geometry: Methods and Applications. — 2011. — Т. 7. — Бібліогр.: 28 назв. — англ. 1815-0659 2010 Mathematics Subject Classification: 35A30; 35K55; 58J70; 74E10 DOI:10.3842/SIGMA.2011.052 http://dspace.nbuv.gov.ua/handle/123456789/147169 en Symmetry, Integrability and Geometry: Methods and Applications Інститут математики НАН України |
| institution |
Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| collection |
DSpace DC |
| language |
English |
| description |
With few exceptions, known explicit solutions of the curve shortening flow (CSE) of a plane curve, can be constructed by classical Lie point symmetry reductions or by functional separation of variables. One of the functionally separated solutions is the exact curve shortening flow of a closed, convex ''oval''-shaped curve and another is the smoothing of an initial periodic curve that is close to a square wave. The types of anisotropic evaporation coefficient are found for which the evaporation-condensation evolution does or does not have solutions that are analogous to the basic solutions of the CSE, namely the grim reaper travelling wave, the homothetic shrinking closed curve and the homothetically expanding grain boundary groove. Using equivalence classes of anisotropic diffusion equations, it is shown that physical models of evaporation-condensation must have a diffusivity function that decreases as the inverse square of large slope. Some exact separated solutions are constructed for physically consistent anisotropic diffusion equations. |
| format |
Article |
| author |
Broadbridge, P. Vassiliou, P. |
| spellingShingle |
Broadbridge, P. Vassiliou, P. The Role of Symmetry and Separation in Surface Evolution and Curve Shortening Symmetry, Integrability and Geometry: Methods and Applications |
| author_facet |
Broadbridge, P. Vassiliou, P. |
| author_sort |
Broadbridge, P. |
| title |
The Role of Symmetry and Separation in Surface Evolution and Curve Shortening |
| title_short |
The Role of Symmetry and Separation in Surface Evolution and Curve Shortening |
| title_full |
The Role of Symmetry and Separation in Surface Evolution and Curve Shortening |
| title_fullStr |
The Role of Symmetry and Separation in Surface Evolution and Curve Shortening |
| title_full_unstemmed |
The Role of Symmetry and Separation in Surface Evolution and Curve Shortening |
| title_sort |
role of symmetry and separation in surface evolution and curve shortening |
| publisher |
Інститут математики НАН України |
| publishDate |
2011 |
| url |
http://dspace.nbuv.gov.ua/handle/123456789/147169 |
| citation_txt |
The Role of Symmetry and Separation in Surface Evolution and Curve Shortening / P. Broadbridge, P. Vassiliou // Symmetry, Integrability and Geometry: Methods and Applications. — 2011. — Т. 7. — Бібліогр.: 28 назв. — англ. |
| series |
Symmetry, Integrability and Geometry: Methods and Applications |
| work_keys_str_mv |
AT broadbridgep theroleofsymmetryandseparationinsurfaceevolutionandcurveshortening AT vassilioup theroleofsymmetryandseparationinsurfaceevolutionandcurveshortening AT broadbridgep roleofsymmetryandseparationinsurfaceevolutionandcurveshortening AT vassilioup roleofsymmetryandseparationinsurfaceevolutionandcurveshortening |
| first_indexed |
2025-07-11T01:31:26Z |
| last_indexed |
2025-07-11T01:31:26Z |
| _version_ |
1837312238881341440 |
| fulltext |
Symmetry, Integrability and Geometry: Methods and Applications SIGMA 7 (2011), 052, 19 pages
The Role of Symmetry and Separation
in Surface Evolution and Curve Shortening?
Philip BROADBRIDGE † and Peter VASSILIOU ‡
† School of Engineering and Mathematical Sciences, La Trobe University,
Melbourne, Victoria, Australia
E-mail: P.Broadbridge@latrobe.edu.au
‡ Faculty of Information Sciences and Engineering, University of Canberra,
Canberra, A.C.T., Australia
E-mail: Peter.Vassiliou@canberra.edu.au
Received January 23, 2011, in final form May 25, 2011; Published online June 01, 2011
doi:10.3842/SIGMA.2011.052
Abstract. With few exceptions, known explicit solutions of the curve shortening flow
(CSE) of a plane curve, can be constructed by classical Lie point symmetry reductions or by
functional separation of variables. One of the functionally separated solutions is the exact
curve shortening flow of a closed, convex “oval”-shaped curve and another is the smoothing of
an initial periodic curve that is close to a square wave. The types of anisotropic evaporation
coefficient are found for which the evaporation-condensation evolution does or does not
have solutions that are analogous to the basic solutions of the CSE, namely the grim reaper
travelling wave, the homothetic shrinking closed curve and the homothetically expanding
grain boundary groove. Using equivalence classes of anisotropic diffusion equations, it is
shown that physical models of evaporation-condensation must have a diffusivity function
that decreases as the inverse square of large slope. Some exact separated solutions are
constructed for physically consistent anisotropic diffusion equations.
Key words: curve shortening flow; exact solutions; symmetry; separation of variables
2010 Mathematics Subject Classification: 35A30; 35K55; 58J70; 74E10
1 Introduction
The standard curve shortening flow is a nonlinear evolution by curvature,
∂γ
∂t
(x, t) = κ(x, t)n(x, t), γ(x, 0) = γ0(x), (1)
where κ(x, t) is the Euclidean curvature of γ(x, t) and n(x, t) is the outward unit normal to the
curve at each point. For the curve-shortening flow parameterized as x 7→ (x, y(x, t)),
n · eyyt = n · ∂γ
∂t
(x, t), yt =
yxx
1 + y2
x
, (2)
which describes the curve shortening flow in Cartesian coordinates. This equation arises in the
practical context of metal surface evolution [24]. More recently, it has been used extensively as
an isotropic curve-smoothing mechanism in image processing [23, 25]; see also [9]. Equation (1)
is well known and has been the subject of several important studies, for instance [16, 17]. Few
explicit and exact curve shortening flows of plane curves are known in the literature. Some
?This paper is a contribution to the Special Issue “Symmetry, Separation, Super-integrability and Special
Functions (S4)”. The full collection is available at http://www.emis.de/journals/SIGMA/S4.html
mailto:P.Broadbridge@latrobe.edu.au
mailto:Peter.Vassiliou@canberra.edu.au
http://dx.doi.org/10.3842/SIGMA.2011.052
http://www.emis.de/journals/SIGMA/S4.html
2 P. Broadbridge and P. Vassiliou
similarity solutions with prescribed-slope boundary conditions have been constructed paramet-
rically in terms of integrals of algebraic functions of elementary functions [5, 20]. King [21] writes
down a number non-invariant solutions. In addition, there are known to be various self-similar
rotating spiral and flower-head solutions [1, 18].
We demonstrate in this paper how for most solutions known to us, Lie symmetries of various
types or separable coordinate systems play a vital role and we conjecture that all currently
known solutions and many new solutions can be obtained via a series of higher order constraints
in the method of functional separation. In Section 2 we briefly review the solutions that can
be obtained by symmetry reduction. In Section 3, we show how the method of functional
separation of variables recovers additional interesting explicit, non-self-similar solutions. Some of
the properties and materials science applications of these solutions are developed in more detail.
In Section 4, an anisotropic version of the second-order Mullins equation is derived for materials
that include a realistic dependence of evaporation coefficient on surface orientation. Equivalence
classes for anisotropic evaporation coefficients are constructed under the Euclidean group. Using
the equivalence classes, in Section 5 we investigate under which types of anisotropy, analogs of the
standard solutions for isotropic diffusion, do or do not exist. Some examples of exact solutions
are produced by direct construction and in Section 6 by functional separation of variables.
2 Symmetry reductions of the curve shortening equation
One of the most widely used techniques for constructing explicit, exact solutions of nonlinear
partial differential equations is Lie symmetry reduction. This has been extensively investigated
in relation to the partial differential equation (2) or its derivative forms
ut = ∂x[D(u)ux], D(u) = 1/
(
1 + u2
)
, (3)
κt = κ2κθθ + κ3, (4)
Rt = ∂θ[R
−2Rθ]−R−1, (5)
where u = yx, κ is curvature, R = 1/κ and θ is the orientation angle along a convex curve [2].
Although the standard nonlinear diffusion equation (3) is rarely used in the context of curve
shortening, it has non-trivial Lie potential symmetries that enable one to construct exact simi-
larity solutions [4]. The reaction-diffusion equation in standard form (5), another equation that
is rarely used in the context of curve shortening, has been fully classified by Lie point symmetry
reductions [15].
In common with all autonomous nonlinear diffusion equations of second order, the curve
shortening equation (2) is invariant under translations in x, y and t, plus the Boltzmann scaling
group generated by x∂/∂x+ y∂y + 2t∂/∂t:
(x̄, ȳ, t̄) =
(
eεx, eεy, e2εt
)
.
This allows the possibility of two types of scale-invariant similarity solution: the expanding
solution of the form y/
√
t− t0 = G(x/
√
t− t0), (t > t0) and the contracting solution of the form
y/
√
t0 − t = G(x/
√
t0 − t), (t < t0). Being Euclidean-invariant, (2) has the rotation group as an
additional symmetry. Some interesting exact solutions to this equation, and consequently to (3)
can be constructed by consecutive symmetry reductions [4]; see also [11].
The expanding homothetic similarity solution with initial-boundary conditions yx(0, t) = m,
y → 0, x→∞ and y(x, 0) = 0, was given in [5]. The solution takes the form of the symmetrized
(upper) curve in Fig. 1, which represents an evolving grain boundary groove [24, 5].
Alternatively, the above solution may be extended smoothly to the interior of an obtuse-
angled wedge, as shown. These and the homothetic solutions in an acute-angled wedge, ob-
tainable from the same type of reduction, make up the “open-angle” solutions [20].
Symmetry and Separation in Surface and Curve Flows 3
Figure 1. Self-similar evolution by evaporation near a grain boundary groove (symmetrized upper curve)
or equivalently by condensation in an obtuse-angled wedge. Output times are at t = 1 and t = 4.
The homothetically shrinking simple closed-curve solution is the circle r =
√
2(t0 − t)1/2.
Self-intersecting closed-curve “flower-head” solutions were given in [1].
Steady-state solutions, invariant under time translations, are simply the straight lines. Like-
wise, the solutions that are invariant under spatial translations in a particular direction are
simply the straight lines in that direction. Another simple known solution is the travelling
wave solution of the form y − ct = G(x), which is invariant under the translation in space and
time generated by ∂/∂t + c∂/∂y. This is the well-known “Calabi grim reaper” [17], taking the
appearance of a shepherd’s crook
y − ct =
−1
c
log cos(cx).
Uniformly rotating solutions of the form r = G(φ−σt) (r and φ being plane polar coordinates)
follow from reduction by the symmetry generated by σ(−y∂/∂x+ x∂/∂y)− ∂/∂t.
Expanding or contracting rotating solutions of the invariant form r/
√
t = G(φ− σ
2a3
log t) or
r/
√
−t = G(φ − σ
2a3
log(−t)) follow from reductions under the one-parameter group generated
by a linear combination of generators of dilatations and rotations
Γ = Γ4 + σΓ5, Γ4 = x
∂
∂x
+ y
∂
∂y
+ 2t
∂
∂t
, Γ5 = −y ∂
∂x
+ x
∂
∂y
.
A classification of the types of solutions of the associated reduced ordinary differential equa-
tions, but without solving them explicitly and without referring to symmetries, was given re-
cently by Halldorsson [18].
2.1 Reciprocal transformations
It is well known that the class of nonlinear heat equations is stabilised by reciprocal transfor-
mations which can sometimes be used to generate new solutions from old [22]. A reciprocal
transformation may be viewed as a map of the graph of a solution u(x, t) of
ut = (D(u)ux)x
4 P. Broadbridge and P. Vassiliou
to that of a solution u′(x′, t′) of
u′t′ =
(
D′(u′)u′x′
)
x′
.
The map (x, t, u)→ (x′, t, u′) is defined by
dx′ = u dx+D(u)ux dt, dt′ = dt, u′ =
1
u
. (6)
It can be shown that this induces the transformation of diffusivities
D′(s) =
1
s2
D
(
1
s
)
. (7)
The diffusivity of interest here
D(u) =
1
1 + u2
is seen to be invariant under the reciprocal transformation,
D′(u′) =
1
1 + u′2
= D(u′).
This raises the possibility of using the known solutions and generating new solutions of (3) and
hence new curve flows by quadrature. However, it turns out that the reciprocal transformation
of solutions of (3) acts geometrically trivially on solutions of the curve shortening equation.
Proposition 1. Suppose u(x, t) is a solution of (3) and y(x, t) the corresponding solution of
the curve shortening equation. Let (x′, t′, u′) be the reciprocal transformation of (x, t, u). Then
the induced transformation on solutions of (2) satisfies (x′, t′, y′) = (y, t, x+ c), c = const. That
is, for each t, the reciprocal transformation of a known solution of (3) induces a reflection (up
to additive constant) of the solution y(x, t) in the line y = x.
Proof. Equation (3) and the reciprocal transformation (6) permits us to write
x′ =
∫ x
0
u(x1, t)dx1 −
∫ t
a
V0(t1)dt,
where a is a constant and V0(t) = −D(u(0, t))ux(0, t). In terms of Cartesian coordinates (x, y),
x′ =
∫ x
0
yx(x1, t)dx1 −
∫ t
a
V0(t1) = y(x, t)− y(0, t) +
∫ t
a
D(yx(0, t1))yxx(0, t1)dt1
= y(x, t)− y(0, t) +
∫ t
a
yt(0, t1)dt1 = y(x, t)− y(0, a).
Then
u′ = 1/u = 1/yx = xy = xx′ .
Now as y′ satisfies y′x′ = u′ just as yx = u, then it follows that y′ = x + H(t′). However, since
y′(x′, t′) is to satisfy the curve shortening equation
y′t′ =
y′x′x′
1 + (y′x′)
2
,
it follows that H = const. �
Symmetry and Separation in Surface and Curve Flows 5
Let f1 be a function such that the reflection of its graph defines a function f2. In general f1
and f2 will satisfy different differential equations and this can be very useful. For example,
the nonlinear diffusion equation with diffusivity D(s) = s−2, seen from (7) to be directly
transformable to the classical linear heat equation (k = 1), is one of many integrable equations
that can be linearized in its potential form by the hodograph transformation [12]: (x′, y′) = (y, x)
given in Proposition 1. However, from a geometric point of view a curve and its reflection are
indistiguishable. Thus invariance of (3) under reciprocal transformations does not increase the
class of known curve flows. This highlights the difficulty of constructing solutions of the curve
shortening equation and partly explains why so few interesting exact solutions are known.
3 Functionally separable nonlinear heat equations
Apart from symmetry methods, one can seek reductions that arise from second or higher order
differential constraints rather than the first order differential constraints implied by classical Lie
reduction. Unfortunately, there is no general procedure for seeking such reductions and much
effort has gone into devising new reduction strategies for partial differential equations. Thus
in Doyle and Vassiliou [14] the authors managed to classify all one-dimensional sourceless heat
equations
ut = (D(u)ux)x (8)
that admit separation of variables in some field variable ū. That is, one asks for a change of
field variable ū = m(u), such that the image of (8) under the change of variable, constrained by
the additively separable condition
ūxt = 0,
is a differential system of finite type that can therefore be solved by ordinary differential equa-
tions. Indeed the resulting differential system has the general form
ūt = f(ū)ūxx + g(ū)ū2
x, ūxt = 0, (9)
for some functions f > 0, g. It is proven in [14] that for any such pair f , g there is a change of
dependent variable that transforms (9)1 to (8). It turns out that system (9) has a maximal 3-
parameter solution space for any given pair f , g. In this manner the authors obtain exactly nine
distinct diffusivities D up to the maximal transformation group that preserves the canonical
form (8) for which the maximal 3-parameter solution space is achieved. In many cases the 3-
parameter solution was constructed. This considerably extends the list of nonlinear diffusion
equations for which explicit solutions are available. One equation on the Doyle–Vassiliou list is
the nonlinear heat equation
ut =
(
ux
1 + u2
)
x
. (10)
Thus the differential 1-form
ω = u dx+
ux
1 + u2
dt
is closed on solutions of (10). It is easy to see that any function y(x, t) satisfying dy = ω is
a solution of the curve shortening equation (2). The solutions of (10) constructed in [14] are the
functions
u(x, t) = U
(
σ(x+ a), σ2(t+ b)
)
,
6 P. Broadbridge and P. Vassiliou
where a, b and σ 6= 0 are arbitrary constants and U is one of the functions
U(x, t) = tanx,
U(x, t) =
x√
−x2 − 2t
,
U(x, t) =
±1√
e2(x−t) − 1
,
U(x, t) =
sinhx√
− cosh2 x− e−2t
, (11)
U(x, t) =
± coshx√
− sinh2 x+ e−2t
,
U(x, t) =
sinx√
cos2 x− e2t
,
U(x, t) =
sinx√
cos2 x+ e2t
.
Each solution (11) provides a 3-parameter family of explicit curve shortening flows except for
the complex-valued, (11)4.
The function y(x, t) may be obtained from U(x, t) simply by integrating in x, then adding
a suitable function of t.
1. (11)1 then leads to the Calabi “grim reaper” travelling wave.
2. (11)2 integrates to the well-known shrinking circle x2 + y2 = 2(t0 − t).
3. Integration of (11)3 merely produces a horizontal version of the vertical grim reaper.
4. (11)4 is complex-valued, not considered further in the current practical context.
5. (11)5 is equivalent to (11)7 by a rotation in the xy-plane.
6. The final two solutions are related to those previously presented by King [21], expressed
in the time-reversed form as examples of finger growth; these deserve closer inspection.
3.1 Exact heat flow of a convex curve
The studies of Gage–Hamilton [16] and Grayson [17] on flow by curvature of embedded plane
curves is a justly celebrated chapter in differential geometry.
Theorem 1 (Gage–Hamilton). Let γ0 : Ix → R2 be a convex curve embedded in the plane.
Let γ0 evolve by the curve shortening flow. That is,
∂γ
∂t
(x, t) = κ(x, t)n(x, t), γ(x, 0) = γ0(x),
where κ(x, t) is the Euclidean curvature of γ(x, t). Then the curve remains convex and becomes
circular as it shrinks in the sense that
1) the ratio of the inscribed radius to the circumscribed radius approaches 1;
2) the ratio of the maximum to the minimum curvature approaches 1;
3) the higher order derivatives of the curvature converge to zero uniformly.
(11)6 provides the only known example of an explicit, non-self similar curve shortening flow
in case the initial curve γ0 is a closed, convex embedded plane curve which is not a circle. We
use results described in Section 2. The sixth function in (11) is
u =
sinx√
cos2 x− e2t
. (12)
Symmetry and Separation in Surface and Curve Flows 7
In the case of (12), the corresponding solution of (2) is
y = t− ln
(
cosx+
√
cos2 x− e2t
)
. (13)
Clearly if y is a solution of (2) then so is −y. The two solutions join smoothly along y = 0 and
can be jointly expressed in the simple implicit form
cosh y − e−t cosx = 0. (14)
This solution is recorded in [21] and shown here to arise from functional separation. Solutions of
this type are also studied in [13] and referred to as ‘Angenent ovals’, although explicit solutions
are not written down in the latter paper. We will now study some properties of this solution. For
each t ∈ (−∞, 0) equation (14) defines a closed, convex “oval-shaped” curve which is symmetric
about the x- and y-axes for
x ∈
(
− cos−1
(
et
)
, cos−1
(
et
))
, y ∈
(
− cosh−1
(
e−t
)
, cosh−1
(
e−t
))
.
By analogy with an ellipse, eccentricity may be defined as
ε(t) =
√
1− (xmax/ymax)2 =
√
1−
[
cos−1 (et)
cosh−1 (e−t)
]2
=
√
2
3
|t|+O
(
|t|3
)
,
showing approach to circularity (ε → 0) as t approaches the extinction time 0. The curvature
at each point (x, t) of the curve is
κ(x, t) =
e−t cosx√
e−2t − 1
.
For each t ∈ (−∞, 0), the maximum curvature occurs at x = 0 with value κmax = e−t/
√
e−2t − 1,
while the minimum occurs at the extremities along the minor axis, x = ± cos−1
(
et
)
with value
κmin = 1/
√
e−2t − 1. Hence
lim
t→0−
κmax
κmin
= lim
t→0−
e−t = 1,
verifying the Gage–Hamilton theorem and more specifically showing that the ratio of maximum
to minimum curvature is an exponential function converging to unity. Of course, the curvature
itself is an unbounded function of time as the flow continues toward extinction.
It transpires that for this curve flow, arclength can be expressed as a function of time in
terms of the incomplete elliptic integral of the first kind.
At very early times t � 0, the solution with x ∈ (−π/2, π/2) is asymptotic to two grim
reapers joined smoothly and approaching each other with constant speed:
y = ±(t− log cos(x)− log 2) +O
(
e4t sec4(x)
)
.
A question that may be posed is: why construct curve shortening flows by first solving (10)
rather than solving (2) directly? Remarkably, it can be shown that
Proposition 2. The image of the curve shortening equation (2) under the transformation y 7→
m(y) does not have a (maximal) 3-parameter family of joint solutions with the linear wave
equation yxt = 0 in any field variable m(y) except m = identity in which case the solution gives
rise to the grim-reaper flow.
8 P. Broadbridge and P. Vassiliou
Figure 2. Evolution by heat shrinking flow of the curve defined by cosh y − 5 cosx = 0.
Proof. Construct the image of the curve shortening equation (2) under the change of variable
y 7→ m(y) and apply the Cartan–Kähler theorem to the differential system consisting of the
transformed curve shortening equation and the constraint yxt = 0. �
This is in sharp contrast to the rich separability properties of nonlinear heat equations (8)
as discussed in Section 2. Thus starting with (10) appears to be an important first step in
constructing non-trivial curve flows1.
3.2 Decaying periodic solution initially close to square wave
For a curve fixed at two end-points, we prescribe the boundary conditions
y = 0, x = 0, `.
We now define dimensionless space and time variables X = x/`, Y = y/`, τ = Bt/`2, in terms
of which the curve shortening equation2 is
Yτ =
YXX
1 + Y 2
X
, (15)
to be solved subject to boundary conditions
Y = 0, X = 0, 1 (16)
and continuous initial conditions Y (X, 0) = Y0(X).
Integrating the seventh member of the list (11) and then applying translational and scaling
invariance transformations, we obtain
Y =
1
K
ln
(√
exp(2K2[τ − τ0]) + cos2(K[X −X0]) + cos(K[X −X0])
exp(K2[τ − τ0])
)
, (17)
1In view of the relationship between u(x, t) satisfying (10) and y(x, t) satisfying (2), one might try for the
higher order constraint yxxt = 0 rather than yxt = 0. This is possible but the calculations rapidly become very
complicated as the order increases.
2Also known in the materials science community as the Mullins equation [24].
Symmetry and Separation in Surface and Curve Flows 9
with K, X0, τ0 arbitrary constants. It may be verified that (17) satisfies (15). In fact, when
τ − τ0 is large and negative,
Y =
ln 2
K
−K[τ − τ0] +
1
K
ln cos(K[X −X0]) +O
(
e−2K2[τ0−τ ]
cos(K[X −X0])
)
.
This shows that away from the singularities of ln cos(X −X0), the solution is asymptotic in the
distant past to the Calabi “grim reaper” solution of the curve-shortening flow [17]. At all times,
the solution (17) is a deformation of the grim reaper solution but now it is extended smoothly
and periodically, without singularities over a domain of any length. Unlike in the grim reaper
solution, there are fixed points so that we may apply the Dirichlet boundary conditions (16),
which lead to
Y = ± 1
K
ln
(√
exp(2K2[τ − τ0]) + sin2(KX) + sin(KX)
exp(K2[τ − τ0])
)
, (18)
with
K =
nπ
`
.
The amplitude of Y (X, τ), which is the value of |Y (X, τ)| at X = ±π/(2K), is approximated by
Ymax = K[τ0 − τ ] + (ln 2)/K +O
(
exp
(
−K2|τ − τ0|
))
.
In terms of dimensional quantities,
ymax =
nπB[t0 − t]
`
+ `
ln 2
nπ
+O
(
exp
(
−n2π2B[t0 − t]/`2
))
, (19)
where ymax = `Ymax and t0 = `2τ0/B. This shows that at early times, the amplitude decreases
linearly as a function of time.
As an example, the solution is graphed for the case n = 5 and with the minus sign preceding
the right hand side of (18). The solution at early times is shown in Figs. 3, 4. Although for
large τ0, the initial condition resembles a periodic square wave, it actually converges pointwise
to a differentiable bounded and periodic grim reaper as τ0 approaches ∞.
When K2[τ − τ0] is large and positive, the solution is smoother, approximated by a single
sine wave,
Y = ±exp(−K2[τ − τ0])
K
sin(Kx) +O
(
exp(−2K2[τ − τ0])
)
, (20)
whose amplitude decreases exponentially in time. This is the sinusoidal solution of the classical
linear diffusion equation that approximates the curvature driven diffusion equation at large
times in the small-slope approximation.
The number of extrema n within the fixed domain [0, 1] may be freely chosen. For large
values of the parameter K2τ0, as in Fig. 3, the initial condition is described well as a periodic
square wave, resembling a diffraction grating. The time scale for decay may be viewed as the
time at which the formal expression (19) for Ymax is zero. This is
t = t0 +
`2 ln 2
n2π2B
.
Fig. 3 evidences considerable smoothing of the initial conditions but at this time, the solution
does not yet resemble a simple sinusoid. At larger times, the solution is close to a single sinusoidal
wave. In the regime of the sinusoidal profile, from (20) decay times are shorter, of the order
of the standard time (`/nπ)2/B that is familiar from exponential decay of Fourier modes in
linear diffusion models (e.g. [10]) when the half-wavelength `/n is the typical distance between
neighboring regions of high and low mass concentration.
10 P. Broadbridge and P. Vassiliou
Figure 3. Filled curve: exact solution with τ − τ0 = −0.07, K = 5π. Dashed lines: approximate bounds
y = ±K[τ0 − τ ] + (ln 2)/K.
Figure 4. Filled curve: exact solution with τ − τ0 = −0.02, K = 5π. Dashed lines: approximate bounds
y = ±K[τ0 − τ ] + (ln 2)/K.
4 Surface evolution by evaporation and condensation
For some metals such as gold, surface diffusion persists for several thousand years as the dom-
inant surface transport mechanism but for others such as magnesium, after a time less than
a day, surface evolution occurs predominantly by evaporation-condensation. As described by
Mullins [24], the Gibbs–Thomson relation for evaporation rate in terms of deficit from equi-
librium pressure over a curved surface, leads to a second-order equation for diffusion by mean
curvature. In terms of two-dimensional Cartesian coordinates (x, y) and time t, the Mullins
equation for points on a material surface is
yt = B
yxx
1 + y2
x
,
Symmetry and Separation in Surface and Curve Flows 11
where B is constant. This equation applies to two dimensional cross sections of solids when
surface nano-scale features such as grooves, ridges and furrows extend rectilinearly into the
third dimension. Because of the nonlinearity, very few useful exact solutions to this equation
are known [8], except in a linear approximation. Two decades ago [5], the exact solution was
constructed in parametric integral form for nonlinear surface evolution near a symmetric grain
boundary, with constant slope at the grain boundary groove, initial flatness and zero displace-
ment at infinity. A similar procedure produces the more general “open angle” solutions for
deposition in a wedge [20]. Subsequently, the fourth-order Mullins equation for curvature-driven
surface diffusion on an almost-isotropic material, was solved with boundary conditions repre-
senting a symmetric grain boundary [6, 28, 7].
In the context of surface evaporation, the parallel asymptotes of the grim reaper solution
represent a long thin metallic foil that is evaporating at the ends. The constant travelling-wave
speed of the grim reaper solution shows that evaporation will take place at a constant rate
that depends on foil thickness as well as the evaporation coefficient B. The thickness ∆x is the
distance between the two asymptotes, ∆x = πB/c. Hence the steady evaporation rate at the end
of a strip of metallic foil will be c = πB/∆x. For example, for a foil of a few microns in thickness
made of the volatile metal Mg, this rate will be of the order of one millimetre per millenium.
In the context of metal surface smoothing, the periodic solution of the previous section
predicts the smoothing of initial conditions that resemble a diffraction grating. For a surface
fixed at two end-points, we prescribe the boundary conditions
y = 0, x = 0, `.
Physically, this corresponds to the surface being clamped and shielded from the surrounding
atmosphere outside of the exposed spatial domain [0, `].
4.1 Evaporation from anisotropic crystals
C. Herring [19], showed that for anisotropic crystals, surface energy is proportional to b(φ)(γ +
γ′′(φ)), where φ is the polar angle arctan yx, γ is surface tension and b is an evaporation
coefficient. An equation of the form
yt = D(yx)yxx, (21)
implying the standard nonlinear diffusion equation (8) with u = yx, may be regarded as an
anisotropic form of the Mullins evaporation-condensation equation
yt = B(yx)
yxx
1 + y2
x
,
with surface slope-dependent anisotropy factor B(yx) = D(yx)(1 + y2
x), which originates in the
physical derivation from a constant multiple of b(φ)(γ+ γ′′(φ)) (e.g. [27]). Crystalline materials
are indeed anisotropic, with the evaporation coefficient minimized when the cut surface is aligned
with crystal planes. The group of equivalence transformations of this class of equations includes
the general linear group GL(2,R). By the polar decomposition theorem, each invertible linear
transformation can be decomposed as an othogonal transformation followed by multiplication by
a positive definite symmetric dilatation matrix. Under a rotation about the origin by angle α,
x̄ = x cosα− y sinα, ȳ = x sinα+ y cosα, t̄ = t,
ū = ȳx̄ =
u+ tanα
1− u tanα
, yx =
− sinα+ ȳx̄ cosα
cosα+ ȳx̄ sinα
, yxx =
ȳx̄x̄
(cosα+ ȳx̄ sinα)3
.
12 P. Broadbridge and P. Vassiliou
The axis x̄ = 0 has been rotated by angle −α/2. Writing θ̄ = arctan ȳx̄,
yt =
n · eȳ
n · ey
ȳt̄ =
cos θ̄
cos(θ̄ − α)
ȳt̄ =
1
cosα
1
1 + ȳx̄ tanα
ȳt̄.
Hence, by rotation, the isotropic nonlinear diffusion equation (21) is equivalent to
ȳt̄ = D̄(ȳx̄)ȳx̄x̄, D̄(s) =
1
(s sinα+ cosα)2
D
(
s cosα− sinα
s sinα+ cosα
)
.
Tritscher [26] used this device of rotational equivalence classes to solve integrable forms of
fourth-order surface diffusion equations. In particular, after rotation by angle α = π/2,
D̄(ȳx̄) =
1
ȳ2
x̄
D
(
−1
ȳx̄
)
.
By following the π/2 rotation by a trivial reflection x̄ → −x̄, we recover the result of the
reciprocal transformation (7).
By the principal axis theorem, a positive symmetric matrix can be written as OtQO, where
O is a rotation matrix and Q is diagonal, Qji = aiδ
j
i with aj > 0. Therefore, to consider the
effect of an additional dilatation, we need only consider the effect of a diagonal rescaling:
x̄ = a1x, ȳ = a2y, ȳt = a2
1D
(
a1
a2
ȳx̄
)
ȳx̄x̄,
which although trivial, allows us to construct simple anisotropic models from the rotationally
invariant isotropic model.
5 Anisotropic analogs of isotropic model solutions
The scale invariance group still applies to the general anisotropic diffusion equation (21). There-
fore both expanding and shrinking types of similarity solution exist but they may be significantly
different from those of the isotropic model. The grain boundary groove solution still exists, as
can be seen from the solvability of the reduced boundary value problem on R+ × R+,
u = F (ρ), ρ = xt−1/2,
−ρ
2
F ′(ρ) =
d
dρ
[D(F )F ′(ρ)],
F (0) = m > 0, F (ρ)→ 0, ρ→∞.
The grain-boundary groove solutions for all models have some common features. Although it
represents a highly anisotropic material, the linear model groove solution approximates that of
an isotropic model for groove slopes of up to 0.5, which was used in Mullins’ original paper [24].
However, whereas the linear model predicts that the groove depth increases in proportion to m,
the groove depth increases very slowly, of order (logm)0.5 at large m for the isotropic model [3].
5.1 Anisotropic homothetically shrinking closed curve
For anisotropic evaporating materials, the closed-curve homothetic solution represents the fixed-
shape cross section of an evaporating wire, which is circular when the material is isotropic.
The homothetic evaporating closed-curve solution satisfies
y = [t0 − t]1/2G(ρ), ρ = x[t0 − t]−1/2,
Symmetry and Separation in Surface and Curve Flows 13
−G
2
+
ρ
2
G′(ρ) = G′′(ρ)D(G′),
G′(0) = 0, G(ρ0) = 0, G′(ρ)→∞, ρ→ ρ0.
This implies
u = F (ρ), F = G′,
ρ
2
F ′(ρ)−D′(F )(F ′)2 = F ′′(ρ)D(F ), (22)
F (0) = 0, F (ρ)→∞, ρ→ ρ0.
For some functions D(u), the homothetic closed-curve solution does not exist. For example,
with the linear model with constant D, the general solution satisfying F (0) = 0 must be
F (ρ) =
√
π
2
A erf (ρ) = A
∫ ρ
0
e−s
2/4ds, A ∈ R.
Since this cannot take an infinite value at any point ρ = ρ0 ∈ R, the closed-curve homothetic
solution does not exist.
Let us now assume that F (ρ) ≈ A0(ρ0 − ρ)ν . If a homothetic closed-curve solution exists,
G(ρ0) ∈ (−∞, 0), and
0 = G(ρ0) = G(0) +
∫ ρ0
0
F (s)ds.
This can hold only if the integral does not diverge, implying
−1 < ν < 0.
Now we suppose that at large u,D(u) is asymptotic to a power law D(u) ≈ D0u
n. Then by
balancing terms at the leading order in ρ0 − ρ, (22) implies
(i) n = 1/ν, which is less than −1, and ρ0 = −2νD0A
1/ν
0 , or,
(ii) n = 1−ν
ν , which is less than −2.
For the linear model, the anisotropy factor is B(u) = (1 + u2)D which diverges when the
curve is vertical. In the current application, we are interested in cases of realistic anisotropy for
which the evaporation coefficient and the anisotropy factor are bounded, and the latter with
a minimum value greater than zero:
∀u ∈ R, 0 < B0 < B(u) < B∞ <∞.
That statement must be true, independent of the orientation of the coordinate axes. For example
after rotation by π/2 it must be true that D̄(ū) → D̄0 ∈ (0,∞) as ū → 0. By rotating back to
the original orientation, this implies
D(u) = u−2D̄(u) ∼ D̄0u
−2.
This physical restriction rules out scenario (ii) in the above leading-order analysis. The simplest
example satisfying the restriction is D(u) = 1/(1 + (βu)2) (β ∈ R constant) for which the
anisotropy factor B(yx) varies between 1 and β2. In this case the homothetic solution is simply
an ellipse
x2
2[t0 − t]
+
y2
2β2[t0 − t]
= 1
elongated in the direction of weakest evaporation.
14 P. Broadbridge and P. Vassiliou
5.2 Grim reaper solution for anisotropic material
The grim reaper solution is a travelling wave constrained between two vertical asymptotes. This
implies a steady state solution of the equation ut = [D(u)ux]x. For arbitrary D(u), there exists
a two-parameter steady state solution
us(x) = K−1(cx+ c2),
where
K(u) =
∫ u
0
D(s)ds
which is an increasing invertible function because D(s) > 0. This integrates formally to
y(x, t) = y0(t) +
∫ x
0
us(x1, t)dx1.
It then follows from (21) that y0(t) can only be a constant-velocity translation, y0 = c(t − t0).
This travelling wave solution does not necessarily have a vertical asymptote. For example, there
is no such asymptote when D is constant. However if D(s) satisfies the physical requirement
D(s) ≈ s−2, it follows that K(u) has a finite limit as u → ∞, therefore K−1 has a vertical
asymptote at some location x = x0, duplicated at x = −x0 if D(s) is an even function, as is
commonly the case when the x-axis denotes the orientation of the crystal planes from where
evaporation is weakest.
6 Anisotropic models allowing functional separation
The classification of Doyle and Vassiliou gives all functions D(s) for which functional separation
of variables is possible in the general form
ū = f(u) = v(x) + w(t), (23)
with f an invertible function.
The simplest example of a physically feasible isotropic model is simply that obtained from
the isotropic diffusion equation by unequally rescaling x and y. For example, Fig. 2 could be
dilated in one direction, displaying a non-homothetic closed curve approaching a homothetically
shrinking ellipse.
The only other member of the Doyle–Vassiliou list with realistic anisotropy is much more
complicated:
D(u) = D0 cos(z(Au)), Au =
∫ z
0
(cos s)−3/2ds, −π/2 < z < π/2. (24)
This model is close to isotropic when A =
√
2. In that case, it is easy to show that when
u(= yx) is small, B(u) = 1 +O(u4) and B → 2 as u→∞. In Fig. 5, the function D(u)/D0 for
this weakly anisotropic model is compared to D(u)/D0 of the isotropic model.
For the sake of completeness, we construct the special solution u(x, t) compatible with (23),
that was not given explicitly by Doyle and Vassiliou [14]. The parameter A may be changed by
rescaling u. For convenience, without loss of generality we now set A to 1. Also we may set D0
to 1 by using D0t as the time coordinate. From the general approach of Doyle and Vassiliou [14],
ū is a sum of separated functions v(x) and w(t) satisfying
v′′′(x) +
1
2
[v′(x)]3 = 0, (25)
Symmetry and Separation in Surface and Curve Flows 15
Figure 5. Nonlinear diffusivity function D(u) for isotropic evaporation (lower) and anisotropic model
(upper).
w′(t) = W (w) = cosw
[
v′′ cos v +
1
2
(v′)2 sin v
]
+ sinw
[
−v′′ sin v +
1
2
(v′)2 cos v
]
.
By construction, the following are first integrals of (25), which can be verified by substitution:
λ = v′′(x) cos v +
1
2
(v′)2 sin v, µ = −v′′(x) sin v +
1
2
(v′)2 cos v.
It then follows that w′(t) = λ cosw + µ sinw, and that
t =
2
R2
∫
dw
sin(w + δ)
, x =
1
R
∫ v
0
ds√
cos(s− δ)
, (26)
with constants x0, t0, R = ±21/2(µ2 + λ2)1/4 and δ = arctan(λ/µ), whose arbitrariness is of
little consequence. Let v′(x) = W (v), so v′′(x) = WW ′(v). Then
WW ′(v) cos v +
1
2
W 2 sin v = λ, −WW ′(v) sin v +
1
2
W 2 cos v = µ
⇒ −λ tan v +
1
2
W 2 sin2 v
cos v
+
1
2
W 2 cos v = µ
⇒ W 2 = 2(µ cos v + λ sin v) = R2 cos(v − δ) ⇒ v′(x) = ±R
√
cos(v − δ),
x =
1
R
∫ v
0
ds√
cos(s− δ)
.
When we include the parameter B0, the general solution for v and w in terms of elementary
functions, the standard elliptic integral F (θ|12) and the standard Jacobi elliptic function sn(X|12),
is
w = 2 arctan eR
2B0t/2 − δ and v = 2 arcsin
(
1√
2
sn
(
R[x− x0]√
2
∣∣∣1
2
))
+ δ,
where
x0 =
√
2
R
F
(
arcsin
(√
2 sin
δ
2
) ∣∣∣1
2
)
.
16 P. Broadbridge and P. Vassiliou
Figure 6. Portion of solution curve for anisotropic model with R = 1 and δ = 0, at t = 0.
Figure 7. Portion of solution curve for anisotropic model at t = −2.0.
Since u is a function of w(t) + v(x), we see that the choice of parameters R and δ has little
consequence on the form of the solution u(x, t). A shift in the phase variable δ has the effect
of a translation in x by x0(δ,R). A change of amplitude R has the effect of rescaling x and t
to Rx and R2t.
Modulo a time-dependent vertical translation, the curve y(x, t) is obtained from u(x, t) by
integration. Since the integrands in (24) and (26) must be real valued, the constructed solution
has a truncated domain. An example is given in Figs. 6–8.
As t approaches −∞, the solution approaches a steady state. Since D(u) is symmetric, this
steady state is a symmetric grim reaper. It is approximated in Fig. 8 by taking t = −8.0.
The curve is not symmetric, as can be seen in Fig. 7. This asymmetry causes its domain to
shift slightly from right to left. For t ≤ 0, the curve has a vertical asymptote at moving location
x = xu =
∫ π/2−w(t)
0
(cos s)−1/2ds
Symmetry and Separation in Surface and Curve Flows 17
Figure 8. Portion of solution curve for anisotropic model at t = −8.0.
and a local minimum at
x = x` = −
∫ w(t)
0
(cos s)−1/2ds.
7 Conclusion
Of the exact solutions to the curve shortening equation known to us, most can be obtained by
Lie point symmetry reductions. The two interesting solutions that cannot be constructed in this
way, can indeed be recovered by functional separation of variables for the standard nonlinear
diffusion equation (3) that is obtained from the curve shortening equation by differentiation.
The classification obtained by Doyle and Vassiliou [14] of nonlinear diffusion equations that
admit functional separation of variables, leads to the two exact non-self-similar solutions with
non-trivial initial conditions that appear to be achievable in an approximate sense in applica-
tions. In addition, it leads to a new separated solution for a physically realistic anisotropic
evaporation-condensation diffusion equation. Although the second order nonlinear surface evo-
lution equations for slope u(x, t) (≡ yx(x, t) admit a number of possibilities for functional se-
paration of variables in Cartesian coordinates, we have proved that this is not possible for the
equation (2) governing y(x, t), nor is it possible in a coordinate system consisting of canonical
variables for a symmetry other than translation. This is in contrast with the point symmetry
analysis, which leads to a richer array of possibilities for the evolution of y(x, t) than for the
evolution of u(x, t).
The invariance of the isotropic equation (10) under the well known reciprocal transformation
was shown (Proposition 1) to lead to no new planar curve heat flows. The group of geometric
equivalence transformations of the class of general anisotropic equations (21) includes not only
the reciprocal transformation in the guise of a reflection in the plane, but the whole general linear
group. The equivalence is shown by carrying out the equivalence transformations explicitly.
Physical evaporation coefficients must have a positive real value when the surface is oriented
along the crystal planes. Since physical restrictions must be independent of orientation of the
coordinate axes, it follows from the equivalence transformations that the nonlinear diffusivity
D(yx) must behave like y−2
x at large-slope. This also allows for the existence of a closed non-
circular homothetic solution which cannot exist unless D(u) decreases faster than u−1. The
Doyle–Vassiliou classification produces another anisotropic model that satisfies this physical
18 P. Broadbridge and P. Vassiliou
requirement. An exact solution has been constructed, involving Jacobi elliptic functions and
other inverse integrals of rational functions.
Exact solutions sometimes have the advantage of leading to concise conceptually simple
relationships. For example, Fig. 3 demonstrates the efficacy of a simple expression for wave
amplitude of a corrugated nano-scale surface in the early stages of smoothing by evaporation-
condensation when the system cannot be adequately described by a linear model. However,
exact solutions can be obtained only in very special cases of initial and boundary conditions, so
that approximate numerical solution methods will continue to be important.
Acknowledgements
This paper is submitted in appreciation of the valuable on-going contributions of Professor
Willard Miller Jr. The first author gratefully acknowledges support by the Australian Research
Council under project DP1095044.
References
[1] Abresch U., Langer J. The normalized curve shortening flow and homothetic solutions, J. Differential Geom.
23 (1986), 175–196.
[2] Angenent S., On the formation of singularities in the curve shortening flow, J. Differential Geom. 33 (1991),
601–633.
[3] Arrigo D.J., Broadbridge P., Tritscher P., Karciga Y., The depth of a steep evaporating grain boundary
groove: application of comparison theorems, Math. Comput. Modelling 25 (1997), no. 10, 1–8.
[4] Bluman G.W., Reid G.J., Kumei S., New classes of symmetries for partial differential equations, J. Math.
Phys. 29 (1988), 806–811, Erratum, J. Math. Phys. 29 (1988), 2320.
[5] Broadbridge P., Exact solvability of the Mullins nonlinear diffusion model of groove development, J. Math.
Phys. 30 (1989), 1648–1651.
[6] Broadbridge P., Tritscher P., An integrable fourth-order nonlinear evolution equation applied to thermal
grooving of metal surfaces, IMA J. Appl. Math. 53 (1994), 249–265.
[7] Broadbridge P., Goard J.M., Temperature-dependent surface diffusion near a grain boundary, J. Engrg.
Math. 66 (2010), 87–102.
[8] Cahn J.W., Taylor J.E., Overview no. 113 surface motion by surface diffusion, Acta Metall. Mater. 42
(1994), 1045–1063.
[9] Cao F., Geometric curve evolution and image processing, Lecture Notes in Mathematics, Vol. 1805, Springer-
Verlag, Berlin, 2003.
[10] Carslaw H.S., Jaeger J.C. Conduction of heat in solids, 2nd. ed., Clarendon Press, Oxford, 1959.
[11] Chou K.-S., Li G.-X., Optimal systems and invariant solutions for the curve shortening problem, Comm.
Anal. Geom. 10 (2002), 241–274.
[12] Clarkson P.A., Fokas A.S., Ablowitz M.J., Hodograph transformations of linearizable partial differential
equations, SIAM J. Appl. Math. 49 (1989), 1188–1209.
[13] Daskalopoulos P., Hamilton R., Sesum N., Classification of compact ancient solutions to the curve shortening
flow, J. Differential Geom. 84 (2010), 455–465, arXiv:0806.1757.
[14] Doyle P.W., Vassiliou P.J., Separation of variables for the 1-dimensional non-linear diffusion equation,
Internat. J. Non-Linear Mech. 33 (1998), 315–326.
[15] Galaktionov V.A., Dorodnitsyn V.A., Elenin G.G., Kurdyumov S.P., Samarskii A.A., A quasilinear equation
of heat conduction with a source: peaking, localization, symmetry, exact solutions, asymptotic behavior,
structures, J. Soviet Math. 41 (1988), 1222–1292.
[16] Gage M., Hamilton R., The heat equation shrinking convex plane curves, J. Differential Geom. 23 (1986),
69–96.
[17] Grayson M. 1987 The heat equation shrinks embedded plane curves to round points, J. Differential Geom.
26 (1987), 285–314.
http://dx.doi.org/10.1016/S0895-7177(97)00070-8
http://dx.doi.org/10.1063/1.527974
http://dx.doi.org/10.1063/1.527974
http://dx.doi.org/10.1063/1.528113
http://dx.doi.org/10.1063/1.528300
http://dx.doi.org/10.1063/1.528300
http://dx.doi.org/10.1093/imamat/53.3.249
http://dx.doi.org/10.1007/s10665-009-9332-9
http://dx.doi.org/10.1007/s10665-009-9332-9
http://dx.doi.org/10.1016/0956-7151(94)90123-6
http://dx.doi.org/10.1137/0149071
http://arxiv.org/abs/0806.1757
http://dx.doi.org/10.1016/S0020-7462(97)00013-9
Symmetry and Separation in Surface and Curve Flows 19
[18] Halldorsson H.P., Self-similar solutions to the curve shortening flow, arXiv:1007.1617.
[19] Herring C., Surface tension as a motivation for sintering, in The Physics of Powder Metallurgy, Editor
W.E. Kingston, McGraw-Hill, 1951, 143–179.
[20] Ishimura N., Curvature evolution of plane curves with prescribed opening angle, Bull. Austral. Math. Soc.
52 (1995), 287–296.
[21] King J.R., Emerging areas of mathematical modelling, R. Soc. Lond. Philos. Trans. Ser. A Math. Phys.
Eng. Sci. 358 (2000), no. 1765, 3–19.
[22] Kingston J.G., Rogers C., Reciprocal Bäcklund transformations of conservation laws, Phys. Lett. A 92
(1982), 261–264.
[23] Malladi R., Sethian J.A., Image processing via level set curvature flow, Proc. Nat. Acad. Sci. USA 92 (1995),
7046–7050.
[24] Mullins W.W., Theory of thermal grooving, J. Appl. Phys. 28 (1957), 333–339.
[25] Olver P.J., Sapiro G., Tannenbaum A., Invariant geometric evolutions of surfaces and volumetric smoothing,
SIAM J. Appl. Math. 57 (1997), 176–194.
[26] Tritscher P., Integrable nonlinear evolution equations applied to solidification and surface redistribution,
PhD Thesis, University of Wollongong, 1995.
[27] Tritscher P., An integrable fourth-order nonlinear evolution equation applied to surface redistribution due
to capillarity, J. Austral. Math. Soc. Ser. B 38 (1997), 518–541.
[28] Tritscher P., Broadbridge P., Grain boundary grooving by surface diffusion: an analytic nonlinear model
for a symmetric groove, Proc. Roy. Soc. Lond. A 450 (1995), no. 1940, 569–587.
http://arxiv.org/abs/1007.1617
http://dx.doi.org/10.1017/S0004972700014714
http://dx.doi.org/10.1098/rsta.2000.0516
http://dx.doi.org/10.1098/rsta.2000.0516
http://dx.doi.org/10.1016/0375-9601(82)90081-0
http://dx.doi.org/10.1063/1.1722742
http://dx.doi.org/10.1137/S0036139994266311
http://dx.doi.org/10.1017/S0334270000000849
http://dx.doi.org/10.1098/rspa.1995.0101
1 Introduction
2 Symmetry reductions of the curve shortening equation
2.1 Reciprocal transformations
3 Functionally separable nonlinear heat equations
3.1 Exact heat flow of a convex curve
3.2 Decaying periodic solution initially close to square wave
4 Surface evolution by evaporation and condensation
4.1 Evaporation from anisotropic crystals
5 Anisotropic analogs of isotropic model solutions
5.1 Anisotropic homothetically shrinking closed curve
5.2 Grim reaper solution for anisotropic material
6 Anisotropic models allowing functional separation
7 Conclusion
References
|