The asymptotic behavior of solutions to an area-preserving curvature flow of planar curves in the upper half plane is investigated. Two endpoints of the curve slide along the horizontal axis with prescribed fixed contact angles. First, by establishing an isoperimetric inequality, we prove the global existence of the solution. We then study the asymptotic behavior of solutions with concave initial data near a traveling wave.
In this paper, we deal with the motion by area preserving curvature flow of planar curves having two moving end points on the x-axis with fixed interior contact angles to this axis. The problem is formulated as follows. For any curve γ lying on the upper half-plane having two end points, we shall denote the left and right interior angles by and . For any given such curve , our problem is to look for a family of curves having two end points on the x-axis with the same fixed contact angles as above evolve by the area preserving curvature flow equation:
Here V is the outward normal velocity of the curve, is the length of and κ is the (signed) curvature of . If the solution curve is represented as a graph, the problem is reduced to the following free boundary problem (P):
where and we assume that , . Note that the arclength parameter s is given by . The purpose of this paper is to prove that the curve remains uniformly bounded and to study the asymptotic behavior of the solution as for the problem (P). Throughout this paper, we also assume and
Note that these conditions also yield the positivity on .
As is easily seen, a direct calculation yields
where (the area enclosed by the graph and x-axis). Here we used the identities
which follow from differentiating by t and using the boundary condition (1.4). For this reason, Eq. (1.1) is so-called area-preserving curvature flow, which arises in application fields [4,8,9,14,16] such as phase transitions, image processing and droplet model, etc. For the curve shortening flow , there is extensive literature on this subject both for simple closed curves and non-simple ones see [2,10,11] for example. As for the above free boundary problem for , we refer to the work [5] and the references cited therein. Recently, the authors of [12] consider the flow with a constant under the same boundary condition. They first concern the classification of solutions. There are three possibility for the behavior of the solution, that is, area expanding, area bounded and area shrinking types. The detail asymptotic behavior of solutions and eventual concavity property are considered. They show that solutions are asymptotically self-similar both in the area expanding and the area shrinking cases. For the area bounded case, solutions converge to either a stationary solution or a traveling wave. Our main result explained in the below is similar to the convergence result for the area bounded case of [12], but the proof in our paper is completely different. The reason is that the intersection argument of [12], which was exceedingly useful, cannot be applied to our problem (P).
Our first main result is related to global existence and boundedness of concave solutions:
(Globally bounded).
Any concave solution of (P) exists globally in time. The solutionand all derivatives are uniformly bounded. Furthermore,remains bounded from above and below by two positive constants as.
To prove this theorem, we first need to show that the length of the curve is uniformly bounded. In order to obtain it, we show that the aspect ratio of the curve remains bounded by using an idea similar to Proposition 6.2 of [12]. Combining this and the local existence theorem, we can prove the solution exists globally in time.
The second main result of this paper is related to the asymptotic behavior of the solution. A traveling wave of the problem (P) is a solution that has the form , where c denotes the wave speed, is called the profile of the wave and a is an arbitrary constant that adjusts the phase. Substituting this form into (1.2)–(1.4), we obtain
where . Here is some constant that shift the center of the support of to the origin. We shall discuss the existence of traveling wave solutions in Section 4. Roughly speaking, for any given area parameter (the area enclosed by the traveling wave and x-axis), there exists concave traveling wave solution and that is unique up to translation. In order to find a traveling wave solution, we analyze the corresponding evolution equation for the curvature, whose stationary solution corresponds to a traveling wave by following the idea of [12].
The second main theorem of this paper is the following.
(Asymptotic shape).
Letbe a solution of (P) with (1.6) and letbe a solution of (1.9) with. Letbe the curvature of the solution curve, where. Letwith. Then there exists a positive constantsuch that for any solution u satisfying the condition, there exists a constantsuch thatasanduniformly on. Here we regard thatoutside the intervalandoutside the interval.
The argument of the proof is based on the spectral analysis of the equation for the curvature, but the linearized operator of the equation for the curvature around the stationary solution is difficult to handle, because of the non-local term arising from the area preserving property. Then we recover the original variable of (P) from κ to conclude the uniform convergence.
The rest of this paper is organized as follows. In Section 2, we review from [12] the local existence and uniqueness of the solution to the problem (P) and also the continuous dependence of solutions on the initial data. In Section 3, we prove Theorem 1.1, by showing the preservation of the concavity and the boundedness of the length. We also provide some derivative estimates. In Section 4, we analyze the stationary problem of the equation for the curvature. In Section 5, we study more detailed asymptotic and prove Theorem 1.2.
Short time existence
The authors of [12] consider the local existence of the solution for , where is a constant. They first rewrite the evolution equation to the angle variable and use the normalized arclength parameter as the space variable. The problem will be converted into a simple semilinear problem, which makes it possible to apply a well-documented general theory for semilinear problems. They assume with for technical reasons, where is the little Hölder space. This makes the assumption slightly stronger than that in [3], but their approach is more elementary and more self-contained; furthermore, the continuous dependence result is considered. As before, we introduce
where , and I is a closed interval. We shall write the statement of the local existence theorem as [12].
Let the initial data satisfyfor somewith. Assume also thatfor some constant. Then there existanddepending only on N such that the problem (P) has a unique classical solutionwhere, and that this solution satisfiesFurthermore,isinandisin, and for anyand, there exists a constantdepending only on N, k andsuch that
Now we shall outline the strategy of the proof, which is similar to that of Section 7 in [12]. We define a function by . Then Θ satisfies
for along with the following free boundary condition:
Next, we normalize the interval to a fixed interval by using the transformation
Then the above problem is reduced to the following problem for , :
where
The initial conditions for (2.3) are given by
where is the inverse function of
By (2.2), we have
This implies that, once the solution of (2.3) is obtained, one can recover the original solution of (P) as follows:
where and are defined through the relations:
Note that determine first, and is determined later from these ones. We apply the standard analytic semigroup theory to (2.3) in the same setting as in [12] and, by using (2.4), we can prove Theorem 2.1. We can also check that the following theorem holds from the proof of Theorem 7.2 of [12] easily. Here is the maximal existence time of the solution to the problem (P).
Let the functionsandsatisfyfor all. Suppose thatwhere we setoutside ofand, respectively. Furthermore, we assume that there exists a constantsuch thatThen for any, it holds thatfor sufficiently large n and that
The problem (2.3) is well-posed in for as in [12]. This implies that the problem (4.1) is well-posed . If two initial datum of the problem (4.1) in Section 4 are sufficiently close uniformly, then the corresponding solutions become sufficiently close in for sufficiently small , by the smoothing effect. This also holds for the generalized problem (5.1) in Section 5.
Some a priori estimates
In this section, we provide some preliminaries for later purposes. By Theorem 2.1, the solution is in x, t for . Hence we are able to calculate higher derivatives of in the following proof. Let be the maximal time interval for the existence of a solution u to (P). In this subsection, we shall use the idea of [12] (or [5]) to derive some a priori bounds for , and , respectively. In the following, we shall introduce the following notations:
The following estimate holds:Furthermore
Let v be a solution of (2.3). By the concavity of the initial data
It follows from the strong maximum principle that in . Thus (3.1) follows. Note that and . Thus, by applying the Hopf lemma, we get
(Preservation of concavity).
Any solution of the problem (P) remains strictly concave for all.
Let , where μ is a constant to be chosen later. If we use the s variable, as in [9], the evolution equation for the curvature is written as
Note that we are able to use this equation even the curve is not necessarily concave. From this equation, we obtain the following differential equation for :
Let denote the closed curve consisting of and its mirror image below the x-axis. The length of equals and the area enclosed by equals . By the isoperimetric inequality (see, for example, [6, p. 33]), we have
If we choose , then
Now fix a time and define
Note that from (3.2). Suppose that at some time , where . We assume is the smallest time such that . Since implies that , we have
This contradicts the fact that W satisfies (3.3) and proves that
The following lower bound follows immediately from the isoperimetric inequality (3.4).
The inequalityholds for all.
Next we prove the upper bound of time derivatives.
Letbe the solution of the problem (P). For alland, the following estimates hold:
Lemmas 3.2–3.3 and (1.2) immediately give us an upper bound (3.5). Combining the identity (1.8) and (3.5), we obtain (3.6). By (1.7) and the Cauchy inequality, we obtain . This and (3.6) yield (3.7). □
Next we shall give a lower bound for and that is related to the height of the solution. Since the problem is non-local, we cannot use intersection number argument as in [12], we assume concavity for the initial data. The proof is similar to that of Lemma 2.2 in [12], except for this difference, but we write its proof for the reader’s convenience.
For any solutionof the problem (P), we havefor,, whereand.
We introduce
where is a positive function given by
Note that is non-decreasing, since is non-increasing. We define a differential operator by
By a simple calculation, we have
where the derivative is to be understood in a generalized sense. Let be the maximum point of the function , then
By substituting the boundary values and initial value, we can check that
This is a consequence of the following calculation. Let . Then, for ,
On the other hand, for , we have
The function is continuous, because of Lemma 3.2 and the regularity of the solution. Thus (3.12), (3.13) and (3.14) imply that on the parabolic boundary of the domain . Hence by (3.11) and the maximum principle,
Hence, by (3.13) and (3.15), we have for , which implies
where is given in (3.10). Similarly, replacing J by the function
and arguing as above in the region , we obtain
Now we introduce a function . Then a direct calculation shows
Note also that . Combining these with (3.16) and (3.17), and the comparison principle, we conclude
Therefore, the desired inequalities (3.8) and (3.9) follow from (1.2) and . The proof of Lemma 3.5 is complete. □
The next geometric lemma is the most technical step to prove Theorem 1.1. Here we again assume the concavity for the solution, which is not assumed in Proposition 6.2 of [12] for the shrinking case.
There exists a constantthat depends on initial curve such thatfor all.
We define the following positive constants:
Let . The inequality holds by the following argument. By the definition of and Lemma 3.3, we get
This implies
Thus, we have . Define
By (1.7) and (3.6), we have
We estimate the last term of (3.18). Let be the maximum point of at time t. By the concavity of the curve, the area of the triangle whose apexes are , and is smaller than the area that is enclosed by and x-axis. That implies
Furthermore, the concavity of the curve also gives us the inequality for . Combining those inequalities, (1.3) and
we obtain
By the Hölder inequality and the inequality
we have
If
then
This fact and (3.20) imply
Combining (3.19), (3.21) and (3.23), we have
Therefore,
This implies that is non-increasing if the right-hand side of this inequality is negative and (3.22) holds. Note that these two conditions are satisfied if is sufficiently large. Hence Lemma 3.6 has proved. □
By (3.9) and Lemmas 3.2 and 3.6, is uniformly bounded from below and negative. Thus and is uniformly bounded. By Theorem 2.1, the solution exists globally in time and its curvature and its length are uniformly bounded as . The uniform boundedness of higher derivatives follows from (2.1). □
Equation for the curvature
The time evolution of the curvature for strictly concave solution for is discussed in Section 2 of [12]. The argument here is quite similar to that of [12], except for the proof of the uniqueness of traveling wave for any given area parameter . From Lemma 3.2 in Section 3, the solution preserve concavity for if initial curve is concave, and the function , where , is well-defined. By repeating a similar calculation as [12], we have
where and the initial data satisfies the following additional condition:
The meaning of this condition is the following. In order that a solution of (4.1) to represents a curve satisfying , the above condition is needed. It is not difficult to show that
By considering a solution of (4.1) whose shape is independent of time, we obtain the following proposition for the existence of traveling wave solution.
For any givenand, there exist unique constants,and a unique concave functionthat satisfy (1.9). Furthermore, c satisfiesand the correspondence betweenandis one-to-one, whereis the area enclosed by the graph of the solutionand x-axis.
We shall consider only for concave curve in this paper. The eventual concavity property, which means that all solution of the problem (P) becomes concave for any large , will be investigated in the forthcoming paper. If we use the eventual concavity, we can prove the uniqueness of traveling wave without concavity assumption for initial data. We also remark that we cannot apply the eventual concavity results of [12], because their zero number arguments do not work for our non-local problem.
Under the assumption (4.2), and
Thus , where is some given positive constant. Therefore, a pair of the stationary solution and the constant c in (4.1)–(4.2) satisfies
Now we find and c that satisfy these two equations. Since , we obtain
The restriction (4.2) again gives us the following necessary condition:
This implies that if and only if so does the other inequality and the equality. Let us define
for any fixed . Without loss of generality, we can assume , thus . We can easily check that is monotone decreasing and is monotone increasing for . Also by the assumption . By taking limit and applying the intermediate theorem, we conclude that there exists a unique c for the pair . Consequently, for any given , there exists a unique constant c for which the solution of (4.3) satisfying the constraint (4.4) exists. Let us also mention that the constant is related with the solution by . Therefore, for any given length , there exists a unique stationary solution of (4.1) satisfying
Now, since and , the solution corresponding to the above stationary solution (4.3) satisfies
Consequently, , which means that is a traveling wave solution of the form , where is a solution of (1.9) and a is some constant. Conversely, if is a traveling wave solution of the form with a strictly concave function , then it is clear that the corresponding curvature function satisfies (4.3) along with the constraint (4.4). It is easy to check that the set of solutions for (4.6) is invariant under the transformation
for any . Note that and by this transformation, where is the area enclosed by the graph of the solution and the x-axis. Hence the correspondence between and is unique. In summary, for any given area , there exists a unique traveling wave solution up to translation. □
Spectral analysis
Our aim of this section is to prove the exponential stability of a traveling wave solution for the problem (P). For this purpose, we consider the following generalized problem of (1.1):
where H is defined by
Note that the problem (5.1) include the problem (1.1) when . If the solution curve is represented as a graph, the problem is reduced to the following free boundary problem:
For this equation we define , then it is easy to check that . Therefore, this equation also satisfies the area preserving property.
By using Eq. (5.8) in below and integrating by part, we can easily check that
Therefore the problem (5.1) has two Hamiltonians. One is the area function and the other is the height difference of the two endpoints.
The problem (5.1) can be also written as
Here we have used the equation
which is obtained by differentiating (5.4) with respect to t, and then using (5.3) and (5.5) along with the identity . The solution of the problem (5.7) exists locally in time. This can be checked by following the argument of Theorem 1.1. Furthermore, by the continuous dependence of semi-group to the parameter H, the solution of (5.7) exists globally in time provided that the parameter H is sufficiently close to zero [1,15].
We note that there is an important difference between (5.1) and (1.1). If we consider the problem (1.1) for , there is no traveling wave, which exists when . On the other hand, the problem (5.1) has a two parameter family of a traveling wave parametrized by H and L. Furthermore, the linearized operator associated to this problem has 0-eigenvalue that is semisimple with multiplicity two. Our strategy is embedding the problem (1.1) into (5.1) and prove the exponential convergence of the curvature. The problem (4.1) is then changed to
By a similar argument as Section 4, the curvature of a traveling wave is represented as
where c is a unique constant depends on H and L. For any initial data , the value is determined. From the structure of a family of traveling waves of the problem (1.1) for , the correspondence between and is one-to-one. In this section, for this given , we denote by to shorten the notation.
Let us explain that the problem (1.1) have traveling wave solution only when . The relations and (5.2) yield
Thus
and is needed for the existence of a traveling wave. On the other hand, for the problem (5.1), the above relations change to
By a similar calculation, the following condition is needed for the existence of a traveling wave
Let us define . Then this equation becomes
where . Here is a parameter that is invariant under the scaling transformation centered at the origin. The function is strictly monotone decreasing for ν and
By taking the limit, we can check that (5.10) has a solution for each as in Section 4. Furthermore, the parameter ν depends only on , and .
Now let us start to consider the linearized operator and its spectrum. We introduce the weighted space:
equipped with the inner product
The linearized operator of (5.8) at becomes the following:
In fact, by substituting to (5.8), and equating the coefficients of , we can check that
Here we used
to get the above operator. The domain of the operator is
where is the associated Sobolev space of . It is easy to check that the kernel of the operator is
by using the relations
By a simple calculation we can see that the following operator is the adjoint operator of :
where
The domain of this operator is
This operator seems to be easier to handle, since the boundary condition is not non-local. Moreover, the linearized operator is related to the linearization of distance function, and that is well analyzed in the field of differential geometry [7,16]. By using the similar argument as Lemma 2 of [8], we can show the operator is sectorial and the spectrum consists of infinite countably many eigenvalues with finite multiplicities. Since is adjoint operator of , the operator is also sectorial and it has infinite countable many eigenvalues. Thus, to consider the spectrum of the linear operator , we analyze the spectrums of the linear operator . Our approach is different from those of [7,16], since we do not need to introduce a technical curvelinear coordinate, which is invented to overcome the difficulty of the free boundary motion.
Letbe the support function of the traveling wave whose curvature is, and whose center is the left end point of the traveling wave. More precisely it is defined byThen, the kernel of the operatoris represented as the following:Furthermore, all non-zero eigenvalues are real and negative.
The boundary condition of for can be checked from the above representation formula (5.11). By the fundamental property of the support function [17], we obtain
By integrating this and using the boundary condition, we have
Therefore, holds. We can also see that the function belongs to , and by a simple calculation, we have
hence . To show , we introduce the following linear operator :
It is not difficult to check that is self-adjoint with respect to the inner product in and . Any must satisfy
Hence if . For the other case , Eq. (5.15) becomes
By the Fredholm alternative, (5.13), (5.14) and , the solution of is written as for any , which implies . Thus in both cases. Therefore, .
Now we prove the second statement. Let λ be a non-zero eigenvalue of , and be the corresponding eigenfunction. By a simple calculation, we can easily check that
Combining this with , we have
Therefore
The operator is self-adjoint with respect to the inner product in , hence we can apply the Sturm–Liouville theorem. Since the function is the simple 0-eigenfunction of and it has only one zero point in , the spectrum of consists of eigenvalues
Furthermore, if we set is the eigenfunction of the operator for the eigenvalue is , then is positive for every point on and .
Our aim is to prove that for any eigenfunction of the operator if . First we shall prove that any function can be represented by , where α, β are some constants and w is a linear combination of . All we need to prove is that . Then and span the same space spanned by . This can be checked as the following:
Therefore, there are constants , and a function , that is a linear combination of such that . Moreover, since the function is the 0-eigenfunction of , we can assume without loss of generality. By (5.16), we get
By applying this, we can prove that , since
Furthermore, by substituting into , we get
Now by (5.16), the following holds:
Finally we conclude that
Here we use the standard property of the support function explained in [17]. Therefore, and the lemma has proved. □
The spectrum ofconsists of countable many negative eigenvalues and zero. Moreover,and both geometric and algebraic multiplicity of the 0-eigenvalue are two.
We already showed is adjoint operator of with respect to the inner product in and has 0-eigenvalue and others are negative. Thus, also has 0-eigenvalue and negative ones. We have already checked that (5.17) holds and the geometric multiplicity of 0-eigenvalue is two. In order to prove the last statement, we assume that . Then there exists but . This condition implies . Thus, by (5.17), there exists such that . By (5.12), we have
The function in that is perpendicular to is proportional to . This can be proved by the following calculation:
Hence for some , but this is impossible since
This contradiction shows that the algebraic multiplicity is two and the lemma has proved. □
Letbe the curvature of a curve whose endpoints are on x-axis and the area enclosed by the curve and x-axis is. Letbe a curvature of a traveling wave solution whose endpoints are on x-axis having the same area. Assume thatis sufficiently close touniformly on, then the corresponding solutionof the problem (4.1) converges toexponentially onuniformly as.
We will apply the general theory Theorem 2.1 of [13] and Proposition 5.2. We set and of [13] by and . The set of the stationary solutions for the problem (5.1) is
where is a function that depends on H and L. Define a function by
The function F is smooth with respect to μ, H and L, and from the proof in Remark 5.2, . By the implicit function theorem, the function that satisfies
is a smooth function with respect to H and L. It is obvious that is smooth in the domain and dimension of the above manifold is two in the space . As is easily seen, the constant L corresponds to the scaling of a traveling wave at some fixed point on x-axis and H corresponds to the parameter which changes the difference of the height of the two endpoints. To calculate the tangent space of the set of stationary solutions of (5.8) at in the space , we differentiate with respect to H and L at the point . Then we have
Now, we check that and are linearly independent. If this does not hold, by comparing the first term of right-hand side of (5.19) and (5.20), we have
By differentiating (5.18) with respect to H and L at , we also have
It is easy to see , and hence combining this and (5.21), we have
By a simple calculation, we also have
Substituting these equations into (5.22) and using , we see that the following equality must hold:
By substituting (5.9) into this, we get
This is a contradiction. Thus, as we mentioned above, the tangent space of the set of stationary solutions of (5.8) is spanned by and , which is the same space spanned by 1 and . Combining this with Proposition 5.2, we have already shown that
the set of all equilibria of (5.8) is a smooth 2-dimensional manifold in ,
the tangent space of at is , and the eigenvalue 0 is semisimple,
.
Combining these and Remark 2.1, we can apply Theorem 2.1 of [13] to conclude that converges to some element in . We have already known that there exist and such that converges to exponentially uniformly on as . Since (4.2) holds and is a preserved function, the following equations must be satisfied:
and the correspondence between and is unique, hence and . Therefore, converges to exponentially on uniformly as , and the theorem has proven. □
The reader might think that, we need to consider the spectrum of the linearized operator around for Eq. (4.1). This is one strategy to prove Theorem 1.2, and it is possible to prove that all eigenvalues are negative or zero. Unfortunately, the zero eigenvalue is not semisimple and we cannot apply the stability theory to get the exponential convergence to . In order to overcome this difficulty, we consider the problem (5.1).
Now we go back to the proof of our main theorem.
By differentiating (1.3) with respect to t, and then using (1.2) and (1.4), we get
Theorem 5.1 implies that κ converges to exponentially fast. Thus and for some as . In particular,
By integrating , we obtain the convergence of as . Hence, we can define . Let us reconstruct the solution curve from the curvature by the relation:
By a similar way, we again construct the curve for the graph of from the curvature , whose shape represents a traveling wave solution:
Therefore, as , we conclude
uniformly for . Thus, boundedness of and the mean value theorem imply
as as long as . This and (5.23)–(5.24) give us pointwise convergence on as , where we consider zero extension as in the statement. Combining these with space derivative bounds, we conclude the uniform convergence. □
Footnotes
Acknowledgement
This work is partially supported by JSPS KAKENHI Grant no. 23740128.
References
1.
H.Amann, Dynamic theory of quasilinear parabolic equations II. Reaction diffusion systems, Differential and Integral Equations3 (1990), 13–75.
2.
S.Angenent, On the formation of singularities in the curve shortening flow, J. Diff. Geom.33 (1991), 601–633.
3.
Y.-L.Chang, J.-S.Guo and Y.Kohsaka, On a two-point free boundary problem for a quasilinear parabolic equation, Asymptotic Analysis34 (2003), 333–358.
4.
X.-L.Chao, X.-R.Ling and X.-L.Wang, On a planar area-preserving curvature flow, Proceedings of the American Mathematical Society141 (2013), 1783–1789.
5.
X.Chen and J.-S.Guo, Motion by curvature of planar curves with end points moving freely on a line, Math. Ann.350 (2011), 277–311.
6.
M.Do Carmo, Differential Geometry of Curves and Surfaces, Prentice-Hall, Englewood Cliffs, NJ, 1976.
7.
S.-I.Ei, M.Sato and E.Yanagida, Stability of stationary interfaces with contact angle in a generalized mean curvature flow, American J. Math.118 (1996), 653–687.
8.
J.Escher and Z.Feng, Exponential stability of equilibria of the curve shortening flow with contact angle, Dynamics of Continuous, Discrete and Impulsive Systems. Series A, Mathematical Analysis14 (2007), 287–299.
9.
M.Gage, On an area-preserving evolution equation for plane curves, in: Nonlinear Problems in Geometry, Contemporary Mathematics, Vol. 51, 1986, pp. 51–62.
10.
M.Gage and R.S.Hamilton, The heat equation shrinking convex plane curves, J. Diff. Geom.23 (1986), 69–96.
11.
M.A.Grayson, The heat equation shrinks embedded plane curves to round points, J. Diff. Geom.26 (1987), 285–314.
12.
J.-S.Guo, H.Matano, M.Shimojo and C.-H.Wu, On a free boundary problem for the curvature flow with driving force, Arch. Ration. Mech. Anal. (2015), to appear.
13.
J.Prüss, G.Simonett and R.Zacher, On convergence of solutions to equilibria for quasilinear parabolic problems, Journal of Differential Equations246 (2009), 3902–3931.
14.
G.Sapiro and A.Tannenbaum, Area and length preserving geometric invariant scale-spaces, IEEE Transactions on Pattern Analysis and Machine Intelligence17 (1995), 67–72.
15.
G.Simonett, Invariant manifolds and bifurcation for quasilinear reaction–diffusion systems, Nonlinear Analysis, Theory, Methods and Applications23 (1994), 515–544.
16.
T.I.Vogel, Sufficient conditions for capillary surfaces to be energy minima, Pacific Journal of Mathematics194 (2000), 469–489.
17.
X.P.Zhu, Lectures on Mean Curvature Flows, American Mathematical Society, Studies in Advanced Mathematics, 2002.