This paper estimates the relaxed area of a discontinuous map from a disk to three points in the plane, using a Plateau-type problem involving three minimal surfaces coupled at a triple point.
Contribution
It introduces a novel method to bound the relaxed area by relating it to a minimal surface problem with a triple junction, without symmetry assumptions.
Findings
01
Relaxed area is bounded above by a Plateau-type problem solution.
02
Construction of smooth approximations via three coupled minimal surfaces.
03
Dependence of the estimate on the choice of target triple point and connection.
Abstract
In this paper we estimate from above the area of the graph of a singular map u taking a disk to three vectors, the vertices of a triangle, and jumping along three C2− embedded curves that meet transversely at only one point of the disk. We show that the relaxed area can be estimated from above by the solution of a Plateau-type problem involving three entangled nonparametric area-minimizing surfaces. The idea is to "fill the hole" in the graph of the singular map with a sequence of approximating smooth two-codimensional surfaces of graph-type, by imagining three minimal surfaces, placed vertically over the jump of u, coupled together via a triple point in the target triangle. Such a construction depends on the choice of a target triple point, and on a connection passing through it, which dictate the boundary condition for the three minimal surfaces. We show that the…
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Full text
On the relaxed area of the graph of discontinuous maps from the plane to the plane taking three values with no symmetry assumptions
Giovanni Bellettini111
Dipartimento di Ingegneria dell’Informazione e Scienze Matematiche, Università di Siena, 53100 Siena, Italy,
and International Centre for Theoretical Physics ICTP,
Mathematics Section, 34151 Trieste, Italy.
E-mail: [email protected]
Alaa Elshorbagy222
Area of Mathematical Analysis, Modelling, and Applications,
Scuola Internazionale Superiore di Studi Avanzati ”SISSA”,
Via Bonomea, 265 - 34136 Trieste, Italy,
and
International Centre for Theoretical Physics ICTP,
Mathematics Section, 34151 Trieste, Italy E-mail: [email protected]
Maurizio Paolini333
Dipartimento di Matematica e Fisica, Università Cattolica del Sacro Cuore, 25121 Brescia, Italy.
E-mail: [email protected]
Riccardo Scala444 Dipartimento di Matematica “Guido Castelnuovo”, Università La Sapienza, Piazzale Aldo Moro 5, 00185 Roma.
E-mail: [email protected]
Abstract
In this paper we estimate from above the area of the
graph of a singular map u taking a disk to three
vectors, the vertices of a triangle, and
jumping along three C2− embedded curves that meet transversely at only
one point of the disk.
We show that the relaxed area can be estimated from above by the
solution of a Plateau-type problem involving three entangled
nonparametric area-minimizing surfaces.
The idea is to “fill
the hole” in the graph of the singular map
with a sequence of approximating
smooth two-codimensional surfaces of
graph-type, by imagining three minimal surfaces, placed vertically over
the jump of u, coupled together via a triple point in the target
triangle. Such a construction depends on the choice of a target triple
point, and on a connection passing through it, which dictate the
boundary condition for the three minimal surfaces.
We show that the singular part of the relaxed area of u
cannot be larger than what we obtain by minimizing over all possible
target triple points and all corresponding connections.
Let Ω⊂R2 be an open set and v=(v1,v2):Ω→R2
a Lipschitz map. It is well known that the area of
the graph of v is given by
[TABLE]
Extending to nonsmooth maps
via relaxation the definition of the area
is a difficult question
[12],
and is motivated by rather natural problems in calculus of variations:
we can mention for example
the use of direct methods to face
the two-codimensional Plateau problem in R4 in cartesian form,
and the study of lower semicontinuous
envelopes of polyconvex functionals with nonstandard growth [3],
[10].
A crucial issue is to decide
which topology one has to consider in order to compute the relaxed
functional of A(⋅,Ω): of course, the weakest the topology,
the most difficult should be the computation of the relaxed functional, but the easiest
becomes the coerciveness.
We recall
that when v is scalar valued, the natural choice is the L1(Ω)-convergence,
and
the relaxation problem is completely solved [8],
[2];
the L1(Ω)-relaxed functional in this case
consists, besides the absolutely continuous part, of a singular part which
is the total variation of the jump and Cantor
parts of the distributional derivative of v in Ω; in particular,
the relaxed functional, when considered as a function of Ω,
is a measure.
The case of interest here, namely when v takes values
in R2, is much more involved,
due to the nonconvexity of the integrand in (1.1),
and to the unilateral linear growth
[TABLE]
Choosing again the
L1(Ω;R2)-convergence
(as we shall do in this paper), the
relaxed functional A(⋅,Ω) of A(⋅,Ω),
i.e.,
[TABLE]
is, for v∈L1(Ω;R2)∖W1,2(Ω;R2),
far from being understood, and exhibits surprising features.
One of the few known facts that
must be pointed out is that, for a large
class of nonsmooth maps v, the function
Ω→A(v,Ω)
cannot be written as an integral [3], [5], [6];
this interesting phenomenon, related to
nonlocality, has at least two sources.
For simplicity, let us focus our attention on nonsmooth functions
with jumps, thus neglecting
the case of vortices.
The first source of nonlocality
has been enlightened answering to
a conjecture in
[9]. Specifically, consider the symmetric triple junction map usymm,
i.e., the singular map from a disk D of RS2=R2 into RT2=R2,
taking only three values – the vertices of an equilateral triangle Teq⊂RT2 – and
jumping along three segments meeting
at the origin in a triple junction
at equal 120∘ angles:
then
A(usymm,⋅) is not subadditive. This result has been proven in
[3]; subsequently
in [4] it is shown that the value
A(usymm,D)
is related to the solution of three one-codimensional
Plateau-type problems in cartesian form suitably entangled together through the Steiner
point in the triangle Teq. Due to the
special symmetry of the map usymm, the three-problems
collapse together to only one one-codimensional Plateau-type problem in cartesian form,
on a fixed rectangle R whose sides are the radius of D and the side of
Teq.
Positioning three copies of
this minimal surface “vertically” (in the space of graphs, i.e.,
in D×R2) over the jump of usymm
allows, in turn, to construct
a sequence {uε} of Lipschitz maps from D
into R2
the limit area of which improves the upper estimate of [3].
Optimality of this construction
has been shown in
the recent paper
[13], on the basis of a symmetrization procedure for currents.
It is one of the aims of the present paper
to inspect solutions of the above mentioned three Plateau-type problems
in more general situations,
in order to provide upper estimates for
A(u,D), for suitable piecewise constant maps u.
A second source of nonlocality for the functional A(u,Ω) is given by the
interaction of the jump set of a discontinuous map u with
the boundary of the domain Ω.
This phenomenon, already observed in [3] for the map
with one-vortex
at the center of a suitable disk, appears also for functions with jump
discontinuities not piecewise constant [6].
More surprisingly, it appears also for piecewise constant
maps taking three values, provided the jump is sufficiently
close to the boundary of Ω, as observed in
[13], taking as Ω a sufficiently thin
tubular neighbourhood of the jump itself. We shall not be concerned here with this second
source of nonlocality.
As already mentioned above,
in this paper we are interested in estimating from above
the area of the graph of a singular map u taking three
(non collinear) values
and jumping along three embedded curves of class C2 that meet transversely at only one point, see Figure 1.
Let us state this in a more precise way, referring to Sections
2 and 3 for all details.
For simplicity, from now on we fix Ω to be an open disk
D containing the origin 0S in the source
plane R2=Rx,y2=RS2.
Take three non-overlapping non-empty two-dimensional connected regions E1,\leavevmodeE2,\leavevmodeE3
of D
such that
[TABLE]
The three regions are separated by three embedded curves of class C2 (up to the boundary) of length r12,\leavevmoder23,\leavevmoder31 respectively, that meet only at Q (source triple junction);
moreover, each curve is supposed to meet the
boundary of D transversely and we assume also that Q is a
transversal intersection for the
three curves, see Figure 1(a).
Let α1,\leavevmodeα2,\leavevmodeα3 be the vertices of a closed triangle T with non empty interior in the target plane.
Set
[TABLE]
We suppose that
T contains the origin 0T in its interior.
Let us introduce the space X of connections (Definition 3.1 and (3.3), (3.4)); a connection
Γ=(Γ1,Γ2,Γ3) consists of three rectifiable
curves in T,
that connect the vertices of T
to some point
inside the triangle (called target triple point). We shall suppose that each curve
can be written as a graph, possibly with vertical parts, over
the corresponding two sides of T. When Γ
consists of three Lipschitz graphs, we write Γ∈XLip,
and we say that Γ is a Lipschitz connection.
We now show how to construct a new functional
G, consisting of
the sum of the areas of three minimal surfaces – graphs
of three suitable area-minimizing
functions m12, m23, m31 defined on certain
rectangles – coupled together
by the connection considered as a Dirichlet
boundary condition, see Definition 3.4. Set
[TABLE]
Assume Γ∈X. Then
Γij:=Γi∪Γj, ij∈{12,23,31} are (generalized) graphs of functions φij of bounded variation over [0,ℓij].
With a small abuse of notation,
set
[TABLE]
The graph of φ12 on R12 is depicted in Figure 2(a).
Let mij=mij(Γ) be
the unique solution of the Dirichlet-Neumann minimum problem
[TABLE]
where
[TABLE]
Notice that
the minimization is taken among all functions having a Dirichlet condition on three of the four sides of the rectangle Rij;
the missing side corresponds to the intersection points of the jump
with the boundary of D.
From (1.6) it follows that the Dirichlet condition is zero on the sides {0}×[0,rij] and {ℓij}×[0,rij] of Rij;
see Figure 2(b).
Set
[TABLE]
The main result of the present paper reads as follows (see Theorem
4.1 and Corollary 5.8).
Theorem 1.1**.**
Let u:D→{α1,α2,α3} be the discontinuous
BV(D;R2) function defined as
[TABLE]
Then
[TABLE]
This theorem says that the singular part of A(u,D) can be
estimated from above by
[TABLE]
and that such an infimum is a minimum. Intuitively, to
“fill the hole” in the graph of u with smooth two-codimensional approximating
surfaces of graph-type, we start to imagine three minimal surfaces,
placed vertically over the jump of u, coupled together via a
triple point in the target triangle T (notice that the union of these
three minimal surfaces, viewed in D×R2, is not smooth in correspondence
of the source triple junction). Such a construction
depends on the choice of a target triple point, and on
a connection Γ passing through it, dictating
the boundary condition for the three minimal surfaces, over
the sides of the triangle T. Theorem 1.1
asserts that the interesting part of the
relaxed area of u, namely its singular part, cannot be
larger
than what we obtain by minimizing over all possible target triple points
and all corresponding
connections.
As a direct consequence of
the results in [4], [13],
when u=usymm (and 0S is the center of D),
the inequality in (1.10) is an equality, and
the infimum in (1.11) is achieved by the
Steiner graph connecting the three vertices of T (the optimal triple point being the Steiner point, i.e., the barycenter of T).
This seems
to be an interesting result that could be stated purely as
a problem of three entangled area-minimizing surfaces
(each of which lies in a half-space of R4, the three
half-spaces having only {0}×R2 in common)
without referring to the relaxation of the functional
A(⋅,D).
We do not know whether, in general, the
Steiner graph is still the solution of the minimization problem in (1.11),
when no symmetry assumptions (the case we are considering here)
are required.
However it is reasonable to expect that, if in the source
we have symmetry, i.e., the source triple junction is positioned at the center
of D and u
jumps along three segments meeting
at equal 120∘ angles, and if the target triangle T is close
to be equilateral, the inequality in (1.10)
to be still an equality. In this respect, it is worthwhile to
observe that showing a lower estimate,
for instance showing that, in certain cases, the inequality
in (1.10) is an equality, seems difficult.
One of the main technical
obstructions is due to
the poor control on the tangential derivative of vε
in proximity of the jump of a discontinuous L1-limit function v
(see [6]), where {vε}
is a sequence of Lipschitz maps converging
in L1(Ω;R2) to v, and satisfying the uniform bound
supεA(vε,Ω)<+∞.
We also notice that
the symmetrization
methods of [13] cannot be applied anymore, in view
of the lackness of symmetry.
It is worth mentioning that the restriction that we assume
on the connections Γ, namely that each Γi is a graph
(possibly with vertical parts) on the corresponding two sides of T,
cannot be avoided in our approach: indeed, only under this graphicality
assumption we can solve the minimum problem in (1.10) in the class
of surfaces which are graphs over the rectangles Rij. In turn,
the graphicality of such minimal surfaces allows to construct
the sequence {uε}, see (4.11). Removing the graphicality
assumption on Γ requires some change of perspective, and needs
further investigation.
The content of the paper is the following.
In Section 2.1 we recall
some properties of functions of bounded variation of one variable,
the definition of generalized graph (formula (2.1)), and the chain rule.
In Section 2.2 we recall some
properties of Cartesian currents carried by a BV-function.
The functional G, appearing on the right hand side of (1.10), is introduced in Definition 3.4.
In Section 4 we show that
[TABLE]
see Theorem 4.1.
The proof is rather involved, mainly due to technical difficulties:
we first start by supposing that the jump of u is piecewise linear
(Proposition 4.4). Some work is required
to define uε on an ε-strip around the jump
of u and avoiding a neighbourhood of
the source triple junction
(formula (4.11)) and to define uε in
the missing neighbourhood of the source triple junction (step 3
of the proof
of Proposition 4.4):
the construction must be done in such a way that uε remains
Lipschitz, and turns out to be rather involved in the three triangles T1ε,T2ε,T3ε, see Figure 4(b).
In Section 5 we prove that the infimum in (1.11)
is a minimum. The proof is achieved by defining
a topology in the space X which allows
to prove the
density of XLip in X (Lemma 5.2),
the continuity of the functional G (Proposition 5.4)
and the sequential compactness of X (Theorem 5.6).
This latter result is also based on a uniform bound on the length
of the connections (Proposition 5.3), which is
a consequence of the graphicality assumptions on the connections.
2 Some preliminaries
In this section we
recall some results on functions of bounded variation of one variable [2],
and on cartesian currents[12], needed in the sequel.
2.1 Functions of bounded variation in the interval
Let (a,b)⊂R be a bounded open interval and φ∈BV((a,b)); then
•
φ is bounded, and it is
continuous up to an at most countable set of
points of (a,b) denoted by Jφ (jump set);
•
the right and left limits φ(s±) of φ
exist at any s∈(a,b);
the right limit φ(a+)
and the left limit φ(b−) exist.
Thus we may define
[TABLE]
•
the distributional derivative φ′ of φ splits as
[TABLE]
where φ˙ds is the absolutely continuous part and φ˙ is the differential of φ [2, p.138 and Cor.3.33], φ˙(j) and φ˙(c) are the jump and the Cantor part respectively.
We shall always assume that φ is a
good representative in its L1 class such that
φ(s)=φ+(s) for all s∈(a,b);
the pointwise variation of φ
is equal to the total variation ∣φ′∣((a,b)).
The generalized graph of φ is defined as
[TABLE]
and the subgraph of φ as
[TABLE]
We recall that, if φ∈L1((a,b)), then φ∈BV((a,b)) if and only if SGφ,(a,b) has finite perimeter in (a,b)×R.
We denote by ∂−SGφ,(a,b) the reduced boundary of SGφ,(a,b).
We conventionally set φ(a−)=0, φ(b+)=0; in this case we can define Γφ as in (2.1) with (a,b) replaced by [a,b], hence the generalized graph will always pass through the end points of the interval (with possibly vertical
parts over a and b).
The following result can be found for instance in [12, p.486].
Theorem 2.1**.**
Let g∈C1(R) and φ∈BV((a,b)). Then
g∘φ∈BV((a,b)) and
[TABLE]
where n(s,Jφ):=∣φ(s+)−φ(s−)∣φ(s+)−φ(s−) and
δs is the Dirac delta at s.
2.2 Cartesian currents
Let I⊂R be a bounded open interval and φ∈BV(I).
We denote by
[TABLE]
the 2-current in R×R defined as the integration over the subgraph SGφ,I.
The current
[[SGφ,I]]\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0ptI×R
can be also identified with an integer multiplicity current in I×R;
moreover SGφ,I has finite perimeter in I×R so, if ∂[[SGφ,I]]\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0ptI×R denotes the 1-current in I×R defined as its boundary, this results of finite mass.
For future purposes we recall the following result, see [12, Section 4.2.4].
Theorem 2.2**.**
Let φ∈BV(I) and T be the current defined by
[TABLE]
Then T∈D1(I×R)
is a Cartesian current, and
[TABLE]
where ∗ is the Hodge operator and ν(⋅,SGφ,I) is the inward generalized unit normal. Moreover T can be decomposed
into three mutually singular currents
[TABLE]
such that
[TABLE]
where ω=ω1ds+ω2dσ.
The current
T
is boundaryless in I×R, namely ∂T=0.
Furthermore,
if Γφ is the generalized graph of φ
as defined in (2.1), it turns
out that ∂−SGφ,I∩(I×R) is a subset of Γφ and they
differ of a H1-negligible set. Namely
[TABLE]
It easily follows that the current T coincides with the integration over the rectifiable set Γφ (with the correct orientation).
From now on, when the interval is clear from the context, we will simply denote SGφ,I by SGφ.
3 The functional G
In order to prove our main result (Theorem 1.1) we need
some preparation.
Take
three open non-overlapping non-empty connected regions E1,\leavevmodeE2,\leavevmodeE3
of an open disk D, each Ei with non empty interior and with E1∪E2∪E3=D, and let Cij be their boundaries in D
as in the introduction.
Let α1,\leavevmodeα2,\leavevmodeα3 be
the vertices of a closed
triangle T as in Section 1;
we suppose that
T contains the origin 0T in its interior,
and let ℓij be as in (1.4).
Definition 3.1** (Connections in T).**
We say that Γ:=(Γ1,Γ2,Γ3) is a BV graph-type (resp. Lip graph-type) connection in T if Γi,\leavevmodei∈{1,2,3}, are subsets of T such that Γ1∩Γ2=Γ2∩Γ3=Γ3∩Γ1 is one point p of T called target triple point of Γ, αi∈Γi for any i=1,2,3, and
[TABLE]
can be written as the generalized graph (resp. graph) of a function of bounded variation (resp. Lipschitz function) over the closed segment αiαj
(see Figure 1(b)).
Note that the case p∈∂T is not excluded. However, by definition,
if πij:T→Rαiαj, ij∈{12,23,31}, is the orthogonal projection on the line Rαiαj containing αiαj, then πij(p)∈αiαj. Set
[TABLE]
If necessary, in the sequel we will often identify Γij with the (generalized) graph Γφij of a function
[TABLE]
of bounded variation. If T is acute, choosing a suitable cartesian coordinate system where the s-axis is the line Rαiαj, we necessarily have φij(0)=φij(ℓij)=0. In contrast, if the angle of T at αi is greater than or equal to 2π then φij might have a vertical part over αi and φij(0+)>0. In this case the generalized graph of φij does not pass through αi.
In the sequel it will be often convenient to consider an extension of φij on (−∞,0)∪(ℓij,+∞). This extension is denoted by φ~ij. In the case of acute triangle φ~ij is always set equal to [math] on (−∞,0)∪(ℓij,+∞).
Remark 3.2**.**
If for any ij∈{12,23,31}, wij in (3.1) is a point of
continuity of φij then the intersection of the generalized graph of φki and the set [wki,ℓki]×R coincides with Γi which is also the intersection of the generalized graph of φij
with the set [0,wij]×R, where ij,ki∈{12,23,31},\leavevmodeij=ki. If wij is a
discontinuity point of φij this is in general not true, as in Figure 14(b), when i=2.
Remark 3.3**.**
Assume that an angle of T is greater than 2π, say for instance the angle at α1; as already said, the generalized graphs composing a connection Γ are allowed to have vertical parts over α1. The target triple point p of any connection Γ belongs to Tint⊂T, the part of the triangle T which is enclosed between the two lines passing through α1 and orthogonal to α1α2 and α1α3 respectively.
Define the classes:
[TABLE]
Obviously XLip⊂X.
3.1 Useful results on one-codimensional area-minimizing cartesian surfaces
Let
Rij be as in (1.5), and Γ∈X. Then Γij,\leavevmodeij∈{12,23,31} are (generalized) graphs of functions φij of bounded variation over [0,ℓij]. Let B⊂R2 be an open disk containing the doubled rectangle Rij
defined as
[TABLE]
We use for simplicity the same notation φij for the
extension of φij to Rij,
defined as
[TABLE]
and for the extension of φij to a W1,1 function on B∖Rij as in [11, Theorem 2.16].
Let mij=mij(Γ), ij∈{12,23,31}, be the solution of following Dirichlet minimum problem:
[TABLE]
where ∫Rij1+∣Df∣2 is the extension of the area functional to BV(Rij) as defined in [11, Definition 14.1].
From [11, Theorem 15.9] and the fact that the restriction of φij to ∂Rij is continuous up to a countable set of points, it follows that mij solves also
[TABLE]
Let mij=mij(Γ) be the restriction of mij to Rij. Then, by the symmetry of φij with respect to the line {t=rij}, mij is the unique solution of the Dirichlet-Neumann minimum problem
(1.7).
From (1.6) it follows that the Dirichlet condition is zero on the sides {0}×[0,rij] and {ℓij}×[0,rij] of the rectangle Rij. Note that mij is analytic in the interior of Rij but not necessarily Lipschitz in Rij [11, Theorem 14.13], see Figure 2(b).
The properties of the functional G will be discussed in Section 5.
4 Infimum of G as an upper bound of A(u,D)
The aim of this section is to provide
the following upper bound for
A(u,D).
Theorem 4.1**.**
Let u∈BV(D;{α1,α2,α3}) be the function
defined in (1.9). Then
[TABLE]
It is not difficult to see,
by truncating the minimal surfaces in (1.7) with the lateral boundary of the prisms [0,ℓij]×T, that the infimum in (4.1) is the same as the infimum obtained without requiring in Definition 3.1 that Γi⊂T, i∈{1,2,3}.
Lemma 4.2**.**
Let ℓ≥0,\leavevmodep≥0,\leavevmodeφ∈Lip([0,ℓ];[0,+∞)) be such that φ(0)=φ(ℓ)=0 and w,∈[0,ℓ] so that φ(w)=p. Then there exists a sequence {φσ} of C∞ equi-Lipschitz functions in [0,ℓ], converging to φ in L1([0,ℓ]) and uniformly on [0,ℓ] as σ→0+, such that
[TABLE]
Proof.
Let us extend φ in R such that φ(s)=0 in R∖[0,ℓ], so that the extension (still denoted by φ) belongs to Lip(R). Let
φσ(s):=ησ∗φ in R,
where {ησ} is a standard sequence of mollifiers. Hence φσ∈C∞(R), Lip(φσ)≤Lip(φ) and the sequence {φσ} converges uniformly to φ on compact subsets of R. Without
loss of generality we may assume
φσ(s)=0 in R∖(−σ/2,ℓ+σ/2) and φσ(ℓℓ+2σw−σ)=p+cσ,\leavevmodecσ=o(1). Let us first suppose p=0. We define
[TABLE]
It is easy to see that
φσ∈C∞([0,ℓ]),
φσ(0)=φσ(ℓ)=0,\leavevmodeφσ(w)=p,φσ are equi-Lipschitz, and {φσ} converges to φ in L1([0,ℓ]) as σ→0. Notice that the obtained approximation is constantly null in a neighborhood of [math] and ℓ.
In the case p=0, we argue differently. We consider the two intervals [0,w] and [w,ℓ] and we repeat the same approximation above in the single intervals; more precisely we choose two points w1∈(0,w) and w2∈(w,ℓ) with φ(w1)>0, φ(w2)>0 (if these points does not exist it means that the functions are constantly [math] and they are already smooth, so there is nothing to prove). Then we approximate the two functions φ\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0pt(0,w) and φ\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0pt[w,ℓ] as before, and we glue them along w. Note that the glued function is smooth in w since both the two smooth approximations are constantly [math] in a neighborhood of w.
∎
To prove Theorem 4.1 we
use the three area-minimizing functions mij,ij∈{12,23,31}, introduced in Section 3.1, to construct a sequence {uε} of Lipschitz functions that converge to u in L1(D;R2). However mij,ij∈{12,23,31}, are only locally Lipschitz so we need the following smoothing lemma.
Lemma 4.3**.**
Let Γ∈XLip, ij∈{12,23,31}. Let φij=φij(Γij)∈Lip([0,ℓij]), mij=mij(Γij)∈W1,1(Rij), be defined as in Section 3.1. Then there exists a sequence {mijσ} of Lipschitz functions such that mijσ:Rij→R, mijσ=φij\leavevmodeon\leavevmode∂DRij, and
[TABLE]
Proof.
This can be easily proved using an argument similar to the one in [4, p.378: p.381], and using also Lemma 4.2 with the choice w=wij and p=φij(wij). ∎
We start to prove Theorem 4.1 in the special case of a piecewise linear jump, as in Figure 3.
Proposition 4.4**.**
Let u∈BV(D;{α1,α2,α3}) be the map
defined in (1.9) and assume that the jump set of u consists of three distinct segments that meet at the origin and reach the boundary of D. Then (4.1) holds.
Proof.
Let Γ∈XLip be a connection passing through p∈T and G(Γ):=A12(Γ)+A23(Γ)+A31(Γ).
To prove the proposition it is sufficient to construct a sequence {uε}⊂Lip(D;R2) converging to u in L1(D;R2) such that
[TABLE]
Case 1. Assume that the segments separating E1,\leavevmodeE2,\leavevmodeE3 meet at the origin with angles less than π, as in Figure 3.
To simplify the computation we may assume that
p=0T, see Figure 1(b). The idea
of the proof is similar to the one used in [4], with however new difficulties, in particular in Tε (step 3). We will specify various subsets of D and define the sequence {uε} on each of these sets.
Let ε>0 be sufficient small and δε>0 be such that δε→0 as ε→0. Define Tε to be the triangle with the origin 0S in its interior, with vertices ζ1=ζε1, ζ2=ζε2, and ζ3=ζε3, and sides of lengths ε12,\leavevmodeε23,\leavevmodeε31, εij:=∣ζi−ζj∣; the sides of Tε are perpendicular to the lines containing r12, r23, r31 (respectively) and their distance from the origin 0S equals δε. Define three cygar-shaped sets S23ε,\leavevmodeS31ε and S12ε as in Figure 4(a): if for instance y is a coordinate on r12 and x is the perpendicular coordinate, then S12ε is defined as
[TABLE]
where
[TABLE]
Let us set
[TABLE]
Step 1. Definition of uε on E1ε∪E2ε∪E3ε.
We define
[TABLE]
Note that
A(uε,E1ε∪E2ε∪E3ε)=∣E1ε∣+∣E2ε∣+∣E3ε∣, hence
[TABLE]
Step 2. Definition of uε on S23ε∪S31ε∪S12ε .
We will start with the construction on S12ε. Set
[TABLE]
where ⟂ denotes the counterclockwise rotation of π/2.
Let ψ12ε:[δε,r12+cε]→[0,r12]
be linear, increasing, surjective, where cε>0 is the smallest number such that
[TABLE]
Note that for any y∈[δε,r12+cε] we have
[TABLE]
Let m12σ be the map defined in Lemma 4.3, whose area on R12 is by construction close to A12, with {σε}⊂(0,+∞) a sequence such that
[TABLE]
We set, with σ=σε for simplicity,
[TABLE]
Observe that uε=(u1ε,u2ε)∈Lip(S12ε;R2),
uε=α1 on {(x,y)∈S12ε:x=ζ11},
and uε=α2 on {(x,y)∈S12ε:x=ζ12}. By the definition of m12σ, it is uniquely defined the point (depending on ε) wa=(w1a,w2a)∈ζ1ζ2 such that uε(w1a,w2a)=0T (see Figure 4(b)).
Write for simplicity
[TABLE]
Using that ∣ξ∣=∣η∣=1, ξ1η1+ξ2η2=0, and ξ1η2−ξ2η1=1, we compute
[TABLE]
where
ms,mt denote, respectively, the partial derivatives
of m with respect to s:=ε12x−ζ11ℓ12 and t:=ψ12ε(y), and
are evaluated
at (ε12x−ζ11ℓ12\leavevmode,\leavevmodeψ12ε(y)).
As a consequence
[TABLE]
where the last equality follows by the change of variables
Hence, employing the same construction in the strips S23ε and S31ε we obtain
[TABLE]
Step 3. Definition of uε on Tε.
We divide Tε into four closed triangles T1ε,
T2ε, T3ε and T0ε as in Figure 4(b).
We set
[TABLE]
We first define uε on ∂T1ε as follows:
(i)
the value of uε at ζ1
is α1;
(ii)
the value of uε on the side wcwa is
0T.
Note that uε is already defined on the edges ζ1wa and ζ1wc and its graph over both edges is given by a rescaled version of the curve Γ1 suitably parametrized.
More precisely, we recall that
π12:Γ1→α1α2 and π31:Γ1→α3α1 are the orthogonal projections onto the edges α1α2 and α3α1. Since Γ1 is, by hypothesis, a part of a Lipschitz graph, the maps π12\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0ptΓ1 and π31\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0ptΓ1 are bi-Lipschitz bijections between Γ1 and the segments α1π12(p) and α1π31(p), respectively. We know that if (s,t) are coordinates on T with respect to the system with s-axis α1α2, then the inverse of π12\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0ptΓ1 is given by Φ12:α1α2→Γ1,
[TABLE]
Let us denote by L12=L12ε:ζ1wa⊂RS2→α1π12(p)⊂RT2 and L31=L31ε:ζ1wa⊂RS2→α1π31(p)⊂RT2 the linear bijective maps
[TABLE]
Then we define
[TABLE]
and
[TABLE]
compare formula (4.11).
Since Φ12 and Φ31 are Lipschitz with Lipschitz constant independent of ε, and the Lipschitz constants of L12 and L31 have order ε1, it follows that the Lipschitz constants of uε over the segments ζ1wa and ζ1wc have order ε1.
Now we want to define uε in the interior of T1ε. First we observe that the map π31∘Φ12:α1π12(p)⊂RT2→α1π31(p)⊂RT2 is a bi-Lipschitz bijection, with constant independent of ε. A direct computation then provides that the map Ψ:ζ1wa⊂RS2→ζ1wc⊂RS2 defined by
[TABLE]
is bi-Lipschitz between ζ1wa and ζ1wc with bi-Lipschitz constants of order 1 as ε→0+. Given Q∈ζ1wa let Q′:=Ψ(Q)∈ζ1wc.
Now we show that T1ε is
foliated by the segments QQ′, i.e., for any R∈T1ε we can find a unique Q∈ζ1wa for which R∈QQ′.
First we notice that
QQ′∩SS′=∅ for any Q=S∈ζ1wa with Q′=Ψ(Q) and S′=Ψ(S).
Indeed, thanks to the fact that Ψ is a homeomorphism and that it keeps ζ1 fixed, it is
easy to see that if S∈ζ1Q then S′∈ζ1Q′, or if Q∈ζ1S then Q′∈ζ1S′. Consider the function
[TABLE]
where τ:=∣wa−ζ1∣wa−ζ1 and ν(q):=∣Ψ(qτ)−qτ∣Ψ(qτ)−qτ. It is clear that the image of f is a closed set and Im(f)={QQ′:Q∈ζ1wa,Q′=Ψ(Q)}. Now we show that
Im(f)=T1ε.
Assume by contradiction there is R∈T1ε∖Im(f) and take a disk B⊂T1ε∖Im(f) centered at R. Let Qr,Ql∈ζ1wa be such that qr:=∣Qr−ζ1∣ (resp. ql=∣Ql−ζ1∣) be the supremum (resp. the infimum) parameter for which B lies on the right (resp. left) of QrQr′ (resp. QlQl′). Note that Qr=Ql due to the injectivity of Ψ, thus for any Q∈QrQl the segment QQ′ must intersect B, a contradiction, see
Figure 6(a).
Hence we may define uε on T1ε as
[TABLE]
We want now to show that on T1ε, uε is Lipschitz continuous with Lipschitz constant of order ε1. To prove this let us fix R∈T1ε. By definition uε(R)=uε(Q) for some Q∈ζ1wa and uε is constant on the segment QQ′∋R.
Let e:T1ε→ζ1wa be the function taking (x,y)∈T1ε to the intersection point of ζ1wa and the line passing through (x,y) parallel to QQ′.
Let g:T1ε→ζ1wc be the function taking (x,y)∈T1ε to the intersection point of ζ1wc and the line passing through (x,y) parallel to QQ′. Let R^∈T1ε a point in T1ε, we want to estimate the ratio
[TABLE]
Consider the two segments Qe(R^) and Qg(R^). By definition R^∈SS′ and uε(R^)=uε(S)=uε(S′) for two points S∈ζ1wa and S′∈ζ1wc. It is straightforward that either S∈Qe(R^) or S′∈Q′g(R^). Without loss of generality suppose the first case holds, see Figure 6(b).
Finally, denote by θ the angle between QQ′ and ζ1wa and by θ′ the angle between QQ′ and ζ1wc. Using the fact that the homeomorphism in (4.19) is bi-Lipschitz with constant of order 1 it is not difficult to see that there is a constant θ0>0 independent of ε such that min{θ,θ′}≥θ0.
This is a consequence of the fact that the bi-Lipschitz constant of Ψ in (4.19) is of order 1. Indeed, if L=lip(Ψ) and 1/L′=lip(Ψ−1), we see that
[TABLE]
hence
[TABLE]
where θζ1 is the angle at ζ1 (here we have used the law of sines and that θ′=π−θζ1−θ). A similar estimate holding for θ′, this readily provides the boundedness from below of min{θ,θ′}.
As a consequence we have
[TABLE]
Thus, we compute
[TABLE]
that is bounded by the Lipschitz constant of Φ12∘L12 which is of order ε1.
Eventually we compute the Jacobian of uε in (4.20). By construction the image of T1ε by uε is exactly the curve Γ1, which has zero Lebesgue measure in R2. By a standard application of the area formula it follows that the Jacobian of uε is vanishes a.e. in T1ε. We have concluded the definition of uε in T1ε.
The constructions on T2ε and on T3ε are similar, and similar estimates of the derivatives and Jacobian hold.
Using that the area of the triangle Tε is of order ε2, we have
[TABLE]
From (4.7), (4.11),
(4.16), (4.20), and the estimates above
it follows that
Case 2. Assume that two of the segments separating E1,\leavevmodeE2,\leavevmodeE3 meet at the origin with an angle greater than or equal to π.
Similar to Case 1, we divide the domain D into a finite number of subsets and define the sequence {uε} on each of these sets. Draw the normal to each segment at the point of distance δ from the origin. The normal lines meet at two points ζ1, ζ2. Divide D
into three cygar-shape subsets S23ε,\leavevmodeS31ε,\leavevmodeS12ε (with widths of order δ=O(ε)) and a quadrilateral Hε as in Figure 7(a).
Let
[TABLE]
Set
[TABLE]
Define uε on S23ε∪S31ε∪S12ε as in
Step 2 in case 1. It remains to define uε on Hε.
Recall that by construction there exist uniquely determined three points wa∈ζ1ζ2, wb∈ζ2ζ3 and wc∈ζ1ζ4 such that
[TABLE]
Divide Hε into six triangles T0ε,\leavevmodeT1ε,\leavevmodeT2ε,\leavevmodeT3ε,\leavevmodeT4ε,\leavevmodeT5ε, as in Figure 7(b), where wd is any point in ζ3ζ4 and wd=ζ3,\leavevmodewd=ζ4.
Set
[TABLE]
We define uε in the triangles T1ε and T2ε as in
Step 3; it remains to define uε on T3ε,\leavevmodeT4ε,\leavevmodeT5ε.
Let us first define uε on the edges wcwd and wbwd. The map uε is already defined on the other edges, and its graph over ζ4wc and ζ3wb is given by a suitable reparametrization of the curve Γ3, whereas uε on ζ4ζ3 is constantly α3. Therefore it suffices to define uε in such a way its graph over wcwd and wbwd coincides with Γ3 as well, and then we can define uε inside T3ε and T4ε using the same construction for T1ε in
step 3. Similarly, using that the graph of uε on wcwd and wbwd is again Γ3, we can repeat the construction in the triangle T5ε.
Following the computation as in case 1 we get (4.4).
This concludes the proof.
We will suitably adapt the construction made in the proof of Proposition 4.4. By hypothesis the regions E1,E2,E3 are enclosed by C2-embedded curves Cij,ij∈{12,23,31}, parametrized by arc length cij:[0,rij]→R2,\leavevmodeij∈{23,31,12}. Moreover such curves meet ∂D transversely and intersect each other (transversely) only at one point Q. Suppose that the angles formed at Q by the three curves are all less than π (the other case is similarly adapted from the corresponding case in the proof of Proposition 4.4).
We will divide the domain D into a finite number of subsets and define the sequence {uε} on each of these sets.
Let δε>0 be such that δε→0 as ε→0+. Let τ∈[0,rij] be an arc lenght parameter on Cij, with orthogonal coordinate d that coincides with the signed distance from Cij negative in Ei and positive in Ej. Let Qij∈Cij be the point with arc distance τ=δε from the origin Q. Consider the three lines normal to Cij at Qij. For δε sufficiently small, since the angles at the origin are less than π and the curves are of class C2 up to the closure, these lines mutually meet at points ζ1, ζ2, and ζ3. Let εij be the length of ζiζj, which are of order ε. The tubular coordinates of the points ζ1 and ζ2 with respect to C12 are (d1,δε,) and (d2,δε), with d2−d1=ε12,d1<0,d2>0. For δε small enough we can consider the cylindrical neighborhood of C12 defined as
[TABLE]
where we have prolonged C12 outside D for convenience.
Similarly we define S23ε and S31ε. Let Tε be the triangle with vertices ζ1, ζ2, and ζ3.
Finally, let E1ε,\leavevmodeE2ε,\leavevmodeE3ε be defined as in (4.6), and uε as in (4.7).
Step 1. Definition of uε on S12ε∪S23ε∪S31ε.
We do the construction on S12ε, and uε will be defined similarly on S23ε and S31ε. We know that c12([δε,r12])=C12∩S12ε. The system of coordinates (d,τ) defines a C1-diffeomorphism h between the rectangle [d1,d2]×[δε,ρ12] and its image N12ε,δ which contains S12ε, namely
[TABLE]
where νˉ(τ) is the unit normal vector pointing toward E2 at c12(τ) and ρ12=ρ12ε≥r12 is the infimum of those ρ for which S12ε⊂N12ε,δ, see Figure 8(b).
Since h is a C1-diffeomorphism we have that
[TABLE]
is the inverse of h and is of class C1. We want to estimate the Jacobian of h−1. To this aim, we first see that ∇d(c12(τ))=νˉ(τ) since c12([δε,ρ12]) is the zero level set of d and, from [1, Rem. 3(1)], we have
[TABLE]
Fix τ∈[δε,ρ12]; by definition of tubular coordinates the segment {c12(τ)+dνˉ(τ):d∈[d1,d2]} is a level set of the function τ(⋅), hence,
[TABLE]
therefore
[TABLE]
Thus the Jacobian of h−1 will be
[TABLE]
since ∣∇d∣=1 in N12ε,δ.
Let us compute ∇τ; fix d∈(d1,d2) and define c12d(τ):=c12(τ)+dνˉ(τ).
Now recall (4.28) and that
(c12d)′(τ) is parallel to νˉ⊥(τ), so that
[TABLE]
Let us recall that C12 is parametrized by arc length, i.e., ∣c12′(τ)∣=1, so that νˉ′(τ)=∣c12′′(τ)∣c12′(τ). Thus (c12d)′(τ)=(1+d∣c12′′(τ)∣)c12′(τ). Since τ∘c12d=Id it follows that ∇τ(c12d(τ))T(c12d)′(τ)=∇τ(c12d(τ))⋅(c12d)′(τ)=1. Therefore, from (4.31), we deduce
[TABLE]
and in particular limd→0∣∇τ∣=1 uniformly in S12ε.
We are ready to define uε in S12ε. We first set ψ12ε as in (4.9) with r12+cε=ρ12, i.e., ψ12ε(τ)=κε(z−δε), setting κε:=ρ12−δερ12. Then we define u~ε on [d1,d2]×[δε,ρ12] as in the right hand side of (4.11) and set
[TABLE]
Explicitly, recalling that ξ=ℓ12α2−α1 and η=ξ⊥, for (x,y)∈S12ε we have
[TABLE]
Observe that uε=(u1ε,u2ε)∈Lip(S12ε;R2),
uε=α1 on {(x,y)∈S12ε:d(x,y)=d1}, uε=α2 on {(x,y)∈S12ε:d(x,y)=d2}, and by construction there exists wa∈h([d1,d2]×{δε}) such that uε(wa)=0T.
Write for simplicity m=m12σ.
We have
[TABLE]
where
ms,mt denote the partial derivatives
of m with respect to s=ε12d(x,y)−d1ℓ12 and t=κε(τ(x,y)−δε) respectively, and
are evaluated
at (ε12d(x,y)−d1ℓ12\leavevmode,\leavevmodeκε(τ(x,y)−δε)).
Hence
[TABLE]
where
we have used ∣ξ∣=∣η∣=1 and ξ1η1+ξ2η2=0. From (4.29) we have
[TABLE]
Moreover
[TABLE]
where again ms,mt are evaluated
at (ε12d(x,y)−d1ℓ12\leavevmode,\leavevmodeκε(τ(x,y)−δε)), and
we have used (4.29), (4.30), and
ξ1η2−ξ2η1=1.
Therefore from (4.35) and (4.36) we obtain
[TABLE]
As a consequence
[TABLE]
where ms,mt in the first integral are evaluated at (ε12d(x,y)−d1ℓ12\leavevmode,\leavevmodeκε(τ(x,y)−δε)), ∇τ in the second integral is evaluated at (x,y)=Φ−1(s,t) and the last equality follows from the change of variables
[TABLE]
and Pε:=R12∖Φ(S12ε)
(see Figure 5). Here one checks that Φ=H∘h−1 with H(d,τ)=(ε12d−d1ℓ12\leavevmode,\leavevmodeκε(τ−δε)) so that, using (4.30), the Jacobian of the change of variable is ∣∇τ(Φ−1(s,t))∣1ℓ12κεε12.
Hence, recalling (4.32) and that κε→1 as ε→0+,
[TABLE]
Now, let us recall that m~=m12σ is the approximating function as in (4.3); it follows that
[TABLE]
Hence, employing the same construction in the strips S23ε and S31ε,
and using (4.39) we obtain from a diagonal argument with σ=σε→0 as ε→0+,
[TABLE]
Step 2. Definition of uε on Tε. This is identical to
Step 3 of the proof of Proposition 4.4 and therefore {uε}⊂Lip(Br;R2) and (4.23) holds. Following the same computations of Proposition 4.4 the conclusion follows.
Step 3.
For the case where two of the curves Cij,ij∈{12,23,31} meet at Q with an angle larger than or equal to π we replace Tε with Hε defined in case 2 of Proposition 4.4, in the above construction.
∎
5 Existence of minimizers for the functional G
Let D be an open disk centered at the origin such that E1,\leavevmodeE2,\leavevmodeE3 are circular sectors with 120∘ angles and let T be an equilateral triangle. Let p be the barycenter of T and Γi be the segment connecting αi and p, i∈{1,2,3}. Hence Γ=(Γ1,Γ2,Γ3)∈XLip so that
[TABLE]
Moreover we have
[TABLE]
where u=usymm (see Section 1), and the equality follows from [13, Section 3]
and the inequality follows from Proposition 4.4. Thus
[TABLE]
Hence in this symmetric situation the optimal connection is obtained through the Steiner graph connecting α1,α2 and α3. This motivates the analysis of this section, which is carried on without symmetry assumptions.
We recall that given a connection Γ=(Γ1,Γ2,Γ3)∈X we denote by φij=φij(Γij):[0,ℓij]→R the function whose graph is Γij=Γi∪Γj (see (3.2)).
Definition 5.1** (Convergence in X).**
We say that a sequence {Γn}⊂X converges to Γ∈X in X, and we write Γn→Γ in X, if
[TABLE]
5.1 Density and approximation
We start to show that a BV connection Γ∈X can be approximated by Lipschitz connections; the difficulty is to keep graphicality of each branch of of the approximating connections
with respect to the two corresponding edges of T at the same time.
Recall that Γi is the branch of the connection Γ connecting αi to p and that by Definition 3.1 we have
[TABLE]
Note that we excluded the vertical parts over the points πij(p),\leavevmodeij∈{12,23,31}, due to Remark 3.2; however we still have
[TABLE]
Lemma 5.2** (Piecewise linear approximation).**
For any Γ∈X with target triple point p∈T there exists a sequence {Γn}⊂XLip of connections with target triple point p such that φij(Γijn),ij∈{12,23,31}, is a
piecewise linear555This means that it is Lipschitz piecewise
linear with at most finitely many points of nondifferentiability.
function,
[TABLE]
and
[TABLE]
Proof.
Let ij=12 and let w12 be defined as in (3.1).
Let n12:=(0,1)∈R2 be the inward unit normal to α1α2, n31:=(α,β) be the inward unit normal to α3α1, and ν(sˉ):=(ν1(sˉ),ν2(sˉ)) be the generalized outward unit normal at the point (sˉ,φ12(sˉ)) to the generalized graph Γφ12\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0pt[0,w12] of φ12\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0pt[0,w12] (for all sˉ where it exists), see Figure 9. Without
loss of generality we may assume Γ1=Γφ12\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0pt[0,w12].
We start to show that φ12 cannot have too negative slope, otherwise Γ1 loses graphicality with respect to α3α1.
Step 1.
We claim that
[TABLE]
in the sense of measures, i.e.,
[TABLE]
From the graphicality with respect to α1w31 we have, for all sˉ where ν(sˉ) exists,
[TABLE]
Set
[TABLE]
note that ν(sˉ)=1+(φ˙12(sˉ))21(−φ˙12(sˉ),1) for any sˉ∈Ir. From [12, Thm. 7 p. 301 and Thm. 5 p. 379], we have
where μ:=(φ12′,−L1)=(−ν1,−ν2)∣μ∣
and the second equality follows from the first formula in (5.6).
For any Borel set B⊆[0,w12] we deduce
[TABLE]
Step 2.
Given ϵ∈(0,1), we choose n=n(ϵ)∈N and points
[TABLE]
such that each ξi,\leavevmodei∈{1,⋯,n−1}, is a point of
continuity of φ12, and if we define φn∈Lip([0,w12]) as the piecewise linear interpolation with
[TABLE]
then
[TABLE]
The graph of φn may still have vertical parts over w31α1.
Indeed from [2, Theorem 3.30], and the fact that ξi are
continuity points of φ12 we have
[TABLE]
and equality may hold, hence the graph of φn over π31(p)α1 may have finitely many vertical parts.
It is now sufficient to repeat the argument with φn in place of φ12, choosing a suitable partition of [w31,ℓ31], so to ensure that (out of finitely many points)
[TABLE]
In this way φn is a Lipschitz graph also with respect to π31(p)α1.
Step 3.
We have
[TABLE]
where Φ12∈BV([0,w12];R2) is defined as Φ12(ξ):=(ξ,φ12(ξ)), and the last equality follows from [2, (3.24), p.136].
Step 4.
Define
[TABLE]
Similarly we define Γ2n and Γ3n, and we set Γijn:=Γin∪Γjn.
Then Γn:=(Γ1n,Γ2n,Γ3n),
satisfies the required properties.
∎
Proposition 5.3** (Uniform estimate of the length).**
There exists c>0 such that for all Γ∈X we have
[TABLE]
Proof.
Let Γ∈X be a connection through p∈T. Without
loss of generality we may assume that p=α1. From (5.8) we have
[TABLE]
Choose a partition
[TABLE]
Let Γ1h be the piecewise linear interpolation connecting (ξi−1,φ12(ξi−1)) and (ξi,φ12(ξi)),i∈{1,⋯,h}. The unit tangent to Γ1h is enclosed in the angle formed by n12 and n31, the unit normals to α1α2 and α3α1 (due to the graphicality condition with respect to α1α2 and α3α1), see Figure 10. It follows that Γ1h is the graph of a function ϕ12h over the segment α1p. Fix a Cartesian coordinate system in which the
t-axis is the line α1p and the origin is α1.
For any t∈[0,∣α1−p∣]
(up to a finite set) let τ(t) be the unit tangent to Γ1h at (t,ϕ12h(t)) and let n=(1,0) and n⊥=(0,1).
Hence ϕ12h′=τ⋅nτ⋅n⊥ satisfies
[TABLE]
Note that one between ∣c1−∣ and ∣c1+∣ might be +∞, since one of the sides α1α2 or α3α1 can be horizontal (this happens only if the point p is on one side of the triangle). However
we always have that c1−≤0, c1+≥0.
Furthermore, when the angle α^1 in α1 is less or equal to 2π, it follows that c~1:=min{∣c1−∣,∣c1+∣}≤∣tan(2π−2α^1)∣. In the case that α^1>2π, thanks to the fact that p∈Tint, we have max{∣c1−∣,∣c1+∣}≤∣tan(π−α^1)∣. Thus the only difficulty to prove that the length of Γ1h is controlled when α^1≤2π. So let us assume this and in addition that ∣c1−∣=c~1 (the other case is similar).
Since ϕ12h(∣α1−p∣)=ϕ12h(0)=0 we have
[TABLE]
where (ϕ12h′)+ and (ϕ12h′)− are the positive and negative parts of the measure ϕ12h′=ϕ˙12hdt,
thus we estimate
[TABLE]
Defining c1 as the right-hand side of the last inequality we see that c1 is a positive constant depending only on the geometry of T.
Similarly we may show that H1(Γ2)≤c2 and H1(Γ3)≤c3 for c2,c3>0 depending only on T. This proves (5.9) with c=c1+c2+c3.
∎
The next lemma shows
continuity of the sum of the three areas of area minimizing surfaces defining G in (3.9), with respect to the L1 convergence of the traces in T.
Proposition 5.4** (Continuity of G).**
Let Γ∈X, and let {Γn}⊂X be a sequence converging to Γ in X.
Then
[TABLE]
Proof.
Since Γ∈X and {Γn}⊂X we have φij∈\leavevmodeBV([0,ℓij]) and {φijn}⊂\leavevmodeBV([0,ℓij]) where φij:=φij(Γij), φijn:=φij(Γijn).
Hence from (1.6) and Section 3.1 it follows that there exist mij,\leavevmodemijn∈W1,1(Rij) such that
[TABLE]
where we recall that Rij is the double rectangle defined in (3.5) and φij,\leavevmodeφijn are extended on a disk B containing Rij as in Section 3.1.
Define mijn and mij as
[TABLE]
so that mijn,\leavevmodemij∈BV(B).
Since mijn is competitor in (5.15) and mij is competitor in (5.14)
we have, recalling also the discussion leading to (3.8),
[TABLE]
Thus
[TABLE]
Recall that mijn (resp. mij) is the restriction of mijn (resp. mij) to Rij. Hence, from (3.9), (5.1) and (5.16),
(5.13) follows.
∎
Corollary 5.5**.**
We have
[TABLE]
5.2 Compactness of the class X
The aim of this section is to show that the infimum in (5.17) is attained. To do this we need the following result.
Theorem 5.6** (Compactness).**
Any sequence {Γn}⊂X admits a subsequence converging in X to some Γ∈X.
Remark 5.7**.**
In Definition 5.1
it is required convergence of {Γn} to Γ in L1. For this
reason, if Γn has target triple point pn, it is not guaranteed
that the point b:=limn→+∞pn (it exists up to subsequences)
still belongs to Γij for all ij, see Figures 12
and 15(a). As a consequence, if {Γn} converges to Γ it is not true, in general, that pn→p, where p is the target triple point of Γ.
Proof.
Let {Γn}⊂X and φijn=φij(Γijn),\leavevmodeij∈{12,23,31}.
From Proposition 5.3{φijn} is
uniformly bounded in BV([0,ℓij]) for any ij∈{12,23,31}. Thus, up to a not relabelled subsequence, there exists φij∈BV([0,ℓij]) such that
[TABLE]
We shall adopt our usual convention
[TABLE]
Denote by \tensor[]Γij⊂R2 the limit graph over (the closed segment) αiαj that we identify with the generalized graph of φij over [0,ℓij].
Since T is closed and convex we have
Γij⊂T; moreover, by construction, αi and αj are the endpoints of Γij.
Notice that if we assume that T is acute,
this excludes the presence of
vertical parts over its vertices.
It remains to prove
that the three obtained curves Γij, ij∈{12,23,31}, form
a BV connection; in particular that they intersect mutually in a unique
well-defined point.
We claim that
there exists a unique p∈⋂ij\tensor[]Γij that divides each \tensor[]Γij into two curves Γijl and Γijr
such that
[TABLE]
Let us denote by φ~ijn
the extension to R of the function φijn
vanishing
in (−∞,0)∪(ℓij,+∞).
Similarly φ~ij is the extension of φij
vanishing
in (−∞,0)∪(ℓij,+∞).
Consider the sequence {[[SGφ~ijn]]}n⊂D2(R2) of 2-currents regarded in R2 and the 2-current [[SGφ~ij]]. Their boundaries are the currents carried by the graphs of φ~ijn and φ~ij, as defined in Theorem 2.2.
The 1-currents carried by the graph of φijn and φij, by convention (5.20), coincide with the restrictions of ∂[[SGφ~ijn]] and ∂[[SGφ~ij]] to the closed set [0,ℓij]×R.
Namely, if we denote by
[TABLE]
then
[TABLE]
where Lij is the 1-current given by integration over the two
halflines (−∞,0)×{0}∪(ℓij,+∞)×{0}.
The curves Γijn and Γij coincide with the support of [[Γijn]] and [[\tensor[]Γij]], respectively.
We now prove our claim in three steps.
Step 1. The currents [[Γijn]] converge
(up to a not relabelled subsequence) weakly in the sense of currents to [[\tensor[]Γij]], i.e.,
Finally (5.23) follows from Lemma 5.3 and the
weak lower semicontinuity of the mass of currents, and the proof of
step 1 is concluded.
It is not restrictive to assume that
wijn=∣αi−πij(pn)∣ is a point of continuity of φijn
for all n∈N and all ij∈{12,23,31}.
Indeed given a sequence {Γn}⊂X converging to Γ,
from Lemma 5.2
for all n
we can
assign a sequence {Γm,n}⊂XLip such that
Γm,n→Γn as m→+∞. Thus by a diagonal argument, we find a sequence {Γm(n),n}⊂XLip which tends to Γ and satisfies
the above requirement
(we can also assume that
Γn is Lipschitz, but this will not
be needed in the proof).
Without loss of generality
(up to a not relabeled subsequence) we may
further assume
[TABLE]
{wijn} is a monotone sequence, and
[TABLE]
Before passing to the second
step, it is convenient to divide the target triangle T into
various regions.
Assume first that T is acute.
The point b, together with the heights
[TABLE]
divides T
into three regions Pi,\leavevmodei∈{1,2,3},
as shown in Figure 11(a); precisely,
if Pi denotes the closed region enclosed by hij, hki, αiπij(b)
and αiπki(b), then Pi is defined by
[TABLE]
Similarly we define hijn and Pin by replacing b with pn in (5.24) and (5.25).
Assume now that T is not acute. Without
loss of generality we may assume that the angle at α1 is greater than 2π. The only difference here is with the definition of P1n and P1 since each \tensor[]Γij has to satisfies the graphicality condition with respect to αiαj; hence we define P1 as the closed quadrilateral bounded by h12, h31, m12 and m31, where m12 and m31 are the normals to α1α2 and α3α1, respectively, passing through α1 (see Figure 11(b)).
Similarly we define P1n. Finally we set P1:=P1∖(h12∪h31) and P1n:=P1n∖(h12∪h31).
Step 2. We will prove that we can decompose Γ12∪Γ23∪Γ31 as three currents meeting at a point b.
It is
easy to see that the sets Pin are converging to Pi with respect to the Hausdorff distance. It is not true in general that Γijn is converging to Γij with respect to the Hausdorff distance (see Figure 12); however, since
[TABLE]
for all ij,ki∈{12,23,31},ij=ki, it is readily seen that
[TABLE]
For i∈{1,2,3}, the
integral 1-current [[Γin]]=[[Γijn]]\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0pt[0,wij]×R has boundary δpn−δαi in R2. By the
compactness theorem for integral currents [12, Theorem 2, p.141] there exists an integral current \tensor[]Ti∈D1(R2), i=1,2,3, such that, up to a not relabeled subsequence,
[TABLE]
Clearly
[TABLE]
From (5.26) and thanks to the convergence of Pin to Pi with respect to the Hausdorff distance, we infer
[TABLE]
where ij,ki∈{12,23,31}.
Note that \tensor[]Ti is not necessarily equal to [[\tensor[]Γij]]\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0ptPi, due to a possible cancellation of a vertical part over πij(b),\leavevmodeij∈{12,23,31} (that is, on hij), see Figure 12.
However from [[Γijn]]=[[Γin]]+[[Γjn]] and (5.28) we have
[TABLE]
as currents in R2. Notice that \tensor[]Ti and \tensor[]Tj have multiplicity one, and in (5.31) they contribute with opposite orientation. This allows, if necessary, to identify \tensor[]Ti,\leavevmodei=1,2,3, with its support.
Note also that \tensor[]Ti may have vertical part over αi, see Figure 13.
Now, since [[Γin]] is Cartesian with respect to both the edges αiαj and αkαi, from (5.28) it follows that \tensor[]Ti\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0ptPi is part of two generalized graphs over the same edges, i.e.,
[TABLE]
Moreover, we infer that Ti cannot have vertical part over hij and hki at the same time; in other words once the current \tensor[]Ti touches one of the heights hij or hki it stays there until it reaches b, and \tensor[]Ti cannot have a nonempty support in more than one height, see Figures 14(b)-15(b).
We conclude the following statement:
(A)
The supports of the three currents Ti, i=1,2,3, have as common point b. Moreover, if there are i=j such that the supports of Ti and Tj intersect in a point different from b, then this intersection occurs on the mutual height hij. Finally, if the supports of Ti and Tj intersect on hij outside b, then they intersect on a closed segment and the intersection of the supports of Ti and Tj with Tk is only the point b.
Step 3.
To conclude the proof of our claim
we now analyse the possible cases arising from (A).
Case (i). Assume that the three supports of the currents Ti, i=1,2,3, intersect only at the point b. This includes the case
[TABLE]
as in Figure 14(a). But it may also happen that Ti has vertical part over hij, provided that Tj does not have vertical part over the same height (see for instance Figure 14(b)).
In any case we may set
[TABLE]
where we have identified the currents Ti with their supports.
By (5.29) and (5.32), the claim is achieved.
Case (ii). The second case which must be discussed is the one considering possible overlapping of the support of the currents Ti. By condition (A)
such overlapping, giving rise to cancellations, can occur only on one height hij. Hence, assume there exists one (and only one) ij∈{12,23,31} such that
[TABLE]
Thus we have \tensor[]Ti\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0pthki=0 and \tensor[]Tj\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0pthjk=0.
First assume that [[\tensor[]Γij]]\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0pthij=0, i.e., φij is continuous at wij. Then \tensor[]Ti\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0pthij=\tensor[]Tj\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0pthij, see Figure 15(a). We set, identifying Ti with its support,
[TABLE]
One checks that the connection built above is a BV graph type connection, addressing the claim.
Now assume that
[TABLE]
i.e., φij jumps at wij. Thus either \leavevmodespt\leavevmode[[\tensor[]Γij]]\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0pthij⊆\leavevmodespt\leavevmode\tensor[]Ti\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0pthij or \leavevmodespt\leavevmode[[\tensor[]Γij]]\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0pthij⊆\leavevmodespt\leavevmode\tensor[]Tj\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0pthij. Without
loss of generality we may assume that \leavevmodespt\leavevmode[[\tensor[]Γij]]\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0pthij⊆\leavevmodespt\leavevmode\tensor[]Ti\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0pthij hence
\tensor[]Ti\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0pt(hij∖\leavevmodespt\leavevmode[[\tensor[]Γij]])=−\tensor[]Tj\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0pthij=0 (note that \leavevmodespt\leavevmode[[\tensor[]Γij]]\vruleheight=7.0pt,width=0.5pt,depth=0.0pt\vruleheight=0.5pt,width=6.0pt,depth=0.0pthij={tφij(wij+)+(1−t)φij(wij−):t∈[0,1]}).
We set
[TABLE]
see Figure 15(b).
Also in this case the conclusion follows.
In the end we define
[TABLE]
The proof is achieved.
∎
From compactness of the space of BV connection, combining with Proposition 5.4, we see that the infimum in (4.1) is attained. As a consequence, we can conclude the proof of Theorem 1.1.
Corollary 5.8**.**
We have
[TABLE]
Acknowledgements
The present paper benefits from the support of the GNAMPA (Gruppo Nazionale per l’Analisi Matematica,
la Probabilità e le loro Applicazioni) of INdAM
(Istituto Nazionale di Alta Matematica).
Bibliography13
The reference list from the paper itself. Each links out to its DOI / PubMed record.
1[1] L. Ambrosio, Geometric evolution problems, distance function and viscosity solutions. Calculus of Variations and Partial Differential Equations (Pisa,1996), 5-93, Springer, Berlin, (2000).
2[2] L. Ambrosio, N. Fusco, D. Pallara, Functions of Bounded Variation and Free Discontinuity Problems, Mathematical Monographs, Oxford University Press, 2000.
3[3] E. Acerbi, G. Dal Maso, New lower semicontinuity results for polyconvex integrals , Calc. Var. Partial Differential Equations 2 (1994), 329–371.
4[4] G. Bellettini, M. Paolini, On the area of the graph of a singluar map from the plane to the plane taking three values , Adv. Calc. Var. 3 (2010), 371–386.
5[5] G. Bellettini, M. Paolini, L. Tealdi, On the area of the graph of a piecewise smooth map from the plane to the plane with a curve discontinuity , ESAIM: Control, Optimization and Calculus of Variations 22 (2015), 29–63.
6[6] G. Bellettini, M. Paolini, L. Tealdi, Semicartesian surfaces and the relaxed area of maps from the plane to the plane with a line discontinuity , Ann. Mat. Pura Appl. 195 (2016), 2131-2170.
7[7] H. F. Cullen, Introduction to General Topology, MA:Health, Boston, 1968.
8[8] G. Dal Maso, Integral representation on B V ( Ω ) 𝐵 𝑉 Ω BV(\Omega) of Γ Γ \Gamma -limits of variational integrals , Manuscripta Mat. 30 (1980), 387-416.