The fundamental solutions of the curve shortening problem via the Schwarz function

Curve shortening in the z-plane in which, at a given point on the curve, the normal velocity of the curve is equal to the curvature, is shown to satisfy StSz=Szz\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S_tS_z=S_{zz}$$\end{document}, where S(z, t) is the Schwarz function of the curve. This equation is shown to have a parametric solution from which the known explicit solutions for curve shortening flow; the circle, grim reaper, paperclip and hairclip, can be recovered.


Introduction
The evolution of a smooth curve γ γ γ (x, t) in the plane with given initial shape γ γ γ 0 in which the normal velocity at a given point on the curve is proportional to the curvature κ at that point, is known as the curve shortening problem: where n is the normal. The problem, and its higher dimensional generalisation, has attracted much attention since the 1980s. Important features of the curve evolution, such as the shrinking of any closed, embedded curve first to a convex curve and then to a round point in finite time, have been established e.g. [5,6]. Curve shortening (1), and its variants, also has practical application; for example in the late-time evolution of Hele-Shaw free boundary flow in the presence of surface tension [3]. In the time reversed sense when the curve lengthens, the solutions have relevance to viscous fingering and crystal growth e.g. [12] where the connection to Saffman-Taylor fingering is made. There are only four known, explicit, time-dependent solutions describing curve shortening: the bounded, closed B N. R. McDonald n.r.mcdonald@ucl.ac.uk 1 Department of Mathematics, University College London, London WC1E 6BT, UK curve solutions of the (i) circle and (ii) paperclip, and the unbounded, (iii) steadily translating grim reaper, and the (iv) hairclip solution-see e.g. [1,12,15]. A straight line is a trivial solution and remains stationary. More exotic, but not explicit, solutions taking the form of spirals which rotate and/or expand are described in [8].
There is interest in deriving these exact solutions systematically using simple methods. One approach has been to show the fundamental solutions can be obtained from a nonlinear diffusion equation [1,12] which can be solved by functional separation of variables in conjunction with symmetry assumptions. More recently, Reference [15] presents two methods for deriving the four fundamental solutions. The first is based on seeking explicit solutions to (1) in which the function being sought is itself a function satisfying the classic heat equation h t = h x x . Assuming h t is either constant or linear in h gives each of the four fundamental solutions. The other approach detailed in [15] is an implicit method in which the solution is sought in the form of separable functions of the three variables x, y, t. An alternative and systematic procedure for obtaining the exact solutions for time-dependent curve shortening is presented here. The key idea is to represent the curve in the complex plane using the Schwarz function [4]. Then, in Sect. 2, well-established results for representing the geometric properties of the curve in terms of the Schwarz function are used to derive a PDE governing the evolution of the curve. Section 3.1 shows that the self-similar, collapsing circle is an exact solution of the PDE, and Sect. 3.2 uses elementary methods to derive a steadily translating parametric solution of the PDE which is equivalent to the grim reaper solution. In Sects. 3.3 and 3.4, the parametric solution is generalised and shown to yield the paperclip and hairclip solutions. Section 4 reconsiders the steadily translating problem, and derives directly the grim reaper solution from the Schwarz function equation cast in a complex plane moving with the solution curve. The same approach is used to find a new exact solution taking the form of a steadily translating parabola to the problem in which the normal velocity of the curve is given by the cube root of the curvature.

The governing equation in terms of the Schwarz function
The Schwarz function S(z) of an analytic curve γ γ γ in the z-plane is the unique function which is analytic in the neighbourhood of γ γ γ and such that S(z) =z on γ γ γ [4]. For example, a circle of radius a, centred at the origin, has S(z) = a 2 /z. Representing curves using S(z) has proved remarkably successful in obtaining exact solutions in fluid mechanical applications where the evolving curve may represent, for example, the interface between fluids of different viscosities such as occurs in a Hele-Shaw cell (e.g. [7,9,11]), or, in 2D vortex dynamics, where the interface separates flows with different constant vorticities (e.g. [2,13,14]). The common theme in these applications is that the dynamics of a planar curve can be cast in terms of an equation satisfied by the Schwarz function which, when solved, determines the shape of the curve. Often the success of this approach is owing to the fact that S(z) is an analytic function, at least in the neighbourhood of γ γ γ , enabling methods of complex analysis, such as conformal mapping, to be used.
where subscripts denote partial derivatives. For an evolving curve S = S(z, t), the normal velocity v n at point z on γ γ γ is (e.g. [11]) Curve shortening equates (2) and (3), giving a nonlinear PDE for S(z, t): The curve shortening problem involves solving (4) subject to the initial shape of the curve S(z, 0). An immediate feature of (4) is, givenz = S(z), that it is invariant to the scaling z → λz (λ ∈ ), which implies z = S → λS, and t → λ 2 t-see also [15]. Moreover, as expected, the governing equation (4) is invariant to rotations through angle σ , since z → e iσ z implies S(z) → e −iσ S, rendering (4) invariant.
Letting s be the arclength parameter of γ γ γ , it can be shown (e.g. [4]) dz/ds = 1/ √ S z , and upon conjugating dS/ds = √ S z . These relations can be used in (4) to establish the known connection between the curve shortening flow and the heat equation (e.g. [15]), which in term of the Schwarz function is 3 Explicit solutions of the curve shortening flow

The collapsing circle
By symmetry, an initially circular curve remains circular under the action of curve shortening: this is a known exact (self-similar) solution. It is simple to demonstrate using (4): a circle of radius a(t) with (without loss of generality) centre at the origin, has Schwarz function S(z, t) = a(t) 2 /z which satisfies (4) exactly provided aȧ = −1, whereȧ is the time derivative of a(t). This is equivalent toȦ = −2π where A = πa 2 is the area enclosed by the circle, and consistent with the well-known result that the area enclosed by an arbitrary closed plane curve decreases at the rate of 2π per unit time e.g. [6].

The grim reaper
Let a curve translate steadily, with unit velocity, in the positive real direction, and write it as where ζ = exp (iθ), −π ≤ θ ≤ π , andf = f . The latter condition implies the curve is symmetric about the real axis. Now It follows from (6) and (7) that Also, Substituting (8) and (9) into (4) and simplifying gives Let h(θ ) = e iθ f (e iθ ) and using ∂ ζ = ie iθ ∂ θ , (10) gives a functional-differential equation for the complex-valued h(θ ): Sincef = f , it follows that h(−θ) = h(θ ) and that the real part of (11) is satisfied by inspection with both sides of the equation vanishing.
Letting h(θ ) = r (θ ) exp(iα(θ)) and substituting into (11) gives But and is the general solution to (11). Now h(θ ) = e iθ f (e iθ ) implies d f /dθ = ih and integrating gives where a constant of integration (which is real sincef = f ) has been set to zero since it only serves to shift the curve along the real axis.
Equations (6) and (13) gives and taking real and imaginary parts of (14) gives x = t − log(cos y) which is the grim reaper solution. Note that (14) does not admit embedded curve solutions; the curve must close at infinity. That is, α(±π) = ±π/2. The choice of function α(θ) is immaterial so long as it is a monotonic function such that −π/2 ≤ α(θ) ≤ π/2. Choosing the simple form α(θ) = θ/2 gives the final parametric form of the grim reaper solution

The paperclip
The grim reaper solution in parametric form (15) suggests a general form of solution to (4): where a(t) > 0 (without loss of generality) and g(t) are realvalued functions to be determined. First, the choice a(t) > 1 is made and it is shown this corresponds to the paperclip solution. In Sect. 3.4, the choice 0 < a(t) < 1 is shown to correspond to the hairclip solution.

The hairclip
In this section, it is shown that the hairclip solution corresponds to the choice 0 < a(t) < 1 in the parametric representation (16) with g(t) chosen to ensure that the solution curve is centred about the the imaginary axis, giving Proceeding as in Sect. 3.3 and substituting (27) into (4) yields an ordinary differential equation for a(t) with solution where an unimportant choice of arbitrary integration constant has been made. Note that −∞ < t < ∞ and the solution (27) and (28) can be regarded as an eternal solution of the curve shortening problem. Again, by considering real and imaginary parts of (27) and following similar steps as in Sect. (27) and (28) is equivalent to the hairclip solution e −t cos y = − sinh x.

Revisiting the grim reaper solution
An alternative solution procedure based on solving (4) in a steadily moving frame of reference gives the grim reaper solution in the explicit formz = S(z, t). Let the Z -frame move in the positive real direction with unit speed so that In terms of f (Z ), (4) becomes the ODE which has general solution where K and C are constants. Without loss of generality the condition f (0) = 0 is imposed, which implies the solution curve passes through Z = 0, and gives C = − log(1 + K ). Further, note that as z → ∞, f → Z + log(K /(1 + K )), and so in this limit the solution curve behaves like y → (i/2) log(K /(1 + K )). Thus |K /(1 + K )| = 1. Demanding that the solution is symmetric about the real Z -axis implies that K /(1 + K ) = −1 (i.e. y → ±π/2 as z → ∞) and so and the explicit Schwarz function representation of the grim reaper solution is By taking real and imaginary parts of (33) it is straightforward to verify that it is equivalent to the standard grim reaper expression x = t − log(cos y).
While it is known that the grim reaper is the only steadily translating solution to the curve shortening problem (e.g. [8]), the demonstration in this section suggests that the search for steadily translating solutions to more general curve evolution problems in which v n = F(κ), where F is a differentiable function, might be fruitfully pursued by this approach. That is, by formulating the ODE version of v n = F(κ) in a moving frame, and solving the ODE to find explicit solutions of the Schwarz function.
To give an example of this suppose F(κ) = κ 1/3 , so that the curve evolution v n = κ 1/3 gives the Schwarz function PDE Seeking a solution to (34) translating with unit speed in the positive real direction implies that S(z, t) = f (Z )+t where, from (34), where Z = z − t. The general solution of (35) is f = Z + √ 4 A − 8Z + B, where A and B are constants. Since f − Z = −2iIm(Z ), squaring and considering the imaginary part gives B = −2. Further, requiring the solution curve passes through Z = 0 gives A = 1. Thus, in the z-plane the solution curve has Schwarz function which upon taking real and imaginary parts implies the curve is the translating parabola y 2 = 2(x − t). It is straightforward to verify that this translating parabolic solution γ γ γ satisfies (34) by directly calculating its curvature and showing κ 1/3 = 1/ 1 + y 2 , and showing this is the same as the normal velocity v n = (∂γ γ γ /∂t).n.

Remarks
The Schwarz function formulation of the curve shortening problem in the complex plane, Eq. (4), together with assumptions on the symmetry of the solutions, enables a relatively simple derivation of the fundamental solutions of curve shortening. The four solutions can be summarised by the relation where ζ = exp(iθ) and with the following realisations: 1. a 2 = 1/(1 − exp(2t)), −∞ < t < 0, −π ≤ θ ≤ π is the paperclip solution; 2. a 2 = 1/(1 + exp(2t)), −∞ < t, θ < ∞ is the hairclip solution; 3. The limit t → −∞ of (37) gives a 2 → 1 + exp(2t) and recovers the grim reaper solution (15). That is, the ancient time limit of the hairclip or paperclip consists of either an ensemble, or pair, of grim reapers which approach each other as t increases forming the hair-and paperclip solutions respectively; 4. The limit t → 0 of the paperclip solution, gives a → ∞ and z → ζ /a; that is, a circle with vanishingly small radius.
Note that solutions of the form (37) with |a(t)| ≤ 1 also arise in the analysis of finger evolution resulting from the instability of an interface separating fluids with different viscosities in the absence of surface tension e.g. [9,10], where superposition of such solutions was used to examine their stability and to construct new solutions for evolving fingers. Owing to the nonlinearity of (4) considering a superposition of solutions does not seem to be a way of generating further exact solutions, but it may offer a useful approach to the numerical study of curve shortening of, say, periodic interfaces.
Finally, as noted in Sect. 2, curve shortening flow satisfies, in terms of the Schwarz function, the heat equation S t = 2S ss (5). This form of the evolution equation is strongly connected to the geometric evolution equation used by, e.g. [15], to find the fundamental curve shortening solutions. It would be of interest to further pursue the connection between the Schwarz function approach of the present work to geometric heat equation-based methods.