We extend to the anisotropic setting the existence of solutions for the Kirchhoff-Plateau problem and its dimensional reduction.
Citation: Antonio De Rosa, Luca Lussardi. On the anisotropic Kirchhoff-Plateau problem[J]. Mathematics in Engineering, 2022, 4(2): 1-13. doi: 10.3934/mine.2022011
[1] | Miyuki Koiso . Stable anisotropic capillary hypersurfaces in a wedge. Mathematics in Engineering, 2023, 5(2): 1-22. doi: 10.3934/mine.2023029 |
[2] | Juan Pablo Borthagaray, Wenbo Li, Ricardo H. Nochetto . Finite element algorithms for nonlocal minimal graphs. Mathematics in Engineering, 2022, 4(2): 1-29. doi: 10.3934/mine.2022016 |
[3] | Sisto Baldo, Van Phu Cuong Le, Annalisa Massaccesi, Giandomenico Orlandi . Energy minimizing maps with prescribed singularities and Gilbert-Steiner optimal networks. Mathematics in Engineering, 2023, 5(4): 1-19. doi: 10.3934/mine.2023078 |
[4] | Weimin Sheng, Shucan Xia . Interior curvature bounds for a type of mixed Hessian quotient equations. Mathematics in Engineering, 2023, 5(2): 1-27. doi: 10.3934/mine.2023040 |
[5] | Joan Mateu, Maria Giovanna Mora, Luca Rondi, Lucia Scardia, Joan Verdera . A maximum-principle approach to the minimisation of a nonlocal dislocation energy. Mathematics in Engineering, 2020, 2(2): 253-263. doi: 10.3934/mine.2020012 |
[6] | Chao Xia, Jiabin Yin . Two overdetermined problems for anisotropic $ p $-Laplacian. Mathematics in Engineering, 2022, 4(2): 1-18. doi: 10.3934/mine.2022015 |
[7] | Barbara Brandolini, Florica C. Cîrstea . Anisotropic elliptic equations with gradient-dependent lower order terms and $ L^1 $ data. Mathematics in Engineering, 2023, 5(4): 1-33. doi: 10.3934/mine.2023073 |
[8] | Stefano Berrone, Stefano Scialò, Gioana Teora . The mixed virtual element discretization for highly-anisotropic problems: the role of the boundary degrees of freedom. Mathematics in Engineering, 2023, 5(6): 1-32. doi: 10.3934/mine.2023099 |
[9] | Catharine W. K. Lo, José Francisco Rodrigues . On an anisotropic fractional Stefan-type problem with Dirichlet boundary conditions. Mathematics in Engineering, 2023, 5(3): 1-38. doi: 10.3934/mine.2023047 |
[10] | Chiara Gavioli, Pavel Krejčí . Deformable porous media with degenerate hysteresis in gravity field. Mathematics in Engineering, 2025, 7(1): 35-60. doi: 10.3934/mine.2025003 |
We extend to the anisotropic setting the existence of solutions for the Kirchhoff-Plateau problem and its dimensional reduction.
The study of minimal surfaces spanning elastic boundaries dates back to Courant [13] and Lewy [33]. They studied the Plateau problem under the assumption that the boundary of the minimal surface is not fixed, but is constrained to lie on a prescribed manifold. The generalization to minimal surfaces spanning non-constrained elastic boundaries has been recently addressed by Giomi & Mahadevan [26]. These results have been complemented with an investigation of the stability of flat circular solutions by Chen & Fried [11], Biria & Fried [9,10], Giusteri, Franceschini & Fried [27], and Hoang & Fried [32]. A similar problem has been treated by Bernatzky & Ye [6] employing the theory of currents, however the elastic energy used therein fails to satisfy the physical requirement of invariance under superposed rigid transformations.
The Kirchhoff-Plateau problem differs from the aforementioned works, because the spanning boundary is assumed to lay on the surface of an elastic loop, referred to as the rod, which is modeled as a deformable manifold. On the contrary, in all of the studies above the boundary of the spanning surface was assumed to coincide with the loop midline. In the Kirchhoff-Plateau problem the filament forming the loop is assumed to be thin enough to be modeled faithfully by a Kirchhoff rod, that is an unshearable inextensible rod which can sustain bending of its midline and twisting of its cross-sections, see Antman's work [5]. To model the flexible rods, some physical constraints are imposed, such as local and global non-interpenetration of matter introduced by Schuricht [34]. The isotropic Kirchhoff-Plateau problem, that is minimizing the area functional, has been investigated by Giusteri, Lussardi & Fried [28] with only one filament, and by Bevilacqua, Lussardi & Marzocchi [7] taking into account a system of linked rods. The authors utilize the boundary condition via linking number introduced by Harrison [30] and further investigated by De Lellis, Ghiraldin & Maggi [15]. Moreover, a dimensional reduction of the aforementioned variational problem has been treated by Bevilacqua, Lussardi & Marzocchi [8].
In view of the works of Almgren [3,4], Taylor [36] and Allard [1,2], a natural question is whether the isotropic results [7,8] generalize to anisotropic surface energies. Indeed, an increasing interest has been recently devoted to the study of the anisotropic Plateau problem: see for instance the results by De Philippis, De Rosa & Ghiraldin [16,17,18], De Rosa [20], De Rosa & Kolasinski [23] and Harrison & Pugh [31]. We also refer the reader to [14,19,21,22,24].
The aim of this paper is to address this question, considering the anisotropic Kirchhoff-Plateau problem for systems of linked rods. The energy functional we minimize is given by the sum of the elastic and the potential energy for the link and the anisotropic surface energy of the film. As for the isotropic Kirchhoff-Plateau problem we prescribe the linking type of the system of rods as well as the non-interpenetration of matter for each rod. Furthermore, each midline has a prescribed knot-type.
To conclude, we perform a dimensional reduction in the spirit of the analysis carried out in the isotropic setting.
In this section we recall notation for the geometry of curves. If x1,x2:[0,L]→R3 are two continuous and closed curves, their linking number is the integer value
Link(x1,x2):=14π∫L0∫L0x1(s)−x2(t)|x1(s)−x2(t)|3⋅x1′(s)×x2′(t)dsdt. |
We say that x1 and x2 are isotopic, and we use the notation x1≃x2, if there exists an open neighborhood N1 of x1([0,L]), an open neighborhood N2 of x2([0,L]) and a continuous map Φ:N1×[0,1]→R3 such that Φ(N1,τ) is homeomorphic to N1 for all τ in [0,1] and
Φ(⋅,0)=Identity,Φ(N1,1)=N2,Φ(x1([0,L]),1)=x2([0,L]). |
Following Gonzalez et al. [29], we define the minimal global radius of curvature of a closed curve x∈W1,p((0,L);R3), with p>1, by
Δ(x):=infs≠σ≠τ≠s∈[0,L)R(x(s),x(σ),x(τ)) |
where R(x,y,z) denotes the radius of the unique circle containing x≠y≠z≠x, with the convention R(x,y,z)=+∞ if x,y,z are collinear. The global radius of curvature determines the self-intersections of the tubular neighborhoods of a curve. More precisely, for every r>0 we define the r-tubular neighborhood of x by
Ur(x)=⋃s∈[0,L]Br(x(s)). |
Accordingly to Ciarlet et al. [12] we say that Ur(x) is not self-intersecting if for any p∈∂Ur(x) there exists a unique s∈[0,L] such that ‖p−x(s)‖=r. It turns out (see Gonzalez et al. [29]) that Δ(x)≥r if and only if Ur(x) is not self-intersecting. In particular, if Δ(x)>0 then x is simple, that is x:[0,L)→R3 is injective.
First we recall that a set \(S \subset {\mathbb R}^3\) is said to be \(2\)-rectifiable if it can be covered, up to an \({\mathcal H}^2\)-negligible set, by countably many 2-dimensional submanifolds of class C1, see [35,Chapter 3]. Given a 2-rectifiable set S, we denote by TxS the approximate tangent space of S⊂R3 at x, which exists for H2-almost every point x∈S [35,Chapter 3]. We also denote by G the Grassmannian of unoriented 2-dimensional planes in R3. The anisotropic integrands considered in the rest of the note will be continuous maps
F:R3×G∋(x,π)↦F(x,π)∈(0,+∞), |
verifying the lower and upper bounds
0<λ≤F(x,π)≤Λ,∀(x,π)∈R3×G. | (3.1) |
We also require that F is elliptic [25,5.1.2-5.1.5], namely its even and positively 1-homogeneous extension to R3×(Λ2(R3)∖{0}) is C2 and it is convex in the π variable. Given a 2-rectifiable set S⊂R3 we define:
F(S):=∫SF(x,TxS)dH2(x). | (3.2) |
Next, we need to define the spanning condition. For any closed set H⊂R3, let C(H) be the class of all smooth embeddings γ:S1→R3∖H. We say that C⊂C(H) is closed by homotopy if for every γ∈C then ˜γ∈C for any ˜γ∈[γ]∈π1(R3∖H). We denote by P(H,C) the family of all 2-rectifiable relatively closed sets S⊂R3∖H such that
S∩γ(S1)≠∅,∀γ∈C. |
We recall the following result, see [14,Theorem 2.7]:
Theorem 3.1. The problem
min{F(S):S∈P(H,C)} |
has a solution S∈P(H,C) and the set S is an (F,0,∞)-minimal set in R3∖H in the sense of Almgren [4].
Let N∈N∖{0} and p∈(1,+∞). For every i=1,…,N, let Li>0 and xi0,ti0,di0∈R3 be such ti0⊥di0 and |ti0|=|di0|=1. Moreover let κi1,κi2,ωi∈Lp(0,Li) such that
wi1:=(κi1,κi2,ωi)∈Lp((0,Li);R3), |
wi:=(wi1,xi0,ti0,di0)∈Lp((0,Li);R3)×R3×R3×R3, |
and
w:=(w11,w2,…,wN)∈Lp((0,L1);R3)×N∏i=2((Lp((0,Li);R3)×R3×R3×R3)=:V. |
We endow V with the natural Lp-norm, that we denote by ‖⋅‖V. For any i=1,…,N and for any w∈V, we denote by xi[w]∈W2,p((0,Li);R3) and ti[w],di[w]∈W1,p((0,Li);R3) the unique solutions (as proved in [29,Lemma 6]) of the Cauchy problem
{xi[w]′(s)=ti[w](s)ti[w]′(s)=κi1(s)di[w](s)+κi2(s)ti[w](s)×di[w](s)di[w]′(s)=ωi(s)ti[w](s)×di[w](s)−κi1(s)ti[w](s)xi[w](0)=xi0ti[w](0)=ti0di[w](0)=di0. |
It is easy to see that ti[w](s)⊥di[w](s) and |ti[w](s)|=|di[w](s)|=1 and consequently that
(ti[w](s),di[w](s),ti[w](s)×di[w](s)) |
is an orthonormal frame in R3, for any s∈[0,Li] and for any i=1,…,N. Let η,ν>0 and consider Ai(s)⊂R2 be compact and simply connected such that
Bη(0)⊂Ai(s)⊂Bν(0),∀s∈[0,Li],i=1,…,N. |
For any i=1,…,N we define
Ωi:={(s,ζ1,ζ2)∈R3:s∈[0,Li],(ζ1,ζ2)∈Ai(s)}, |
Λi[w]:={xi[w](s)+ζ1di[w](s)+ζ2ti[w](s)×di[w](s):(s,ζ1,ζ2)∈Ωi}, | (4.1) |
and
Λ[w]:=N⋃i=1Λi[w]. |
The system of closed rods is subjected to some constraints on w, enumerated below, which will identify the admissible subset W⊂V: First of all we assume that the midlines are closed and sufficiently smooth, that is
(C1) xi[w](Li)=xi[w](0)=xi0, for any i=1,…,N
and
(C2) ti[w](Li)=ti[w](0)=ti0, for any i=1,…,N.
To prescribe how many times the ends of the rods are twisted before being glued together, we prescribe the linking number between the midline and a closed curve close to the midline. More precisely, for any i=1,…,N we close up the curve xi[w]+τdi[w], for τ>0 fixed and small enough, defining as in Schuricht [34]
˜xiτ[w](s):={xi[w](s)+τdi[w](s)if s∈[0,Li]xi[w](Li)+τ(cos(φi(s−Li))di[w](Li)+sin(φi(s−Li))ti[w](Li)×di[w](Li))if s∈[Li,Li+1] | (4.2) |
where φi∈[0,2π) is the unique angle between di0 and di[w](Li) such that φi−π has the same sign as di0×di[w](Li)⋅ti0. We trivially identify xi[w] with its extension xi[w](s)=xi(Li) for any s∈[Li,Li+1] and therefore we require that for any i=1,…,N there is some li∈Z such that
(C3) Link(xi[w],˜xiτ[w])=li.
To encode the knot type of the midlines, for any i=1,…,N we fix a continuous mapping ℓi:[0,Li]→R3 such that ℓi(Li)=ℓi(0) and we require that
(C4) xi[w]≃ℓi.
Finally, in order to prevent the interpenetration of matter, following Ciarlet et al. [12] we require that
(C5)
∫Ωi(1−ζ1ki2(s)+ζ2ki1(s))dsdζ1dζ2≤|Λi[w]|∀i=1,…,N,andN⋂i=1int(Λi[w])=∅. |
We now require that our system of rods has a prescribed chain structure. We fix Lij∈ZN×N, with the property that |Li(i+1)|=1 for every i=1,…,N−1 and we assume that:
(C6) Link(xi[w],xj[w])=Lij.
We finally denote by W the set of all constraints, namely
W:={w∈V:(C1)–(C6) hold true}. |
It turns out that W is weakly closed in V (see Gonzalez et al. [29] and Schuricht [34]).
In what follows we will prescribe an elastic energy of the system of rods, which is a proper function
Eel:W→R∪{+∞},satisfyingEel(w)≥c‖w‖V, | (4.3) |
for some c>0. The second energy contribution we want to take into account is the weight of the rods. Let ρi∈L∞(Ωi) with ρ≥0 be the mass density functions and g be the gravitational acceleration. Let us define Eg:W→R∪{+∞} as
Eg(w):=N∑i=0∫Ωiρi(s,ζ1,ζ2)g⋅(xi[w](s)+ζ1di[w](s)+ζ2ti[w](s)×di[w](s))dsdζ1dζ2. |
The last contribution we want to take into account is the surface energy. Let Cw⊂C(Λ[w]) be the class of all γ∈C(Λ[w]) such that there exists i=1,…,N with
|Link(γ,xi[w])|=1,Link(γ,xj[w])=0,∀j≠i. |
Cw is closed by homotopy, see [30]. We define Esf:W→R∪{+∞} as
Esf(w):=inf{F(S):S∈P(Λ[w],Cw)}. |
We define the energy functional of our variational problem as
E(w):=Eel(w)+Eg(w)+Esf(w),w∈W. | (4.4) |
The first main result of the paper is given by the following existence theorem.
Theorem 4.1. Let ¯Eel be the lower semicontinuous envelope of Eel with respect to the weak topology of V. Assume that infWE<+∞. Then the problem
minw∈W¯Eel(w)+Eg(w)+Esf(w) |
has a solution w0∈W and there exists S∞∈P(Λ[w0],Cw0) which is an (F,0,∞)-minimal set in R3∖Λ[w0] in the sense of Almgren such that
Eel(w0)+Eg(w0)+F(S∞)=minw∈W¯Eel(w)+Eg(w)+Esf(w)=infw∈WE(w). |
First of all we prove that the weight and the soap film energy are weakly continuous.
Lemma 4.2. The functional Eg is weakly continuous on W.
Proof. Let (wh) be a sequence in W with wh⇀w in W for some w∈W. Then xi[wh]⇀xi[w] in W2,p and ti[wh]⇀ti[w], di[wh]⇀di[w] in W1,p. Then by Sobolev embedding we deduce that xi[wh]→xi[w] in C1,α and ti[wh]→ti[w], di[wh]→di[w] in C0,α for some α∈(0,1). This is enough to pass to the limit under the integral and get the claim.
The continuity of the soap film energy follows from the next theorem.
Theorem 4.3. Let (wh) be a sequence in W with wh⇀w in W for some w∈W. Assume that
(a) Sh∈P(Λ[wh],Cwh), for every h∈N;
(b) suph∈NF(Sh)=suph∈Ninf{F(S):S∈P(Λ[wh],Cwh)}<+∞.
Let μh:=FH2⌊Sh. Then the following three statements hold true:
μh⇀∗μ(up to subsequences), | (4.5) |
μ≥FH2⌊S∞,where S∞=(sptμ)∖Λ[w] is 2−rectifiable; | (4.6) |
S∞∈P(Λ[w],Cw). | (4.7) |
Proof. We first observe that the classes P(Λ[wh],Cwh) and P(Λ[w],Cw) are good classes in the sense of De Lellis et al. [14,Def. 2.2], as proved in [14,Thm. 2.7(a)]. Then the proof of (4.5) and (4.6) follows verbatim the proof of Theorem 2.5 of [14]. It is sufficient to observe that the convergence of {Λ[wh]} ensures that, whenever x∈S∞, we have d(x,Λ[wh])>0 for h large enough. We are left to prove (4.7), namely that S∞∩γ(S1)≠∅ for any γ∈Cw. Assume by contradiction that there exists γ∈Cw with S∞∩γ(S1)=∅. Since γ is compact and contained in R3∖Λ[w] and S∞ is relatively closed in R3∖Λ[w], there exists a positive ε such that the tubular neighborhood U2ε(γ) does not intersect S∞ and is contained in R3∖Λ[w]. Hence μ(U2ε(γ))=0, and consequently
limh→∞H2(Sh∩Uε(γ))=0. | (4.8) |
Denote by Bε the open disk of R2 with radius ε and centered at the origin of R2, and consider a diffeomorphism Φ:S1×Bε→Uε(γ) such that Φ|S1×{0}=γ. Let y belong to Bε and set γy:=Φ|S1×{y}. Then γy in [γ] represents an element of π1(R3∖Λ[w]). Since wh⇀w in W then (xi[wh]) converges to xi[w] strongly in W1,p((0,L);R3) for every i=1,…,N. In particular, (xi[wh]) converges to xi[w] uniformly on [0,Li] for every i=1,…,N, which implies the existence of δ>0 such that, for h sufficiently large, Λ[wh] is contained in Uδ(Λ[w]) with Uδ(Λ[w])∩Uε(γ)=∅. Hence, for such h and ε it follows that, for any y∈Bε, γy(S1)⊂R3∖Uδ(Λ[w]). This implies that ‖xi[wh]−γy‖∞≥δ for any y∈Bε and for every i=1,…,N. This estimate, together with the W1,p convergence of xi[wh] to xi[w], implies that
limh→+∞Link(xi[wh],γy)=Link(xi[w],γy),∀y∈Bε,∀i=1,…,N. |
As a consequence, for h large enough, γy∈Cwh which, combined with Sh∈P(Λ[wh],Cwh), yields Sh∩γy(S1)≠∅. Take now ˜π:S1×Bε→Bε as the projection on the second factor and let ˆπ:=˜π∘Φ−1. Then, ˆπ is Lipschitz-continuous and Bε is contained in ˆπ(Sh∩Uε(γ)), which entails that
πε2=H2(Bε)≤H2(ˆπ(Sh∩Uε(γ))≤(Lipˆπ)2H2(Sh∩Uε(γ)). |
We thus conclude that
H2(Sh∩Uε(γ))≥πε2(Lipˆπ)2 |
which contradicts (4.8).
Proof of Theorem 4.1. Thanks to the weak continuity of Eg and Esf, proved in Lemma 4.2 and Theorem 4.3, we deduce that ¯Eel(w)+Eg(w)+Esf(w) is the lower semicontinuous envelope of E, from which we get
infw∈W¯Eel(w)+Eg(w)+Esf(w)=infw∈WE(w). |
Let {wh} be a minimizing sequence for Eel+Eg+Esf. Since infWE<+∞ we can say that E(wh)≤c for some c>0. In particular, Eel(wh)≤c and, by coercivity of Eel, we have wh⇀w0 in W. We deduce, using again Lemma 4.2 and Theorem 4.3, that
¯Eel(w0)+Eg(w0)+Esf(w0)≤lim infh¯Eel(wh)+Eg(wh)+Esf(wh)≤lim infhE(wh)=infWE=infW¯Eel+Eg+Esf. |
Moreover, since Esf(w0)<+∞, applying Theorem 2.7 of [14] we deduce the claim.
The second main result of the paper concerns the dimensional reduction. In this section we consider cross sections with vanishing diameter. The set of constraints is almost the same, but in order to prevent the non-selfintersection in the limit configurations (otherwise the knot-type is not well defined) we replace the constraint (C5) by (C5)'. Precisely, we require that:
(C5)′ Δ(xi[w])≥Δ0 for some prescribed Δ0>0.
We denote by W′ the set of all constraints, namely
W′:={w∈Lp([0,L];R3):(C1)-(C2)-(C3)-(C4)-(C5)'-(C6) hold true}. |
It turns out that W′ is weakly closed in V (see again Gonzalez et al. [29] and Schuricht [34]). For every i=1,…,N, every ε>0 small enough and every w∈W′ we let
Λiε[w]:={xi[w](s)+ζ1di[w](s)+ζ2ti[w](s)×di[w](s):(s,ζ1,ζ2)∈Ωiε} | (5.1) |
where
Ωiε:={(s,ζ1,ζ2)∈R3:s∈[0,Li],(ζ1,ζ2)∈εAi(s)}. |
We also let
Λε[w]:=N⋃i=1Λiε[w]. |
The main goal of this section is to prove that as ε approaches 0, we recover by Γ-convergence the anisotropic Plateau problem with elastic one dimensional boundary. The first two energy contributions to take into account are the elastic energy Eel as in (4.3) and the scaled weight
Egε(w):=1ε2N∑i=1∫Ωiερi(s,ζ1,ζ2)g⋅(xi[w](s)+ζ1di[w](s)+ζ2ti[w](s)×di[w](s))dsdζ1dζ2 |
where ρi∈L∞(Ωi1) and ρi≥0. Concerning the soap film energy, similarly to the previous section, we define Cε,w⊂C(Λε[w]) as the class of all γ∈C(Λε[w]) such that there exists i=1,…,N with
|Link(γ,xi[w])|=1,Link(γ,xj[w])=0∀j≠i. |
We define Esfε:W′→R∪{+∞} as
Esfε(w):=inf{F(S):S∈P(Λε[w],Cε,w)}. |
Finally, Eε:W′→R∪{+∞} is given by
Eε(w):=Eel(w)+Egε(w)+Esfε(w). |
We now describe the Γ-limit functional. For any i=1,…,N, let ρi0:[0,Li]→R be given by
ρi0(s):=lim(ξ1,ξ2)→(0,0)ρi(s,ξ1,ξ2) |
and let
E0(w):=¯Eel(w)+N∑i=1∫Li0|Ai(s)|ρi0(s)g⋅xi[w](s)ds+inf{F(S):S∈P(Hw,Cw)}, |
where
Hw:=N⋃i=1xi[w]([0,Li]), |
and Cw is the class of all γ∈C(Hw) such that there exists i=1,…,N with
|Link(γ,xi[w])|=1,Link(γ,xj[w])=0∀j≠i. |
We are ready to state our second main result.
Theorem 5.1. Let (εh) be a positive and infinitesimal sequence and let (wh) be a sequence in W′ with suph∈NEεh(wh)≤c for some c>0. Then, up to a subsequence, wh⇀w in V and w∈W′. Moreover, the family {Eε}ε>0 Γ-converges to E0 as ε→0+ with respect to the weak topology of V, namely:
(a) for any sequence (εh) with εh→0, for any w∈W′ and for any sequence (wh) in W′ with wh⇀w in V we have
E0(w)≤lim infh→+∞Eεh(wh); | (5.2) |
(b) for any w∈W′ there is a sequence (εh) with εh→0 and a sequence (ˉwh) in W′ with ˉwh⇀w in V such that
E0(w)≥lim suph→+∞Eεh(ˉwh). | (5.3) |
As a standard consequence of Theorem 5.1 we have the next result.
Corollary 5.2. Let (εh) be a positive and infinitesimal sequence. For any h∈N and for any σh→0 let wh∈W′ be such that
Eεh(wh)≤infW′Eεh+σh. | (5.4) |
Then up to a subsequence wh⇀w0 in V and
E0(w0)=minW′E0. |
Here we give some preliminary propositions and then prove Theorem 5.1. Fix a positive and infinitesimal sequence (εh).
Proposition 5.3. Let (wh) be a sequence in W′ with suph∈NEεh(wh)≤c for some c>0. Then, up to a subsequence, wh⇀w in V and w∈W′.
Proof. The conclusion follows from the coercivity of Eel.
The study of the weight term is easy, since the weak convergence wh⇀w implies the uniform convergence of the midlines.
Proposition 5.4. For any w∈W′ and for any sequence (wh) in W′ with wh⇀w in V we have
limh→+∞Egεh(wh)=N∑i=1∫Li0|Ai(s)|ρi0(s)g⋅xi[w](s)ds. | (5.5) |
Proof. By the change of variables ζj=εhηj, j=1,2, we obtain that for any i=1,…N,
=1ε2h∫Ωiεhρi(s,ζ1,ζ2)g⋅(xi[wh](s)+ζ1di[wh](s)+ζ2ti[wh](s)×di[wh](s))dsdζ1dζ2=∫Ωi1ρi(s,εhη1,εhη2)g⋅(xi[wh](s)+εhη1di[wh](s)+εhη2ti[wh](s)×diwh](s))dsdη1dη2. |
Passing to the limit as h→+∞, using the fact that xi[wh]→xi[w] uniformly on [0,Li] for any i=1,…,N and applying the Dominated Convergence Theorem we conclude.
Now we pass to the limit in the soap film part of the energy. First of all we need the following Theorem whose proof requires minor modifications of the proof of Theorem 4.3.
Theorem 5.5. Let (wh) be a sequence in W′ with wh⇀w in W′ for some w∈W′. Assume that
(a) ∀h∈N,Sh∈P(Λεh[wh],Cεh,wh);
(b) suph∈NF(Sh)=suph∈Ninf{F(S):S∈P(Λεh[wh],Cεh,wh)}<+∞.
Let μh:=FH2⌊Sh. Then the following three statements hold true:
μh⇀∗μ(up to subsequences), | (5.6) |
μ≥FH2⌊S∞,where S∞=(sptμ)∖ Hw is 2−rectifiable, | (5.7) |
S∞∈P(Hw,Cw). | (5.8) |
Now we prove the existence of a recovery sequence.
Proposition 5.6. Consider w∈W′ and (wk)⊂W′ such that wk⇀w in W′. There exists (wkh) subsequence of (wk) such that
inf{F(S):S∈P(Hw,Cw)}≥lim suph→+∞Esfεh(wkh). | (5.9) |
Proof. Since wk⇀w in W′, xi[wh]→xi[w] uniformly on [0,Li] for any i=1,…,N. Then for every h∈N there exists kh∈N such that
‖xi[wkh]−xi[w]‖∞≤εh2,∀i=1,…,N. | (5.10) |
Since we can assume without loss of generality that
inf{F(S):S∈P(Hw,Cw)}<+∞, |
again applying Theorem 2.7 of [14], we find S∞∈P(Hw,Cw) such that
F(S∞)=min{F(S):S∈P(Hw,Cw)}. |
Now we set
Sh:=S∞∖Λεh[wkh]. |
For any γ∈C(Λεh[wkh]) not homotopic to a point in R3∖Λεh[wkh] we have
(S∞∖Λεh[wkh])∩γ(S1)≠∅. |
As a consequence,
lim suph→+∞Esfεh(w)≤lim suph→+∞F(Sh)≤F(S∞)=min{F(S):S∈P(Hw,Cw)}, |
which concludes the proof.
Proof of Theorem 5.1. The compactness statement is Proposition 5.3. Inequality (5.2) follows combining (5.5) and (5.7) with the subadditivity of the liminf operator. Next, for any w∈W′, we consider the constant sequence wh≡w. Applying Proposition 5.6, for every εh→0, the (unique) subsequence ˉwh≡w of (wh) satisfies obviously ˉwh⇀w in V and (5.9). Inequality (5.3) follows easily combining (5.5) and (5.9) with the superadditivity of the limsup operator.
A. D. R. has been partially supported by the NSF DMS Grant No. 1906451 and the NSF DMS Grant No. 2112311. L. L. has been supported by Gruppo Nazionale per l'Analisi Matematica, la Probabilità e le loro Applicazioni of the Istituto Nazionale di Alta Matematica (INdAM, Roma, Italy), project GNAMPA 2019 Analysis and optimization of thin structures.
The authors declare no conflict of interest.
[1] | W. K. Allard, A characterization of the area integrand, Symp. Math., XIV (1974), 429–444. |
[2] |
W. K. Allard, An a priori estimate for the oscillation of the normal to a hypersurface whose first and second variation with respect to an elliptic integrand is controlled, Invent. Math., 73 (1983), 287–331. doi: 10.1007/BF01394028
![]() |
[3] |
F. J. Jr. Almgren, Existence and regularity almost everywhere of solutions to elliptic variational problems among surfaces of varying topological type and singularity structure, Ann. Math., 87 (1968), 321–391. doi: 10.2307/1970587
![]() |
[4] | F. J. Jr. Almgren, Existence and regularity almost everywhere of solutions to elliptic variational problems with constraints, Mem. Am. Math. Soc., 4 (1976), 165. |
[5] | S. S. Antman, Nonlinear problems of elasticity, In: Applied Mathematical Sciences, 2 Eds., New York: Springer, 2005. |
[6] |
F. Bernatzki, R. Ye, Minimal surfaces with an elastic boundary, Ann. Global Anal. Geom., 19 (2001), 1–9. doi: 10.1023/A:1006734619701
![]() |
[7] | G. Bevilacqua, L. Lussardi, A. Marzocchi, Soap film spanning an elastic link, Quart. Appl. Math., 77 (2019), 507–523. |
[8] |
G. Bevilacqua, L. Lussardi, A. Marzocchi, Dimensional reduction of the Kirchhoff-Plateau problem, J. Elasticity, 140 (2020), 135–148. doi: 10.1007/s10659-020-09763-y
![]() |
[9] | A. Biria, E. Fried, Buckling of a soap film spanning a flexible loop resistant to bending and twisting, Proc. R. Soc. Lond. A, 470 (2014), 20140368. |
[10] |
A. Biria, E. Fried, Theoretical and experimental study of the stability of a soap film spanning a flexible loop, Int. J. Eng. Sci., 94 (2015), 86–102. doi: 10.1016/j.ijengsci.2015.05.002
![]() |
[11] |
Y. C. Chen, E. Fried, Stability and bifurcation of a soap film spanning a flexible loop, J. Elasticity, 116 (2014), 75–100. doi: 10.1007/s10659-013-9458-x
![]() |
[12] |
P. G. Ciarlet, J. Nečas, Injectivity and self-contact in nonlinear elasticity, Arch. Rational Mech. Anal., 97 (1987), 171–188. doi: 10.1007/BF00250807
![]() |
[13] |
R. Courant, The existence of minimal surfaces of given topological structure under prescribed boundary conditions, Acta Math., 72 (1940), 51–98. doi: 10.1007/BF02546328
![]() |
[14] | C. De Lellis, A. De Rosa, F. Ghiraldin, A direct approach to the anisotropic Plateau's problem, Adv. Calc. Var., 12 (2017), 211–223. |
[15] | C. De Lellis, F. Ghiraldin, F. Maggi, A direct approach to Plateau's problem, JEMS, 9 (2017), 2219–2240. |
[16] |
G. De Philippis, A. De Rosa, F. Ghiraldin, A direct approach to Plateau's problem in any codimension, Adv. Math., 288 (2016), 59–80. doi: 10.1016/j.aim.2015.10.007
![]() |
[17] |
G. De Philippis, A. De Rosa, F. Ghiraldin, Rectifiability of varifolds with locally bounded first variation with respect to anisotropic surface energies, Commun. Pure Appl. Math., 71 (2018), 1123–1148. doi: 10.1002/cpa.21713
![]() |
[18] |
G. De Philippis, A. De Rosa, F. Ghiraldin, Existence results for minimizers of parametric elliptic functionals, J. Geom. Anal., 30 (2020), 1450–1465. doi: 10.1007/s12220-019-00165-8
![]() |
[19] |
G. De Philippis, A. De Rosa, J. Hirsch, The Area Blow Up set for bounded mean curvature submanifolds with respect to elliptic surface energy functionals, Discrete Cont. Dyn. A, 39 (2019), 7031–7056. doi: 10.3934/dcds.2019243
![]() |
[20] |
A. De Rosa, Minimization of anisotropic energies in classes of rectifiable varifolds, SIAM J. Math. Anal., 50 (2018), 162–181. doi: 10.1137/17M1112479
![]() |
[21] |
A. De Rosa, S. Gioffrè, Quantitative stability for anisotropic nearly umbilical hypersurfaces, J. Geom. Anal., 29 (2019), 2318–2346. doi: 10.1007/s12220-018-0079-2
![]() |
[22] | A. De Rosa, S. Gioffrè, Absence of bubbling phenomena for non convex anisotropic nearly umbilical and quasi Einstein hypersurfaces, 2018, arXiv: 1803.09118. |
[23] |
A. De Rosa, S. Kolasinski, Equivalence of the ellipticity conditions for geometric variational problems, Commun. Pure Appl. Math., 73 (2020), 2473–2515. doi: 10.1002/cpa.21890
![]() |
[24] |
A. De Rosa, S. Kolasinski, M. Santilli, Uniqueness of critical points of the anisotropic isoperimetric problem for finite perimeter sets, Arch. Rational Mech. Anal., 238 (2020), 1157–1198. doi: 10.1007/s00205-020-01562-y
![]() |
[25] | H. Federer, Geometric measure theory, New York: Springer-Verlag New York Inc., 1969. |
[26] | L. Giomi, L. Mahadevan, Minimal surfaces bounded by elastic lines, Proc. R. Soc. Lond. A, 468 (2012), 1851–1864. |
[27] |
G. G. Giusteri, P. Franceschini, E. Fried, Instability paths in the Kirchhoff–Plateau problem, J. Nonlinear Sci., 26 (2016), 1097–1132. doi: 10.1007/s00332-016-9299-4
![]() |
[28] |
G. G. Giusteri, L. Lussardi, E. Fried, Solution of the Kirchhoff-Plateau problem, J. Nonlinear Sci., 27 (2017), 1043–1063. doi: 10.1007/s00332-017-9359-4
![]() |
[29] |
O. Gonzalez, J. H. Maddocks, F. Schuricht, H. von der Mosel, Global curvature and self-contact of nonlinearly elastic curves and rods, Calc. Var., 14 (2002), 29–68. doi: 10.1007/s005260100089
![]() |
[30] |
J. Harrison, H. Pugh, Existence and soap film regularity of solutions to Plateau's problem, Adv. Calc. Var., 9 (2016), 357–394. doi: 10.1515/acv-2015-0023
![]() |
[31] |
J. Harrison, H. Pugh, General Methods of Elliptic Minimization, Calc. Var., 56 (2017), 123. doi: 10.1007/s00526-017-1217-6
![]() |
[32] | T. M. Hoang, E. Fried, Influence of a spanning liquid film on the stability and buckling of a circular loop with intrinsic curvature or intrinsic twist density, Math. Mech. Solids, 23 (2016), 43–66. |
[33] |
H. Lewy, On mimimal surfaces with partially free boundary, Commun. Pure Appl. Math., 4 (1951), 1–13. doi: 10.1002/cpa.3160040102
![]() |
[34] |
F. Schuricht, Global injectivity and topological constraints for spatial nonlinearly elastic rods, J. Nonlinear Sci., 12 (2002), 423–444. doi: 10.1007/s00332-002-0462-8
![]() |
[35] | L. Simon, Lectures on geometric measure theory, Canberra: Australian National University, Centre for Mathematical Analysis, 1983. |
[36] |
J. E. Taylor, Crystalline variational problems, Bull. Am. Math. Soc., 84 (1978), 568–588. doi: 10.1090/S0002-9904-1978-14499-1
![]() |
1. | Dario Mazzoleni, Benedetta Pellacci, Calculus of variations and nonlinear analysis: advances and applications, 2023, 5, 2640-3501, 1, 10.3934/mine.2023059 | |
2. | Chiara Lonati, Alfredo Marzocchi, On self-contact and non-interpenetration of elastic rods, 2024, 1081-2865, 10.1177/10812865231226311 | |
3. | Guido De Philippis, Antonio De Rosa, The anisotropic min‐max theory: Existence of anisotropic minimal and CMC surfaces, 2024, 77, 0010-3640, 3184, 10.1002/cpa.22189 | |
4. | Antonio De Rosa, Reinaldo Resende, Boundary regularity for anisotropic minimal Lipschitz graphs, 2024, 49, 0360-5302, 15, 10.1080/03605302.2023.2294335 | |
5. | Giulia Bevilacqua, Luca Lussardi, Alfredo Marzocchi, Soap films: from the Plateau problem to deformable boundaries, 2024, 15, 2038-0909, 137, 10.2478/caim-2024-0019 | |
6. | Chiara Lonati, Alfredo Marzocchi, Self-Contact Preventing Energy for Tubular Rods, 2025, 35, 0938-8974, 10.1007/s00332-025-10128-9 |