Solution of the Kirchhoff–Plateau Problem

The Kirchhoff–Plateau problem concerns the equilibrium shapes of a system in which a flexible filament in the form of a closed loop is spanned by a liquid film, with the filament being modeled as a Kirchhoff rod and the action of the spanning surface being solely due to surface tension. We establish the existence of an equilibrium shape that minimizes the total energy of the system under the physical constraint of noninterpenetration of matter, but allowing for points on the surface of the bounding loop to come into contact. In our treatment, the bounding loop retains a finite cross-sectional thickness and a nonvanishing volume, while the liquid film is represented by a set with finite two-dimensional Hausdorff measure. Moreover, the region where the liquid film touches the surface of the bounding loop is not prescribed a priori. Our mathematical results substantiate the physical relevance of the chosen model. Indeed, no matter how strong is the competition between surface tension and the elastic response of the filament, the system is always able to adjust to achieve a configuration that complies with the physical constraints encountered in experiments.


Introduction
Liquid films spanning rigid frames have been of longstanding interest to physicists and mathematicians, thanks to the sheer beauty of the countless observable shapes. After the experimental investigations of Plateau (1849), anticipated by Lagrange's (1867) theoretical treatment of the minimal surface problem, the first satisfactory proofs of the existence of a surface of least area bounded by a fixed contour were provided only in the twentieth century by Garnier (1928), Radó (1930), and Douglas (1931). This formed a basis for a wealth of mathematical investigations regarding minimal surfaces, concerning various aspects and generalizations of the classical Plateau problem. The interested reader is referred to the treatises by Dierkes et al. (2010a, b) for a comprehensive review of the formative contributions.
An important generalization to the situation in which the boundary of the minimal surface is not fixed but is constrained to lie on a prescribed manifold was initially treated by Courant (1940) and Lewy (1951), whose work spurred a number of important mathematical contributions, as reviewed by Li (2015). An existence theorem for a complementary generalization, in which part of the boundary is fixed and the remaining part is an inextensible but flexible string, was later proved by Alt (1973). In the present article, we introduce a problem which combines those generalizations. We consider situations in which the boundary of the minimal surface lies on a deformable manifold, namely the surface of an elastic loop. The filament forming the loop is assumed to be thin enough to be modeled faithfully by a Kirchhoff rod, an unshearable inextensible rod which can sustain bending of its midline and twisting of its cross sections [see Antman (2005, Chapter VIII)]. This is a mathematically one-dimensional theory that describes a three-dimensional object, endowed with a nonvanishing volume, since the material cross sections have nonvanishing area. The Kirchhoff-Plateau problem concerns the equilibrium shapes of a system composed by a closed Kirchhoff rod spanned by an area-minimizing surface.
In recent years, some attention has been drawn to the Kirchhoff-Plateau problem following a paper by Giomi and Mahadevan (2012), and the stability properties of flat circular solutions have been investigated, under various conditions regarding the material properties of the rod, by Chen and Fried (2014), Fried (2014, 2015), Giusteri et al. (2016), and Hoang and Fried (2016). An existence result for a similar problem was given by Bernatzki and Ye (2001), but the elastic energy used therein fails to satisfy the basic physical requirement of invariance under superposed rigid transformations and thereby implicitly entails the appearance of unphysical forces. Moreover, a strong hypothesis is used to avoid the issue of self-contact.
Importantly, in all of these studies the boundary of the spanning surface is assumed to coincide with the rod midline and not to lie on the surface of the rod. This amounts to a slenderness assumption. Moreover, the surface is viewed as diffeomorphic to a disk, except by Bernatzki and Ye (2001), who employ the theory of currents. We relax both of these assumptions. Regarding the filament, we retain its three-dimensional nature for two physically well-justified reasons. First, there is a significant separation of scales between the typical thickness of the liquid film and the cross-sectional thickness of the filaments used in experimental investigations: a minimum of two orders of magnitude. Consequently, while the liquid film is still represented as two-dimensional, it seems appropriate to treat the bounding loop as a three-dimensional object. Second, the possibility of generating nontrivial shapes due to the interaction between the film and the bounding loop relies on the anisotropic material properties of the filament, which are often associated with how its cross sections are shaped.
Furthermore, the physical presence of the bounding loop requires a proper treatment of the constraint of noninterpenetration of matter, which is clearly at play in real experiments and even becomes essential, since the bounding loop can sustain large deflections but remains constrained when self-contact occurs. If, in particular, the relative strength of surface tension with respect to the elastic response of the filament becomes large, then the compliance of the mathematical solution with physical requirements can only be guaranteed by including the noninterpenetration constraint. To include all these properties in a variational framework, we base our treatment of rod elasticity on Schuricht's (2002) elegant approach. We introduce a minor simplification in the presentation of the constraint of local noninterpenetration of matter, obtaining equivalent results specialized to the case of Kirchhoff rods.
In developing our variational approach, the most delicate point involves the treatment of the spanning surface. The physical phenomenon at play is the minimization of the liquid surface area due to the presence of a homogeneous surface tension. Various mathematical models have been proposed with different characteristics. Representing the surface via a mapping from a manifold into the ambient space, despite being the first successful approach, poses severe and completely unphysical limitations on the topology of the surface. To cope with this issue, the theories of integral currents and of varifolds were applied to the Plateau problem by Federer and Fleming (1960) and Almgren (1968), respectively. However, their approaches also fail either to cover all the physical soap film solutions to the Plateau problem or to furnish a sufficiently general existence result. An alternate route was initiated by Reifenberg (1960), who treated the surface as a point set that minimizes the two-dimensional Hausdorff measure. This purely spatial point of view, adopted also by De Pauw (2009) and David (2014), deals nicely with the topology of solutions, but makes it difficult to handle a generic boundary condition. A more complete treatment, which covers all the soap film solutions to the Plateau problem, has been developed by Harrison (2014) and Harrison and Pugh (2016a) using differential chains. In that setting, it is possible to consider generic bounding curves due to the presence of a well-defined boundary operator.
In relation to the problem at hand, the approach to the Plateau problem due to De Lellis et al. (2014) proves to be superior in many respects. First and foremost, De Lellis et al. (2014) treat the surface as the support of a Radon measure, adopting a spatial point of view, and thereby obtain the optimal soap film regularity defined by Almgren (1976) and Taylor (1976). Moreover, their definition of the spanning conditions, built on ideas of Harrison (2014) [further developed by Harrison and Pugh (2016a, b)], allows for an apt treatment of the free-boundary problem, whereas all prior approaches become rather difficult to use when the boundary of the spanning surface is not prescribed. Finally, their strategy has the physically relevant property of being insensitive to changes in the topology of the spanning surface, which can easily be observed during the relaxation to equilibrium of the system with an elastic bounding loop. Another interesting approach that allows for treatment of the free-boundary problem is provided in a recent paper by Fang (2016), but its generalization to the case of a deformable bounding loop seems to require a more sophisticated apparatus than that employed in our treatment.
We present the formulation of the Kirchhoff-Plateau problem in Sect. 2, where we introduce the energy functionals pertaining to the parts of our system and a suitable expression of the spanning condition. Of particular importance is also the introduction of various physical and topological constraints regarding the noninterpenetration of matter and the link and knot type of the bounding loop. Indeed, those natural requirements necessitate a variational approach to the problem, since they entail an inherent lack of smoothness, leaving open the question of the validity of the Euler-Lagrange equations at energy-minimizing configurations. Furthermore, they hinder the convexity properties of a seemingly innocuous functional. This becomes even more evident when the surface energy is added to the picture, so that the presence of multiple stable and unstable equilibria cannot in general be precluded.
We prove our main existence result in Sect. 3, subsequent to establishing some preliminary facts. The essential feature of our treatment is a dimensional reduction which is performed by expressing the total energy of the system (bounding loop plus spanning surface) as a functional of the geometric descriptors of the bounding loop only. This is done by introducing a strongly nonlocal term entailing the minimization of the surface energy for a fixed shape of the bounding loop. Clearly, this step is justified by the existence of a solution to the problem of finding an area-minimizing surface spanning a three-dimensional bounding loop. We prove this (Theorem 3.2) as a direct application of a result established by De Lellis et al. (2014, Theorem 4). Subsequently, it is necessary to adapt the arguments of De Lellis et al. (2014) to our situation, in which the set that is spanned by the surface changes along minimizing sequences for the Kirchhoff-Plateau energy. This is accomplished in Lemmas 3.4 and 3.5, and based on these results, we establish the lower semicontinuity property that is needed to establish, in Theorem 3.6, the existence of a solution to the Kirchhoff-Plateau problem.

Formulation of the Problem
We seek to study the existence of a stable configuration of a liquid film that spans a flexible loop. Such a configuration is in reality metastable, since the liquid film will eventually break after becoming sufficiently thin, but we confine our attention to what happens before any such catastrophic event.
The flexibility of the bounding loop represents a major difference between our problem and the Plateau problem, in which the liquid film spans a fixed boundary. It is therefore essential to model the elastic behavior of the loop in response to deformations, requiring a physical description much more sophisticated than that sufficient for a fixed boundary. In particular, we consider a loop formed by a slender filament with a nonvanishing cross-sectional thickness and subject to the physical constraint of noninterpenetration of matter. As discussed below, it is still reasonable to approximate the liquid film by a surface, since its thickness is typically much smaller than the cross-sectional thickness of the bounding loop, but an appropriate definition of how the surface spans the bounding loop is necessary.
We next introduce the precise mathematical definitions needed to formulate the Kirchhoff-Plateau problem in a way that takes into account the physical requirements mentioned above. In so doing, we follow a variational approach by defining the energy of the different components of our system.

Preliminary Considerations and Assumptions
The main assumption we impose on the bounding loop is that its length is much larger than the characteristic thickness of its cross sections, which allows us to employ the theory of rods in our description. Within that theory, as presented for example by Antman (2005), a rod is fully described by a curve in the three-dimensional Euclidean space, called the midline, a family of two-dimensional sets, describing the material cross section at each point of the midline, and a family of material frames, encoding how the cross sections are "appended" to the midline. Such a family of material frames corresponds also to a curve in the group of rotations of the three-dimensional space.
We also assume that the rod is unshearable (namely that the material cross section at any point of the midline lies in the plane orthogonal to the midline at that point) and that its midline is inextensible. Together, these assumptions amount to choosing a Kirchhoff rod as a model for the filament from which the bounding loop is made. Under these assumptions, the shape of the loop is uniquely determined by assigning the shape of the cross sections and three scalar fields: two flexural densities κ 1 and κ 2 together with a twist density ω. From these fields we can reconstruct the midline x and a director field d, orthogonal to the tangent field t := x , that gives the material frame as {(t(s), d(s), t(s) × d(s)) : s ∈ [0, L]}, where s is the arc-length parameter and L is the total length of the midline. Indeed, once suitable conditions at s = 0 are assumed, the fields x and d are the unique solution of the system of ordinary differential equations (2.2) It is not a priori granted that the solution to (2.1)-(2.2) describes a closed loop. This property, being essential to the treatment of the Kirchhoff-Plateau problem, is imposed later as a constraint on the variational problem.
For the variational problem that we plan to investigate, it is convenient to assume that each of the densities κ 1 , κ 2 , and ω belongs to the Lebesgue space L p ([0, L]; R) for some p in (1, ∞). This, by a classical result of Carathéodory (see, for instance, Hartman 1982), ensures that (2.1)-(2.2) have a unique solution, with We further assume |t 0 | = |d 0 | = 1. On this basis, we can use the structure of (2.1) 2,3 to prove that |t(s)| = |d(s)| = 1 for every s in [0, L].
The material cross section at each s is given by a compact simply connected domain A(s) of R 2 such that the origin 0 2 of R 2 belongs to int(A(s)). A rod of finite cross-sectional thickness can then be described as the image in the three-dimensional Euclidean space of the set (2.3) By our assumptions, there exists an R > 0 such that |ζ 1 | < R and |ζ 2 | < R for any (s, ζ 1 , ζ 2 ) in . We remark that, for the rod model to be a faithful representation of the mechanics of the filament from which the bounding loop is made, it is necessary that the maximum thickness R be small compared to the length L of the loop. Once the family of material cross sections is assigned, any configuration of the rod in the (three-dimensional) ambient space corresponds to an element belonging to the Banach space Whereas all of the information regarding the shape of the rod is encoded in the flexural and twist densities, namely in the component w 1 = (κ 1 , κ 2 , ω) of w belonging to L p ([0, L]; R 3 ), the clamping parameters x 0 , t 0 , and d 0 determine how the rod is translated and rotated in the ambient space. We denote the midline, the director field, and the rod configuration computed from (2.1)-(2.2) and (2.3) for a given w in V as x

Individual Contributions to the Energy of the Bounding Loop
We are now positioned to define the energy of the rod as a functional on V . We consider three contributions entering this functional in an additive way: (1) the stored elastic energy, related to shape modifications; (2) the noninterpenetration constraint; and (3) the potential energy of an external load, such as the weight of the rod. The first term, being related only to shape modifications, can be expressed as the integral of an elastic energy density which depends only on s and w 1 . We then introduce f : R 3 × [0, L] → R ∪ {+∞} and define the shape energy of the bounding loop as and we assume that f (·, s) is continuous and convex for any s in [0, L], that f (a, ·) is measurable for any a in R 3 , and that f (a, s) is uniformly bounded below by a constant. These assumptions guarantee that E sh is weakly lower semicontinuous on the reflexive Banach space V -a key property in applying the direct method of the calculus of variations. To ensure the necessary coercivity, we also impose the natural growth condition To include the local noninterpenetration constraint, we use the characteristic function of the closed subset N of V containing those elements w such that for almost every s in [0, L]. We thus add to the energy of the loop the term Condition (2.4) can be derived by relaxing the standard requirements that p[w] be orientation preserving and locally injective [see Antman (2005, Chap. VIII.7); Schuricht 2002]. We will prove, in Theorem 2.3, that our penalization strategy (which differs slightly from that of Schuricht 2002) is sufficient to guarantee the local injectivity of p[w] on int( ) for configurations with finite energy. The concept of local interpenetration of matter is illustrated in Fig. 1.   Fig. 1 Local injectivity fails when adjacent material cross sections interpenetrate. This is illustrated in the encircled region (magnified on the right) where two material cross sections, traced on the surface, overlap, due to the excessive curvature of the midline. Global injectivity fails also when interpenetration occurs between cross sections that lie far apart along the midline. This situation, in which local injectivity is not hindered, is depicted in the framed region Finally, we account for the effects of the weight of the rod by considering the potential energy where ρ > 0 represents the mass density at each point of the rod and the vector g represents the constant gravitational acceleration. From the above definitions, we obtain a weakly lower semicontinuous functional representing the total energy of the bounding loop.

Energy of the Liquid Film
A fundamental tenet in physical chemistry is that the process of building an interface between two immiscible substances is accompanied by an energetic cost. Roughly speaking, each portion of each given substance prefers to be surrounded by the same substance and some interfacial energy is developed whenever this is not the case. For this reason, a droplet of water surrounded by air tends to assume a spherical shape: for a fixed droplet volume, it minimizes the area of the interface.
To produce a liquid film, it is necessary to counteract this tendency, since doing so entails stretching the droplet, necessarily increasing the area of the liquid/air interface. We operate in two ways. On the one hand, the solid bounding loop provides a third substance to which the liquid is attracted, since the energy density per unit area of the liquid/solid interface is lower than that associated with the liquid/air interface. On the other hand, it is possible to tamper a bit with the liquid to further reduce the energetic cost of the liquid/air interface.
The second objective is accomplished by adding a small amount of surfactant to the liquid. Since the energy density per unit area of the liquid/air interface is lower for a higher surfactant concentration, surfactant molecules migrate toward the interface, leaving water in the bulk. In the liquid film configuration this produces two leaflets of surfactant phase (that lower the interfacial energy) covering a thin water layer (that provides a significant cohesion to the structure). In this context, we can define the surface tension σ of the liquid as the energy density per unit area of the liquid/air interface. It is physically reasonable to assume that σ is a homogeneous positive quantity, representing the ratio between the total interfacial energy and the surface area of the interface.
To arrive at a mathematical model for the liquid film, we now combine a geometric approximation with the notion of interfacial energy discussed above. We assume that the thickness of the film (two surfactant leaflets plus the water layer) is negligible and we represent it as a two-dimensional object S, but we keep track of the fact that it is built with two surfactant leaflets. We thus define the energy of the liquid film as where H d represents the d-dimensional Hausdorff measure. Fig. 2 The thickness of the liquid film is orders of magnitude smaller than the cross-sectional thickness of the bounding loop. The white arrow points to the bright region where the thickness of the soap film slightly increases just before touching the bounding loop, realized using commercial fishing line with cross-sectional diameter of approximately 0.2 mm. Since most of the filament surface is not covered by the liquid film, the thickness of the latter must be considerably less than 0.2 mm Two remarks are in order. First, there is no evidence in (2.6) of the energy associated with the liquid/solid interface along which the film is in contact with the bounding loop. This is, in principle, a significant contribution, but the energy barrier that must be overcome to detach the film from the bounding loop is so high that its effect can be replaced by the spanning condition, discussed in the next section, designed to prevent detachment, encoding in essence the infinite height of the aforementioned barrier. The small corrections to that energy due to the size and shape of the liquid/solid interface are negligible as far as detachment is concerned. The influence of those corrections on the shape energy of the film is also negligible if the film is assumed to be of vanishing thickness. Indeed, they would influence the shape of the liquid/solid interface at length scales that are smaller than its typical thickness and, hence, not captured by our model.
Second, it is important to justify the vanishing-thickness approximation, especially in view of the nonvanishing thickness that we attribute to the rod modeling the bounding loop. The typical thickness of a soap film is on the submicron scale, while a bounding loop made of a strand of human hair would have a cross section with characteristic thickness of some tens of microns: this indicates that in many practical examples of liquid films bounded by flexible loops the characteristic thickness of the loop is at least two orders of magnitude greater than the thickness of the film (see Fig. 2).

The Spanning Condition
We now provide a precise mathematical formulation of the conditions stipulating that the liquid film spans the bounding loop. In so doing, we borrow an elegant idea, that exploits notions of algebraic topology, introduced by Harrison (2014) [and further developed by Harrison and Pugh (2016a, b)]. For our application it is convenient to present that idea in the form provided by De Lellis et al. (2014), specialized to the particular setting in which the ambient space is three-dimensional.
Definition 2.1 Let H be a closed set in R 3 and consider the family C H := γ : S 1 → R 3 \ H : γ is a smooth embedding of S 1 into R 3 . Fig. 3 An appropriate choice of homotopy classes determines which holes of a bounding loop with points of self-contact are covered. For the particular loop depicted here, which is subject to self-contact without interpenetration at the central crossing point (black cross), if we seek a spanning set relative to the homotopy class of the loop a or b, spanning surfaces that cover only the hole on the left or on the right will be allowed, respectively. If, instead, we consider the homotopy class of the loop c, both holes must be covered by the spanning set Then, any subset C of C H is closed by homotopy (with respect to H ) if C contains all the elementsγ of C H belonging to the same homotopy class [γ ] of any γ in C. Such homotopy classes are elements of the first fundamental group π 1 of R 3 \ H .
Given a subset C of C H closed by homotopy, a relatively closed subset The previous definition offers several advantages in comparison with more classical definitions. Most importantly, the set that is spanned by the surface representing the liquid film need not be a one-dimensional structure, and it can be any closed set in the ambient space. This allows for the possibility that the surface spans the bounding loop [w] of finite cross-sectional thickness. Another useful feature of this definition is that it allows for choice regarding the number of "holes to be covered" by the spanning set. Namely, when the set H has a somewhat complex topology, it is permissible to restrict attention to surfaces that span only significant subregions of H . To give a trivial example, illustrated in Fig. 3, where H is a loop subject to self-contact, it is possible to seek surfaces that span one, the other, or both holes. Notice that, as shown below, the noninterpenetration constraint still allows for points on the surface of the bounding loop to come into contact, since it entails only a superposition of points on the surface of the rod.
In particular, the subset D [w] of C [w] containing all γ that are not homotopic to a constant ([γ ] = 1 π 1 ) is closed by homotopy. We will then seek a surface that is a D [w] -spanning set of the bounding loop [w]. This is a maximal choice in the sense that we cannot include paths homotopic to a constant in the spanning condition. Indeed, since [w] is a compact set, it is easy to see that any {1 π 1 }-spanning set of [w] required to intersect all constant paths would fill all of R 3 \ [w]. This would certainly not represent the behavior of a liquid film.

The Kirchhoff-Plateau Problem
The basic step in connecting our mathematical model to the experimental observations consists in finding a minimizer w for the functional where w belongs to a suitable subset U of V that encodes the additional physical and topological constraints on the rod modeling the bounding loop. These constraints are: (1) the closure of the midline; (2) the global gluing conditions; (3) global noninterpenetration of matter; (4) the knot type of the midline.
When we combine the closure constraint with the definition of the Kirchhoff-Plateau functional E KP , a strong competition between the action of the spanning film and the elastic response of the bounding loop can arise. In a typical situation, the curvature of the loop tends to be minimized, producing somewhat wider configurations which, in turn, require spanning films with larger surface areas. The equilibrium shape is strongly influenced by the relative strength of surface tension with respect to the properties of the filament, but the proper inclusion of the global noninterpenetration constraint guarantees the physical relevance of the solutions the existence of which we establish.
We will now present, following Schuricht (2002), the precise formulation of the aforementioned constraints and a lemma, proving that the set U is weakly closed in V , as this is the essential property needed to establish the tractability of those constraints within our variational approach.
The fact that the midline is a closed curve can be readily expressed by (2.9)

x[w](L) = x[w](0)
The clamping conditions (2.8)-(2.9) are important in view of the preferred direction associated with the gravitational acceleration g. Indeed, the weight term breaks the invariance of our problem under rigid rotations. Moreover, (2.8)-(2.9) effectively describe the physical operation of holding the flexible structure at one point with tweezers. The right-hand sides of (2.8)-(2.9), namely x 0 , t 0 , and d 0 , are referred to as the clamping parameters.

Note that we do not require that d[w](L) = d[w](0), since the rod can be glued fixing an arbitrary angle between d[w](L) and d[w](0) (local gluing condition) while
simultaneously respecting the clamping conditions. To prescribe how many times the ends of the rod are twisted before being glued together, we define the integer link type of the closed rod as the linking number of the closed midline x[w] and the curve x[w]+ d [w]. Although a sufficiently small can always be chosen, the curve may need to be closed as indicated by Schuricht (2002, Sect. 4.1 (2.10) Note that although condition (2.10) implies the global injectivity of p[w] only on int( ), it is sufficient for our purposes, since we wish to allow for self-contact of the bounding loop, preventing only interpenetration. Clearly, global injectivity implies local injectivity, but the converse is not true. The local property is only a necessary condition for the global injectivity to hold. Indeed, as depicted in Fig. 1, global injectivity may fail due to the overlapping of cross sections belonging to regions of the rod that lie far apart along the midline and in which local injectivity holds true.
To encode the knot type of the midline, we invoke the notion of isotopy class for closed curves.
By means of isotopy classes we can encode the knot type of the bounding loop as follows. We fix a continuous mapping : [0, L] → R 3 such that (L) = (0) and we say that an element w of V , for which the closure and clamping conditions (2.8)-(2.9) hold, has the knot type of if x[w] . (2.11) We can now prove a theorem about the noninterpenetration constraint and our basic lemma about the weak closure of the set U of competitors for our minimization problem.

. This means that t[w](s) = t[w](a) for every s in
which is (2.16).
We have therefore established the local injectivity of p[w] and moreover that, being continuous, p [w] is an open mapping on int( ). Given the local injectivity, the global injectivity follows from condition (2.10), as proved by Schuricht (2002, Theorem 3.8). Proof If U = ∅ the assertion is true. If U = ∅, then the proof is a mere restatement of those of Lemma 3.9, Lemma 4.5, and Lemma 4.6 of Schuricht (2002).
We stress that the global injectivity condition (2.10) is crucial in ensuring the closure proved in Lemma 2.4. Indeed, just as the physical noninterpenetration of matter makes it impossible to change the type of a knot without tearing the loop that forms it, global injectivity entails that the knot type is preserved when passing to the limit in a sequence within the constrained set U . In particular, the fact that constrained sets defined by different knot types are well separated in V stems from the nonvanishing thickness of the cross sections and condition (2.10) and would be inevitably lost in the vanishingthickness limit, as illustrated in Fig. 4.  Fig. 4 A nonvanishing cross-sectional thickness is essential for distinguishing knot types in the presence of self-contact. As a trivial example, the trefoil knot (left) and the unknot (right) are clearly distinct when the cross sections have a nonvanishing thickness, even in the presence of self-contact. On the contrary, in the vanishing-thickness limit (center), reaching self-contact from each of the two knotted configurations produces topologically equivalent structures and the distinction between a knot and an unknot is lost E KP , under appropriate physical and topological conditions. As necessary steps toward this goal, we first prove the existence of an energy-minimizing configuration for the bounding loop in the absence of the liquid film and then we demonstrate the existence of an area-minimizing spanning surface for a rigid bounding loop. This constitutes a somewhat more physical version of the Plateau problem, since the bounding loop is treated as a three-dimensional object.
Those two results are quite straightforward, given the results of the previous section and the general theorem proved by De Lellis et al. (2014, Theorem 4). On the other hand, the proof of our main existence result requires establishing two lemmas in which the arguments of De Lellis et al. (2014) are adapted to our situation, wherein the set that is spanned by the surface changes along minimizing sequences for the Kirchhoff-Plateau energy. Based on these results, we are able to verify the lower semicontinuity property that is needed to establish the existence of a solution of the Kirchhoff-Plateau problem. Moreover, Lemma 3.4 provides an interesting extension of the main compactness result of De Lellis et al. (2014), which could be used also in other contexts when studying the convergence of minimal surfaces induced by the convergence of the structure spanned by the surfaces. Proof Consider a minimizing sequence {w k } such that E loop (w k ) < M for some M in R, with w k belonging to the subset U of the reflexive Banach space V introduced in Lemma 2.4. The existence of a competitorw guarantees that U is not empty. By the coercivity of E loop , U is bounded in V . Hence, it is possible to extract a weakly converging subsequence w k i w. Since, by Lemma 2.4, U is weakly closed, w must belong to U and, by the weak lower semicontinuity of E loop , we conclude that  D [w] defined as in Sect. 2.3. Moreover, for any γ : S 1 → R 3 and for any r > 0 we denote by U r (γ ) the tubular neighborhood of radius r around γ . Lemma 3.4 Consider a sequence of closed nonempty subsets k of R 3 converging in the Hausdorff topology to the closed set = ∅. Assume that, for every k in N, we have countably H 2 -rectifiable sets S k belonging to P( k ), where P( k ) is a good class in the sense of De Lellis et al. (2014, Definition 1), and such that Then the measures μ k := H 2 S k constitute a bounded sequence, μ k * μ up to the extraction of a subsequence, and the limit measure satisfies the open disk of R 2 with radius ε and centered at the origin of R 2 , and consider a diffeomorphism : Let y belong to B ε (0 2 ) and set γ y := | S 1 ×{ y} . Then γ y in [γ ] represents an element of π 1 (R 3 \ [w]). Let x k and x denote the midlines corresponding, respectively, to w k and w. Since {w k } converges weakly to w in U , {x k } converges to x strongly in W 1, p ([0, L]; R 3 ). In particular, {x k } converges to x uniformly on [0, L], which implies that, for k suffi- Hence, for such k and ε it follows that, for any y in B ε (0 2 ), γ y belongs to Take π : S 1 × B ε (0 2 ) → B ε (0 2 ) as the projection on the second factor and letπ := π • −1 . Then,π is Lipschitz-continuous and B ε (0 2 ) is contained inπ(S k ∩ U ε (γ )), which entails that We thus conclude that which establishes the inequality (3.1).
Theorem 3.6 (Main existence result) Let a continuous mapping : [0, L] → R 3 satisfying (L) = (0), the global gluing conditions, and the clamping parameters x 0 , t 0 , and d 0 belonging to R 3 be given. If there existsw = (w 1 , x 0 , t 0 , d 0 ) in V such that E KP (w) < +∞ and which complies with the constraints (2.8)-(2.11), then there exists a solution w = (w 1 , x 0 , t 0 , d 0 )  Proof Consider a minimizing sequence {w k } for E KP such that E KP (w k ) < M for some M in R. In particular, we then have that E loop (w k ) < M and we can choose w k in the subset U of V introduced in Lemma 2.4. As in the proof of Theorem 3.1, we can extract a weakly convergent subsequence w k i w with w ∈ U . To complete the proof, it remains to establish that E KP is weakly lower semicontinuous on V . Given the weak lower semicontinuity of E loop , this is tantamount to proving that the functional ( 3.2) is weakly lower semicontinuous. Fix a weakly convergent sequence w k w in U . Let S k belonging to S[w k ] be given by Theorem 3.2 such that Without loss of generality, we can assume that H 2 (S k ) ≤ C for some C > 0. For any k in N, let μ k := H 2 S k . Then, up to the extraction of a subsequence, we have μ k * μ on R 3 and we can set S 0 := spt(μ) \ [w]. On applying Lemma 3.4 with We next show that S 0 belongs to S [w]. Assume by contradiction that there exists γ in D [w] with γ ∩ S 0 = ∅ and pick ε as given by Lemma 3.5. We then find that μ(U ε (γ )) = 0 and, therefore, that which contradicts thesis (3.1) of Lemma 3.5. Hence, we obtain the chain of inequalities which establishes the lower semicontinuity of the functional (3.2).

Concluding Remarks
We introduced a mathematical model for experiments in which a thin filament in the form of a closed loop is spanned by a liquid film. In a variational setting, the model is defined by the sum of the energies of the different components of this system. While the loop is modeled as a nonlinearly elastic rod which is inextensible and unshearable, namely a Kirchhoff rod, the elasticity of the liquid film is described by a homogeneous surface tension. Following Schuricht (2002), we required that the loop satisfies the physical constraints of noninterpenetration of matter, fixed global gluing conditions, and fixed knot type. A crucial point in this treatment is that the cross-sectional thickness of the loop is nonzero, implying that it occupies a nonvanishing volume. This led us to consider a somewhat more physical version of the Plateau problem where the liquid film is represented by a surface with a boundary that is not prescribed, but is free to move along the lateral surface of the three-dimensional bounding loop. Our choice to retain a nonvanishing cross-sectional thickness of the loop, while attributing a vanishing thickness to the liquid film, is justified by the typically large separation between those two length scales.
By exploiting the framework recently proposed by De Lellis et al. (2014) for the Plateau problem, we established the existence of a global minimizer, namely a stable solution, for the coupled Kirchhoff-Plateau problem, in which the boundary of the liquid film lies on the lateral surface of the deformable bounding loop. This was achieved by means of a dimensional reduction, performed in expressing the total energy of the system as a functional of the geometric descriptors of the rod only. To this end, a strongly nonlocal term, entailing the minimization of the surface energy for a fixed shape of the loop, was introduced and the effectiveness of this strategy guaranteed by the auxiliary proof of existence of a surface realizing such a minimization. A key step toward the final result is the adaptation, presented in Lemma 3.4, of the main compactness argument of De Lellis et al. (2014) needed to deal with the deformability of the bounding loop.
Combining the approaches of Schuricht (2002), for the bounding loop, and of De Lellis et al. (2014), concerning the liquid film, we established the existence of a solution that complies with important physical constraints on the topology of the bounding loop and of the surface. Indeed, the latter enjoys the soap film regularity identified by Almgren (1976) and Taylor (1976). An important outcome of our analysis is that the existence of a physically relevant solution is obtained irrespective of the relative strength of surface tension compared to the elastic response of the filament from which the bounding loop is made, which instead influences the multiplicity and the qualitative properties of solutions, as discussed below.
Based on the framework introduced above, further investigations can be addressed in two main directions. First, it would be interesting to study the existence of unstable solutions for the Kirchhoff-Plateau problem. We expect that classical techniques used to establish similar results in the context of the Plateau problem will not be easily applicable here, while suitably adapted tools from nonsmooth critical point theory may prove to be required and effective. It would be also important to investigate the unstable solutions generated by bifurcation phenomena, which are expected to be a feature of our model. Indeed, the competition between the action of the spanning film and the elastic response of the filament can trigger the transition between different stability regimes, meanwhile affecting the multiplicity of solutions to the Kirchhoff-Plateau problem.
Secondly, studying a dissipative dynamics of the system considered above represents a challenging task both from the analytical point of view, again due to the lack of smoothness inherent to our setting, and from the mechanical point of view, since a physically consistent representation of the dissipative phenomena at play might not be straightforward. Nevertheless, a dynamical strategy of the proposed kind would likely provide an excellent basis for the implementation of numerical schemes aimed at finding local minima of the Kirchhoff-Plateau functional.