Home Global existence for the Willmore flow with boundary via Simon’s Li–Yau inequality
Article
Licensed
Unlicensed Requires Authentication

Global existence for the Willmore flow with boundary via Simon’s Li–Yau inequality

  • Manuel Schlierf ORCID logo EMAIL logo
Published/Copyright: January 29, 2025

Abstract

It is well known that the Willmore flow of closed spherical immersions exists globally in time and converges if the initial datum has Willmore energy below 8 π , exactly the Li–Yau energy threshold below which all closed immersions are embedded. Extending the Li–Yau inequality for closed surfaces via Simon’s monotonicity formula also for surfaces with boundary, given Dirichlet boundary conditions, one obtains an energy threshold C LY below which surfaces with this boundary are embedded. By a slight modification, one obtains a threshold C LY rot below which surfaces of revolution satisfying the boundary data have no self-intersections on the rotation axis. With a new argument, using this modified Li–Yau inequality and tools from geometric measure theory, we show that the Willmore flow with Dirichlet boundary data starting in cylindrical surfaces of revolution exists globally in time if the energy of the initial datum is below C LY rot . Moreover, given Dirichlet boundary data, we also obtain the existence of a Willmore minimizer in the class of cylindrical surfaces of revolution if the corresponding infimum lies below C LY rot , which improves previous results for the stationary problem.

MSC 2020: 53E40; 49Q10; 35B40; 35K41; 35J35

A Immersions into H 2 with Sobolev regularity

The above definitions of the curvature, tangent vector and elastic energy also make sense for immersions u W 2 , 2 ( I , H 2 ) , where I R is a compact interval. In this article, W 2 , 2 ( I , H 2 ) : = { v W 2 , 2 ( I , R 2 ) : v ( I ) H 2 } is to be understood as an open subset of the Sobolev space W 2 , 2 ( I , R 2 ) C 1 ( I , R 2 ) .

A curve u W 2 , 2 ( I , H 2 ) is called immersion if | x u | > 0 . For such 𝑢, s = 1 / | x u | x denotes the (weak) derivative with respect to the arc-length parameter 𝑠. By the chain rule for Sobolev functions,

(A.1) s 2 u = ( u ( 2 ) | x u | ) 2 x 2 u ( u ( 2 ) | x u | ) 2 x 2 u , x u | x u | x u | x u | + u ( 2 ) x u ( 2 ) | x u | x u | x u | L 2 ( I , R 2 )

exists and, as in [8, (12)], one computes that

(A.2) κ = s u s u = s 2 u + 1 u ( 2 ) ( 2 s u ( 1 ) s u ( 2 ) ( s u ( 1 ) ) 2 ( s u ( 2 ) ) 2 ) .

Moreover, as already computed in [12, (2.9)],

(A.3) E ( u ) = 0 1 ( x 2 u ( 2 ) x u ( 1 ) u ( 2 ) x 2 u ( 1 ) x u ( 2 ) u ( 2 ) + x u ( 1 ) | x u | 2 ) 2 1 u ( 2 ) | x u | 5 d x .

Lemma A.1

Consider immersions u n , u W 2 , 2 ( I , H 2 ) for n N with u n u in W 2 , 2 ( I ) . Then E ( u n ) E ( u ) .

Proof

The claim is due to (A.3) and the compact embedding W 2 , 2 ( I ) C 1 ( I , R 2 ) . ∎

Lemma A.2

For immersions u n , u W 2 , 2 ( I , H 2 ) for n N with u n u in W 2 , 2 ( I ) , one has

E ( u ) lim inf n E ( u n ) .

Proof

After passing to a subsequence, without loss of generality, u n u in C 1 ( I ) and

lim n E ( u n ) = lim inf n E ( u n ) .

It is a well-known fact that the product of a uniformly convergent sequence and a sequence that weakly converges in L 2 is again weakly convergent in L 2 to the product of the respective limits. Thus, the claim is a consequence of (A.3) and the weak lower semi-continuity of L 2 ( I ) . ∎

B Energy estimates

We use the following two technical results from [31].

Proposition B.1

Consider immersions u n : I H 2 with u n W 2 , 2 ( I ) defined on a compact interval I R such that u n ( 2 ) L ( I ) 0 and L R 2 ( u n ) > 0 . Then lim sup n E ( u n ) = .

Proof

Proposition 2.7 in [31] yields the claim for smooth immersions u n . Using [22] to suitably approximate each u n W 2 , 2 ( I ) and using Lemma A.1, one obtains the claim. ∎

Theorem B.2

Theorem B.2 ([31, Theorem 2.5])

Consider smooth immersions u n : I H 2 on a compact interval I R such that 0 < | x u n | L < . If ( | x u n | ) n N converges uniformly on 𝐼, u n u in W 1 , ( I ) and u n u uniformly on 𝐼 and

sup n N E ( u n ) E 8 + η for some E N 0 and η < 8 ,

then Z = { u ( 2 ) = 0 } I consists of at most 𝐸 points.

C Monotonicity formula and estimates for C LY ( p , τ )

As it turns out, the underlying monotonicity formula in [25, Appendix B] from which Proposition 2.3 is concluded gives very useful exact identities for surfaces contained in spheres. Since the integral in (2.5) only depends on the Dirichlet boundary conditions and as the Dirichlet boundary conditions of cylindrical surfaces of revolution are parallel circles with prescribed co-normal field, C LY ( p , τ ) can be estimated using suitable cut-outs of spheres.

First, recall Simon’s monotonicity formula in [25, Appendix B] for varifolds with boundary, cf. [32, (1.2)] and [18, Appendix A] for the versions without boundary. Consider a 2-dimensional integral varifold 𝑉 in R 3 with generalized mean curvature H and W ( V ) , μ V ( R 3 ) < . For any z R 3 such that

B 1 ( z ) | x z , η V ( x ) | | x z | 2 d σ V ( x ) < ,

writing for the orthogonal projection onto the orthogonal complement of the approximate tangent space, one has

(C.1) A z ( σ ) + B ρ ( z ) B σ ( z ) | H ( x ) 2 + ( x z ) | x z | 2 | 2 d μ V ( x ) = A z ( ρ )

for all 0 < σ < ρ , where, for t > 0 ,

A z ( t ) = μ V ( B t ( z ) ) t 2 + 1 4 B t ( z ) | H | 2 d μ V + B t ( z ) H ( x ) , x z t 2 d μ V ( x ) + 1 2 B t ( z ) ( 1 | x z | 2 1 t 2 ) x z , η V ( x ) d σ V ( x ) .

For y = 0 , 1 , consider points p y H 2 and unit tangent vectors τ y S 1 which are not vertical, i.e. τ y ( 1 ) 0 . The result below remains true in the cases where τ y = ± ( 0 , 1 ) is allowed and can be proved similarly by using cut-outs of planes instead of spheres. Recall the energy threshold T ( τ ) = 4 π 2 π τ y ( 2 ) | y = 0 1 in [30].

Proposition C.1

One has that C LY ( p , τ ) T ( τ ) .

To this end, we make the following construction. First, consider the parallel circles

C y = { ( p y ( 1 ) , p y ( 2 ) cos ( θ ) , p y ( 2 ) sin ( θ ) ) : θ R } for y = 0 , 1 .

Consider now an integer rectifiable 2-varifold V y with generalized mean curvature H y , finite mass and Willmore energy such that, for

η y ( p y ( 1 ) , p y ( 2 ) cos θ , p y ( 2 ) sin θ ) = ( 1 ) 1 + y ( τ y ( 1 ) , τ y ( 2 ) cos θ , τ y ( 2 ) sin θ ) for θ R ,

one has σ V y = H 1 C y and η V y = η y . Using [25, (B.8)] for the case z C y and [25, (4.2)] for z C y , one finds

lim t 0 A z ( t ) = θ 2 ( μ V y , z ) π for all z R 3 .

Moreover, estimating B t ( z ) H y , x z t 2 d μ V y as in [25, p. 555] and using μ V y ( R 3 ) < , one finds

lim t A z ( t ) = W ( V y ) 4 + 1 2 R 3 η V y ( x ) , x z | x z | 2 d σ V y ( x ) .

Using (C.1), one concludes

(C.2) θ 2 ( μ V y , z ) π + R 3 | H y 2 + ( x z ) | x z | 2 | 2 d μ V y = W ( V y ) 4 + 1 2 R 3 η V y , x z | x z | 2 d σ V y .

We now make a specific choice for V y . There exist R y > 0 , x y R × { 0 } × { 0 } and a smooth 2-dimensional submanifold M y S R y , x y 2 : = x y + R y S 2 R 3 such that V y : = v ( M y ) satisfies σ V y = H 1 C y , η V y = η y . Using (1.6) and (C.2), one finds

(C.3) C LY ( p , τ ) = 8 π + inf z R 3 y = 0 , 1 ( W ( V y ) 4 θ 2 ( μ V y , z ) π 4 R 3 | H y 2 + ( x z ) | x z | 2 | 2 d μ V y ) = 8 π + W ( V 0 ) + W ( V 1 ) 4 sup z R 3 y = 0 , 1 ( θ 2 ( μ V y , z ) π + R 3 | H y 2 + ( x z ) | x z | 2 | 2 d μ V y ) .

Lemma C.2

One has W ( V 0 ) + W ( V 1 ) = T ( τ ) .

Proof

First, consider any τ S 1 with τ ( 1 ) 0 and the spherical cut-out M S 2 as above with outer co-normal induced by rotating 𝜏 as in the construction above.

Case 1: τ ( 1 ) < 0 . Choosing θ ( 0 , π ) such that

( sin ( θ ) cos ( θ ) ) = τ ,

we can parametrize 𝑀 by the surface of revolution associated to γ ( x ) = ( cos ( x ) , sin ( x ) ) with x [ 0 , θ ) . Since the mean curvature of S 2 has length 1 everywhere,

(C.4) W ( M ) = H 2 ( M ) = 2 π 0 θ γ ( 2 ) ( x ) | γ ( x ) | d x = 2 π ( 1 cos ( θ ) ) = 2 π ( 1 τ ( 2 ) ) .

Case 2: τ ( 1 ) > 0 . With similar arguments, one also obtains (C.4).

Using that the Willmore energy is invariant with respect to scaling and translations, that τ 0 induces the outer co-normal of V 0 and that τ 1 induces the outer co-normal of V 1 , (C.4) yields

W ( V 0 ) + W ( V 1 ) = 2 π ( 1 + 1 τ 1 ( 2 ) + τ 0 ( 2 ) ) = 4 π 2 π τ y ( 2 ) | y = 0 1 = T ( τ ) .

Lemma C.3

It holds that

sup z R 3 y = 0 , 1 ( θ 2 ( μ V y , z ) π + R 3 | H y ( x ) 2 + ( x z ) | x z | 2 | 2 d μ V y ( x ) ) 2 π .

Proof

To this end, simply note that (C.2) applied to S R y , x y 2 , a manifold without boundary, yields, for y = 0 , 1 ,

θ 2 ( μ V y , z ) π + R 3 | H y ( x ) 2 + ( x z ) | x z | 2 | 2 d μ V y ( x ) θ 2 ( H 2 S R y , x y 2 , z ) π + S R y , x y 2 | H S R y , x y 2 ( x ) 2 + ( x z ) | x z | 2 | 2 d H 2 ( x ) = W ( S R y , x y 2 ) 4 = π ,

using that the support of V y is contained in the sphere S R y , x y 2 . ∎

Proof of Proposition C.1

The claim is concluded by combining (C.3) with both lemmata above. ∎

D Computation of C LY rot ( p , τ ) and comparison to [11]

Returning to (1.6), for z = h e 1 , where h R and e 1 = ( 1 , 0 , 0 ) , by direct computation,

(D.1) C η τ ( x ) , b p ( x ) z | b p ( x ) z | 2 d μ b p , ( x ) = 2 π [ p y ( 2 ) τ y ( 1 ) ( p y ( 1 ) h ) + τ y ( 2 ) p y ( 2 ) ( p y ( 1 ) h ) 2 + ( p y ( 2 ) ) 2 | y = 0 1 ] ,

and plugging this into (1.6) yields formula (1.7). In particular, if one considers the case of horizontal boundary data, i.e. τ 0 = τ 1 = ( 1 , 0 ) , p 0 = ( 1 , α ) and p 1 = ( 1 , α + ) for α , α + > 0 , one finds

C LY rot ( p , τ ) = 4 π [ 2 sup h R ( α + ( 1 h ) ( 1 h ) 2 + ( α + ) 2 + α ( 1 + h ) ( 1 + h ) 2 + ( α ) 2 ) ] = 8 π + 4 π inf h R ( α + ( h 1 ) ( 1 h ) 2 + ( α + ) 2 α ( 1 + h ) ( 1 + h ) 2 + ( α ) 2 ) ,

also cf. [11, Theorem 1.2]. Next, adapting our methodology from the previous section, using Simon’s monotonicity formula, we show that the energy thresholds in (1.9) and [11, Theorem 1.1] are equivalent for all possible boundary conditions.

To this end, we rewrite the main players in defining the threshold in [11, Theorem 1.1] with this article’s notation. If p 0 , p 1 H 2 and τ 0 , τ 1 S 1 , for h R , denote by c y h : [ 0 , 1 ) H 2 for y = 0 , 1 a parametrization of a segment of the profile curve of a sphere, a disk or of a Möbius transformed inverted catenoid such that, for y = 0 , 1 ,

(D.2) c y h ( 0 ) = p y , lim x 1 c y h ( x ) = ( h , 0 ) and x c y h ( 0 ) | x c y h ( 0 ) | = ( 1 ) y τ y .

We then have the following.

Proposition D.1

Writing Γ y = { y } × S 1 , one finds that

(D.3) W ( f c y h ) = 4 π 2 Γ y η τ ( x ) , b p ( x ) h e 1 | b p ( x ) h e 1 | 2 d μ b p , ( x )

for y = 0 , 1 , using (1.4), where f c y h is the immersion one obtains by rotating c y h .

Especially, summing over y = 0 , 1 in (D.3), one finds

W ( f c 0 h ) + W ( f c 1 h ) = 8 π 2 C η τ ( x ) , b p ( x ) h e 1 | b p ( x ) h e 1 | 2 d μ b p , ( x ) .

That is, recalling (1.6), the energy thresholds of Theorem 1.4 and [11, Theorem 1.1] coincide.

Proof of Proposition D.1

Suppose first of all that c y h parametrizes a segment of the profile curve of a Möbius transformed inverted catenoid. Moreover, denote by V c y h the varifold one obtains by rotating c y h and by V c y h cl the associated inverted catenoid without boundary whose support contains the support of V c y h . Using (D.2), one finds θ 2 ( μ V c y h cl , h e 1 ) π = 2 π , i.e. the closed inverted catenoid has density 2 at the point h e 1 where it intersects the rotation axis. Moreover, as W ( V c y h cl ) = 8 π , proceeding as in (C.2) yields

2 π + R 3 | H V c y h cl ( x ) 2 + ( x h e 1 ) | x h e 1 | 2 | 2 d μ V c y h cl ( x ) = 2 π .

In particular,

0 R 3 | H V c y h ( x ) 2 + ( x h e 1 ) | x h e 1 | 2 | 2 d μ V c y h ( x ) R 3 | H V c y h cl ( x ) 2 + ( x h e 1 ) | x h e 1 | 2 | 2 d μ V c y h cl ( x ) = 0 .

By (D.2), we have θ 2 ( μ V c y h , h e 1 ) π = π , i.e. the varifold obtained by rotating c y h only has density equal to 1 at the point where it intersects the rotation axis. Therefore, (C.2) applied to V c y h yields (D.3), using (1.4). With the very same argument, one proves the claim also in the case where c y h parametrizes a segment of the profile curve of a sphere.

Finally, if c y h parametrizes the profile curve of a disk, then W ( f c y h ) = 0 . Moreover, using (D.2), τ y = ( 1 ) y ( 0 , 1 ) and p y ( 1 ) = h . Thus, as in (D.1), using (1.4),

Γ y η τ ( x ) , b p ( x ) h e 1 | b p ( x ) h e 1 | 2 d μ b p , ( x ) = 2 π ( 1 ) y + 1 ( p y ( 2 ) τ y ( 1 ) ( p y ( 1 ) h ) + τ y ( 2 ) p y ( 2 ) ( p y ( 1 ) h ) 2 + ( p y ( 2 ) ) 2 ) = 2 π ( 1 ) y + 1 τ y ( 2 ) = 2 π

and (D.3) also follows in this case. ∎

Acknowledgements

The author would like to thank Anna Dall’Acqua for helpful discussions and Sascha Eichmann for his comments on comparing Theorem 1.4 and [11, Theorem 1.1]. Moreover, the author is grateful to the referees for their valuable comments on the original manuscript.

  1. Communicated by: Francesca Da Lio

References

[1] W. K. Allard, On the first variation of a varifold, Ann. of Math. (2) 95 (1972), 417–491. 10.2307/1970868Search in Google Scholar

[2] R. Bryant, A duality theorem for Willmore surfaces, J. Differential Geom. 20 (1984), no. 1, 23–53. 10.4310/jdg/1214438991Search in Google Scholar

[3] R. Bryant and P. Griffiths, Reduction for constrained variational problems and 1 2 k 2 d s , Amer. J. Math. 108 (1986), no. 3, 525–570. 10.2307/2374654Search in Google Scholar

[4] P. Canham, The minimum energy of bending as a possible explanation of the biconcave shape of the human red blood cell, J. Theoret. Biol. 26 (1970), no. 1, 61–81. 10.1016/S0022-5193(70)80032-7Search in Google Scholar PubMed

[5] A. Dall’Acqua, K. Deckelnick and H.-C. Grunau, Classical solutions to the Dirichlet problem for Willmore surfaces of revolution, Adv. Calc. Var. 1 (2008), no. 4, 379–397. 10.1515/ACV.2008.016Search in Google Scholar

[6] A. Dall’Acqua, S. Fröhlich, H.-C. Grunau and F. Schieweck, Symmetric Willmore surfaces of revolution satisfying arbitrary Dirichlet boundary data, Adv. Calc. Var. 4 (2011), no. 1, 1–81. 10.1515/acv.2010.022Search in Google Scholar

[7] A. Dall’Acqua, M. Müller, R. Schätzle and A. Spener, The Willmore flow of tori of revolution, Anal. PDE 17 (2024), no. 9, 3079–3124. 10.2140/apde.2024.17.3079Search in Google Scholar

[8] A. Dall’Acqua and A. Spener, The elastic flow of curves in the hyperbolic plane, preprint (2017), https://arxiv.org/abs/1710.09600. Search in Google Scholar

[9] A. Dall’Acqua and A. Spener, Circular solutions to the elastic flow in hyperbolic space, RIMS Kôkyûroku 2082 (2018), 109–124. Search in Google Scholar

[10] S. Eichmann, Nonuniqueness for Willmore surfaces of revolution satisfying Dirichlet boundary data, J. Geom. Anal. 26 (2016), no. 4, 2563–2590. 10.1007/s12220-015-9639-xSearch in Google Scholar

[11] S. Eichmann, Improved long time existence for the Willmore flow of surfaces of revolution with Dirichlet data, preprint (2024), https://arxiv.org/abs/2402.05580. Search in Google Scholar

[12] S. Eichmann and H.-C. Grunau, Existence for Willmore surfaces of revolution satisfying non-symmetric Dirichlet boundary conditions, Adv. Calc. Var. 12 (2019), no. 4, 333–361. 10.1515/acv-2016-0038Search in Google Scholar

[13] S. Eichmann and R. Schätzle, The rotational symmetric Willmore boundary problem, preprint (2024), https://www.math.uni-tuebingen.de/user/schaetz/publi/eichmann-schaetzle-24.pdf. Search in Google Scholar

[14] S. W. Hawking, Gravitational radiation in an expanding universe, J. Math. Phys. 9 (1968), no. 4, 598–604. 10.1063/1.1664615Search in Google Scholar

[15] W. Helfrich, Elastic properties of lipid bilayers: Theory and possible experiments, Z. Naturforschung C 28 (1973), no. 11–12, 693–703. 10.1515/znc-1973-11-1209Search in Google Scholar PubMed

[16] E. Kuwert and R. Schätzle, The Willmore flow with small initial energy, J. Differential Geom. 57 (2001), no. 3, 409–441. 10.4310/jdg/1090348128Search in Google Scholar

[17] E. Kuwert and R. Schätzle, Gradient flow for the Willmore functional, Comm. Anal. Geom. 10 (2002), no. 2, 307–339. 10.4310/CAG.2002.v10.n2.a4Search in Google Scholar

[18] E. Kuwert and R. Schätzle, Removability of point singularities of Willmore surfaces, Ann. of Math. (2) 160 (2004), no. 1, 315–357. 10.4007/annals.2004.160.315Search in Google Scholar

[19] J. Langer and D. A. Singer, The total squared curvature of closed curves, J. Differential Geom. 20 (1984), no. 1, 1–22. 10.4310/jdg/1214438990Search in Google Scholar

[20] J. M. Lee, Introduction to Riemannian Manifolds, 2nd ed, Grad. Texts in Math. 176, Springer, Cham, 2018. 10.1007/978-3-319-91755-9Search in Google Scholar

[21] P. Li and S. T. Yau, A new conformal invariant and its applications to the Willmore conjecture and the first eigenvalue of compact surfaces, Invent. Math. 69 (1982), no. 2, 269–291. 10.1007/BF01399507Search in Google Scholar

[22] N. G. Meyers and J. Serrin, H = W , Proc. Natl. Acad. Sci. USA 51 (1964), 1055–1056. 10.1073/pnas.51.6.1055Search in Google Scholar PubMed PubMed Central

[23] G. Mingione, Topics in Modern Regularity Theory, Centro di Ric. Mat. Ennio De Giorgi (CRM) Ser. 13, Edizioni della Normale, Pisa, 2012. 10.1007/978-88-7642-427-4Search in Google Scholar

[24] M. Müller and A. Spener, On the convergence of the elastic flow in the hyperbolic plane, Geom. Flows 5 (2020), no. 1, 40–77. 10.1515/geofl-2020-0002Search in Google Scholar

[25] M. Novaga and M. Pozzetta, Connected surfaces with boundary minimizing the Willmore energy, Math. Eng. 2 (2020), no. 3, 527–556. 10.3934/mine.2020024Search in Google Scholar

[26] M. Pozzetta, On the Plateau–Douglas problem for the Willmore energy of surfaces with planar boundary curves, ESAIM Control Optim. Calc. Var. 27 (2021), Paper No. S2. 10.1051/cocv/2020049Search in Google Scholar

[27] F. Rupp and C. Scharrer, Li-Yau inequalities for the Helfrich functional and applications, Calc. Var. Partial Differential Equations 62 (2023), no. 2, Paper No. 45. 10.1007/s00526-022-02381-7Search in Google Scholar PubMed PubMed Central

[28] C. Scharrer, Some geometric inequalities for varifolds on Riemannian manifolds based on monotonicity identities, Ann. Global Anal. Geom. 61 (2022), no. 4, 691–719. 10.1007/s10455-021-09822-0Search in Google Scholar

[29] R. Schätzle, The Willmore boundary problem, Calc. Var. Partial Differential Equations 37 (2010), no. 3–4, 275–302. 10.1007/s00526-009-0244-3Search in Google Scholar

[30] M. Schlierf, On the convergence of the Willmore flow with Dirichlet boundary conditions, Nonlinear Anal. 241 (2024), Article ID 113475. 10.1016/j.na.2023.113475Search in Google Scholar

[31] M. Schlierf, Singularities of the hyperbolic elastic flow: convergence, quantization and blow-ups, Calc. Var. Partial Differential Equations 63 (2024), no. 8, Paper No. 207. 10.1007/s00526-024-02815-4Search in Google Scholar

[32] L. Simon, Existence of surfaces minimizing the Willmore functional, Comm. Anal. Geom. 1 (1993), no. 2, 281–326. 10.4310/CAG.1993.v1.n2.a4Search in Google Scholar

[33] A. Volkmann, A monotonicity formula for free boundary surfaces with respect to the unit ball, Comm. Anal. Geom. 24 (2016), no. 1, 195–221. 10.4310/CAG.2016.v24.n1.a7Search in Google Scholar

Received: 2024-02-19
Accepted: 2024-11-08
Published Online: 2025-01-29
Published in Print: 2025-04-01

© 2025 Walter de Gruyter GmbH, Berlin/Boston

Downloaded on 8.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/acv-2024-0018/html
Scroll to top button