This paper investigates the steady creep flow in fiber-reinforced viscoplastic solids with high viscosity contrast, revealing new anisotropic effects influenced by the Norton exponent.
Contribution
It introduces a novel analysis of how high viscosity contrast affects the homogenization and anisotropic behavior in Norton-Hoff fibered composites.
Findings
01
Identification of new anisotropic effects due to high contrast
02
Influence of Norton exponent on flow behavior
03
Enhanced understanding of composite creep flow
Abstract
We study the steady creep flow of a perfectly viscoplastic solid reinforced by fibers with high viscosity contrast. Our study unveils new effects related to anisotropy and conditioned by the Norton exponent.
Equations1069
(Pε):⎩⎨⎧uε∈Wb1,p(Ω;R3)infFε(uε)−∫Ωf.uεdx,Fε(uε):=∫Ω∖Trεf(e(uε))dx+lε∫Trεg(e(uε))dx,e(uε)=21(∇uε+∇Tuε),f∈Lp′(Ω;R3),p1+p′1=1,Wb1,p(Ω;R3):={ψ∈W1,p(Ω;R3),ψ=0 on Ω′×{0}},
(Pε):⎩⎨⎧uε∈Wb1,p(Ω;R3)infFε(uε)−∫Ωf.uεdx,Fε(uε):=∫Ω∖Trεf(e(uε))dx+lε∫Trεg(e(uε))dx,e(uε)=21(∇uε+∇Tuε),f∈Lp′(Ω;R3),p1+p′1=1,Wb1,p(Ω;R3):={ψ∈W1,p(Ω;R3),ψ=0 on Ω′×{0}},
\displaystyle{\mathcal{W}}^{p}(\boldsymbol{a},{\boldsymbol{\alpha}};U,V):=\Big{\{}{\boldsymbol{\psi}}\in W^{1,p}_{0}(V;\mathbb{R}^{3}),{\boldsymbol{\psi}}(y)=\boldsymbol{a}+\tfrac{2}{{\rm diam}U}{\boldsymbol{\alpha}}\wedge\boldsymbol{y}\hbox{ in }U\Big{\}}.
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.
Taxonomy
TopicsAdvanced Mathematical Modeling in Engineering · Composite Material Mechanics · Advanced Numerical Methods in Computational Mathematics
Full text
Homogenization of Norton-Hoff fibered composites with high contrast
Michel Bellieud
LMGC (Laboratoire de Mécanique et de Génie Civil de Montpellier)
UMR-CNRS 5508, Université Montpellier II, Case courier
048, Place Eugène Bataillon, 34095 Montpellier Cedex 5, France
We study the steady creep flow of a
perfectly viscoplastic
solid reinforced by
fibers with high viscosity contrast.
Our study unveils new effects
related to anisotropy and conditioned by the Norton exponent.
Key words and phrases:
Homogenization, fibered structure, Norton-Hoff materials, visco-plasticity, anisotropic linear elasticity
Given a bounded Lipschitz cylindrical domain Ω:=Ω′×(0,L) of R3 and two strictly convex functions f, g satisfying a growth condition of order p∈(1,+∞), we consider the problem
[TABLE]
where Trε is
a distribution of disjoint cylinders of very small volume fraction
and lε
a large parameter.
The solution uε to (1.1) represents
the
Eulerian velocity field in a
Norton-Hoff material
of Norton exponent p−11 undergoing a steady
creep flow
under the influence of a density f of
applied body forces [38, 42, 46].
Composites comprising a small volume fraction of fibers with strong properties have been studied in various contexts
[6, 12, 16, 18, 23, 24, 36, 40, 42, 43, 48].
They are characterized by the interaction of concentration phenomena
in the fibers and in a small region of space surrounding them.
The equilibrium problem in linear elasticity, a special case of Problem 1.1 corresponding to positive definite quadratic forms f and g,
has been investigated in the isotropic case for fibers of circular cross-sections
in [15].
Our study unveils new effects
related to anisotropy and conditioned by the growth parameter p and the shape of the cross-sections of the fibers, described by some bounded connected open subset S of R2.
The distinctive feature of Problem 1.1 lies in the interaction of concentrations of stress in the fibers and rate of deformation in their close outer neighborhood.
For ease of exposition, we assume that the fibers are ε-periodically distributed. Our analysis goes through in the non-periodic case provided the fibers are well separated.
We show that the asymptotic behavior of uε depends on the order of magnitude of lε, characterized
in terms of the size rε of the cross-sections of the fibers
by the parameters
[TABLE]
and on the effective p-capacity of their cross-sections in Ω′, represented
by
[TABLE]
When 0<κ<+∞, we establish that the fibers locally behave like rigid bodies rotating around their principal axes (parallel to e3)
with an angular velocity δ,
the axes moving at the velocity v.
The field v approximates as ε→0 to the local average vε of uε over the cross-sections of the fibers, and δ
to their mean angular velocity δε.
We show that rεvε3 converges to some function w and demonstrate that
the effective contribution of the fibers
is described by a functional of v, δ, and w
of
the form
[TABLE]
If 0<k<+∞, we prove that the fibers possibly display a larger angular velocity, characterized
by a function θ approximating to θε:=rεδε. Tangential and normal velocities
are then of the same order at the surface of the fibers.
We show that
the concentration of stress in the fibers leads to a contribution of the type
[TABLE]
In the other cases, we establish that Φfibers vanishes. Setting
[TABLE]
we demonstrate that
the limit problem associated with (1.1) takes the form
[TABLE]
for some suitable domain D.
The second term of Φ in (1.7) stems from a concentration
of rate of deformations
in the close outer neighborhood of the fibers, resulting from the interaction between the matrix and the fibers.
We prove that
[TABLE]
for some convex function cf
satisfying, if p=2, the growth condition
[TABLE]
with c,C>0.
A
phenomenon
related to the Stokes’ paradox
induces a different behavior when p=2: we show that then,
[TABLE]
This means that the interaction between the matrix and the fibers precludes large angular velocities
of the fibers.
If γ(p)=+∞,
which always holds when p>2,
the limit problem is simply obtained by setting u=v
and θ=0 in (1.7).
We turn to a more detailed description of the mappings ghom and cf, using standard notations recalled in the next section.
We establish that the functions
ghom in (1.4)
and (1.5)
are, up to a multiplicative constant, the infima with respect
to q over W1,p(S;R3) of
[TABLE]
and
[TABLE]
respectively, where g0,p denotes a p-homogeneous approximation of g near [math].
On the other hand, assuming without loss of generality that ∫Sydy=0, we prove that
for every p∈(1,+∞),
[TABLE]
for some Rε≫rε, where D is the unit ball of R2. The mapping capf is the variant of the notion of capacity introduced in [8] (see also [51]) defined, for any (a,α)∈(R3)2 and any couple (U,V) of open subsets of R2
with U connected, bounded, ∫Uydy=0, and
U⊂V, by
[TABLE]
The extended real cf(a,ζ) can be seen as a capacity density
approximating to the sum of the images of the sections of the fibers under
capf(a,ζe3;.,Ω′) per unit surface.
Our investigation into
the properties of capf leads to the formula
[TABLE]
where f∞,p is a p-homogeneous approximation of f at ∞. If
p=2, we show that cf is independent of the shape of the cross-sections of the fibers.
An interesting feature of our results lies in the dependence of the limit problem on the effective rescaled angular velocity θ.
A non-vanishing θ may only arise when γ(p)<+∞, p<2, and 0<k<+∞.
This explains why this dependence was not detected in [15].
It arose in another context in [7].
One can guess, from (1.14) and (1.15), that θ is conditioned
by the shape of the cross-sections of the fibers:
the matrix is more likely to induce a rotating motion in fibers of pear-shaped cross-sections than circular ones.
Formulae (1.12), (1.14) and (1.15) suggest that
such rotations can also be brought about by the anisotropy of either matrix or fibers.
Besides, we prove that
θ can be influenced by
large twisting body forces applied on the fibers.
Another distinctive aspect of our work lies in
the dependence of Φfibers on the effective angular velocity δ and microscopic longitudinal velocity w when 0<κ<+∞. This was not perceived in the setting of linear isotropic elasticity
because, as we show, these functions vanish when the material constituting the fibers is linear isotropic, whatever the shape of their cross-sections.
By contrast, we prove that δ does not vanish,
in general, for anisotropic fibers. The same is likely to hold for w.
We also establish that, like for θ, large twisting body forces applied on the fibers may have an effect on δ and w and, in particular,
lead to a non-vanishing couple (w,δ) for isotropic fibers.
Notice that,
unlike θ,
the matrix exerts no influence on (δ,w).
Taking advantage of the above study, we next examine the problem
[TABLE]
when Tε is an ε-periodic
distribution of disjoint cylinders of volume fraction of order 1.
High-contrast homogenization problems of this type have been and are intensively investigated in many contexts
[1, 3, 4, 5, 19, 20, 28, 26, 25, 31, 39, 47, 50, 52, 53, 54].
Problem 1.16
has been studied
in detail in [7] in the setting of
linear isotropic elasticity for fibers of circular cross-sections, correcting results previously obtained in [13] where the influence of θ had failed to be taken into account.
We show that the effective problem associated to (1.16) takes the form
[TABLE]
The second term of Φsoft is common with
(1.7). The first one, which unlike (1.8) emanates from large rates of deformation arising in the entire matrix,
is given by
[TABLE]
where W♯1,p(Y;R3) denotes the set of Y-periodic members of Wloc1,p(R2;R3).
The paper is organized as follows: our main results are presented in Section 3 in the periodic case.
In Section 4,
we discuss their extension to a non-periodic or random setting
and
the case of large applied body forces.
Section 5 is devoted to the asymptotic analysis of the sequence of the solutions to (1.1)
and of some auxiliary sequences characterizing the behavior of the fibers. Section 6
comprises a detailed study of the mapping
capf on which our proofs crucially rely.
The demonstrations of our main results, based on the Γ-convergence method [33], are situated in Section 7.
The appendix comprises two technical lemmas relating to the lower bound and the proof of the convergence (1.13) in the case p=2.
Our results were partially announced in [9].
2. Notations
In this paper, {e1,e2,e3} stands for the canonical basis of R3.
Points in R3 and real-valued functions are represented by symbols beginning with a lightface minuscule (example x,i,detA...), vectors and vector-valued functions by symbols
beginning with a boldface minuscule (examples: i, u,
f, g,
divΨ,…),
matrices and matrix-valued functions by symbols beginning with a boldface majuscule
with the following exceptions:
∇u (velocity gradient), e(u) rate of strain tensor). The symbol In represents the n×n identity matrix.
We denote by ui or
(u)i the components of a vector u and by Aij or (A)ij those of a matrix A (that is u=∑i=13uiei=∑i=13(u)iei; A=∑i,j=13Aijei⊗ej=∑i,j=13(A)ijei⊗ej). We do not employ the usual repeated index
convention for summation. We denote
by A:B=∑i,j=13AijBij the inner product of two matrices,
by εijk the three-dimensional alternator, by
u∧v=∑i,j,k=13εijkujvkei the exterior product in R3, by Sn
the set of all real symmetric matrices of order
n, by ♯A the cardinality of a finite set A, by B (resp. D) the open unit ball of R3 (resp. R2), by suppf the support of a function f,
by R the space of rigid motions in dimension 2 or 3.
For any two vectors a, b in R3, we set a⊙b:=21(a⊗b+b⊗a). For any two symmetric matrices A, B, we write A≤B if B−A is semi-definite positive.
Given x∈R3, we write x=(x′,x3) for x′=(x1,x2).
For any two weakly differentiable fields ψ:R2→R3 and u:R3→R3, we set
[TABLE]
Given a topological space X, the symbol B(X) represents the σ-algebra of the Borel subsets of X.
For any Radon measure ν on X and any Banach space E, we denote by
Lνp(X;E) the set of E-valued Borel fields ψ on X such that ∫X∣ψ∣Epdν<+∞.
The letter
C (resp. the symbol C(a) if a dependence on some variable a is indicated) stands for different positive constants whose precise values may vary.
Let Ω:=Ω′×(0,L) be a bounded Lipschitz domain of R3, S a bounded Lipschitz domain of R2 verifying
[TABLE]
and (rε)ε>0 a sequence of real numbers such that
[TABLE]
We consider the ε-periodic distribution of fibers defined by (see fig. 1)
[TABLE]
where for any subset A of R2 and any bε>0,
[TABLE]
We are concerned with the homogenization of Problem 1.1 when f,g:S3→R
are strictly convex and satisfy
[TABLE]
for some p∈(1,+∞) and some positive constants c,C.
We prove that the limit problem depends on the parameters k, κ and γ(p) defined by (1.2), (1.3).
We focuse on the case
[TABLE]
Setting
[TABLE]
we suppose that
[TABLE]
for some ς, ϖ verifying
[TABLE]
Under these hypotheses, we show that the solution uε to (1.1)
converges in a sense defined below to u, the sequence (vε,θε) given
by
[TABLE]
where
[TABLE]
converges, up to a subsequence to
(v,θ), and, if κ>0, (rε1vε3,rε1θε), converges up to a subsequence, to (w,δ),
where
the couple (u,vtuple) defined by (1.6) is a solution to
[TABLE]
and u is the unique solution to
[TABLE]
The functional Φ is defined by
[TABLE]
where, if κ=0,
[TABLE]
if κ>0,
[TABLE]
and cf is defined by (1.13) in terms of capf
given by (1.14), and of any sequence (Rε) verifying
[TABLE]
We show that
cf is well defined, convex, independent of (Rε) and, if p≥2,
of S, and satisfies
Assume (3.2), (3.3), (3.5), (3.6), (3.8), then the solution uε to (1.1)
weakly converges in Wb1,p(Ω;R3) to u, the couple
(vε,θε) defined by (3.10) weakly converges in Lp(Ω;R3)×Lp(Ω), up to a subsequence, to
(v,θ), and, if κ>0, (rε1vε3,rε1θε) weakly converges
in (Lp(Ω))2, up to a subsequence, to (w,δ), where (u,vtuple) given by (1.6)
is a solution to (3.12) and u is the unique solution to (3.13).
Remark 1**.**
*(i)
The functional F in (3.13) is the Γ-limit of the sequence (Fε) in the weak topology of Wb1,p(Ω;R3) (see Remark 11).
(ii) If γ(p)=0 or k=0,
the solution uε to (1.1) weakly converges in Wb1,p(Ω;R3) to the unique solution to: infu∈Wb1,p(Ω;R3)∫Ωf(e(u))dx−∫Ωf⋅udx.
(iii) If f∞,p and g0,p are strictly convex, which the strict convexities of f and g do not ensure,
and p=2, then Φ
is strictly convex, the solution
to (3.12) unique, and all convergences stated in Theorem 1 hold for the whole sequences.
(iv) Linear isotropic elasticity. If f=fλ0,μ0 and g=fλ1,μ1 for some λ0,λ1≥0, μ0,μ1>0, where*
*where m:=infφ∈H1(S)∫−S(∂y1∂φ−y2)2+(∂y2∂φ+y1)2dy.
This infimum being attained, it results from the Schwarz theorem that m>0. Noticing that, by (3.23), the infimum (3.12) is achieved by
(w,δ)=(0,0) and that, by (3.18), θ=0, we
recover the
formulae
obtained for S=D in [15, Remark 2.2].
(v) The following example shows that δ
doesn’t vanish, in general, when anisotropic fibers are considered: assume that 0<κ<+∞ and*
[TABLE]
Then, a solution to the minimization problem (3.16) is q:=(ζ1φ(1)+ζ2φ(2)+aφ(3)+β~φ(4))e3, where β~:=diamS1β and
φ(1),..,φ(4) are solutions in H1(S)/R
to the Neumann problems (n: outer unit normal to ∂S)
[TABLE]
respectively.
We deduce that ghom(ζ1,ζ2,a,β)=ζ1ζ2aβ~⋅Cζ1ζ2aβ~,
where
[TABLE]
If (u,vtuple) is a solution to (3.12), then (w,δ) is a solution to
inf(w,δ)Φ(u,v,0,w,δ).
We infer
[TABLE]
If S=D, we obtain φ(4)=0,
C14=C34=0, C44=4C24>0, and deduce
δ=−2diamS∂x3∂v2.
we prove that the effective problem associated to (1.16) is
[TABLE]
where csoftf, ghom and D are defined by (1.17), (3.15), (3.16) and (3.19), respectively.
Theorem 2**.**
Assume (3.24) and (3.25), then
the solution uε to (1.16)
weakly converges in Lp(Ω;R3), up to a subsequence, to u, and the convergences
of (vε,θε) and (ε1vε3,ε1θε) stated in Theorem 1 hold,
where (u,vtuple)
is a solution to (3.26).
Remark 2**.**
(i) Theorems 1 and 2 can be generalized to multiphase media comprising non-intersecting families of parallel fibers distributed in different directions (for Theorem 2, case p=2, see [7, Sec. 4]).
(ii) The condition (3.25)
is guaranteed when κ>0 or when suitable multiphase media described in **[7, Sec. 5]** are considered, for instance
when k>0 and the fibers are distributed in at least three independent directions.
(iii) Most of Remark 1 applies to Problem 1.16.
4. Variants
4.1. Non-periodic case
The study of Problem 1.1
can be extended to a non-periodic setting:
we then parametrize
it by
the size r of the cross-sections of the fibers.
The collection of fibers is described in terms of the image Gr of the set of their principal axes under the orthogonal projection onto Ω′,
by setting
[TABLE]
We assume that (Gr)r>0 is a family of finite subsets of Ω′ satisfying, as r→0,
[TABLE]
We focuse on the case p∈(1,2].
To take advantage of the notations introduced in the periodic case, we fix an arbitrary positive real number γ(p) and introduce the small parameter εr defined by
[TABLE]
By (1.3) we have γεr(p)(r)=γ(p) for every r>0.
Given a sequence of positive real numbers (lεr)r>0, we consider the sequence of problems (Pr(Gr))r>0 formally deduced from (1.1) by substituting Tr(Gr) for Trε and lεr for lε.
The function defined by
[TABLE]
where Iεr and Yεri are given by (3.3),
locally approximates to the number of fibers crossing a square of size εr in the cross-section of Ω.
We suppose that (nGr)r>0 is bounded in L∞(Ω′) and
weakly⋆ converges to some n as r→0. Setting μ:=nL⌊Ω′2,
by combining the argument of the proof of Theorem 1 with the one developed in [16],
one can prove that
the effective problem takes the form
[TABLE]
where cf and ghom are given by (3.15), (3.16) and (3.18) in terms
of γ(p) introduced above and k, κ defined by
substituting εr for ε in (1.2). The domain
Dμ is deduced from (3.19) by substituting
Lμ⊗L1p(Ω;R3) for Lp(Ω;R3), the spaces
Lμp(Ω′;E) for Lp(Ω′;E) for any Banach space E, and by replacing the homogeneous Dirichlet conditions on
Ω′×{0} by homogeneous Dirichlet conditions μ⊗δ{0}-a.e. on
Ω′×{0}.
Theorem 3**.**
Under the assumptions stated above,
the solution ur to (Pr(Gr))
weakly converges
in Wb1,p(Ω;R3) as r→0 to the unique solution u to (4.5).
Remark 3**.**
(i) Under the same assumptions,
one can show that the measure m~r, defined by substituting Tr(Gr) for Srε×(0,L) and (r,εr) for (rε,ε)
in
(5.9), weakly⋆ converges to μ⊗L1 in M(Ω),
and the
sequences (v~rm~r), (θ~rm~r) and, if κ>0, (εr1v~r3m~r) and (εr1θ~rm~r),
where
[TABLE]
*weakly⋆ converge up to a subsequence in M(Ω;R3) and M(Ω) to vμ⊗L1, θμ⊗L1, wμ⊗L1, and δμ⊗L1, respectively,
where vtuple, defined by (1.6), is a solution to the minimization problem defining Fμ(u) in (4.5).
(ii) Theorem 3 holds, in some cases, when (nGrL⌊Ω′2) is
bounded in M(Ω′) and
weakly⋆ converges to some singular
measure μ which vanishes on all Borel subset of Ω′ of p-capacity zero. For instance,
assume that
Ω′=(0,L)2 and set*
[TABLE]
*The set Tr(Gr) defined by (4.1) represents a family of fibers whose principal axes are εr2-periodically distributed on the surface
Σ:={L/2}×(0,L)2. The sequence (nGrL⌊Ω′2) defined by (4.4) is bounded in M(Ω′)
and weakly⋆ converges to μ:=H⌊{L/2}×(0,L)1.
One can show, by
combining the argument of the proof of Theorem 1 with that developed in [14],
that Theorem 3 holds, as well as
the convergences stated in (i).
(iii) Dirichlet problems in varying domains. Under the assumptions of Theorems 1 or 3 (or Remark 3 (ii)), the solution to*
[TABLE]
weakly converges in Wb1,p(Ω;R3) to the unique solution to
[TABLE]
deduced from (Phom,μ)
by substituting [math] for vtuple.
The second term is analogous to the so-called ”strange term”
[30, 43].
If (nGr)r>0 is only assumed to be bounded in L1(Ω′) and to weakly⋆ converge to some arbitrary
measure μ, we expect the limit problem associated with (4.6) to depend, not on μ, but, up to a subsequence, on some
measure
μ0
defined
through a variant of the γ-convergence introduced in [34].
Many references on this subject can be found in
[35].
In particular, compactness results ensuring the existence of a γ-converging subsequence have been established in various contexts.
Under (4.2), the limit problem associated with (4.6) is likely to be deduced from (4.7), formally, by substituting such a measure μ0 for μ.
This suggests that the limit problem associated with (Pr(Gr)) might possibly be (Phom,μ0).
(iv) The assumption (4.3) on the choice of εr was omitted in **[16]**.
In order the results stated in **[16]** to be correct, one should either assume in **[16, (17)]** that γ(p) is finite (hence 1<p≤2), or in **[16, (5)]** that nε is bounded from below by a positive constant, otherwise the asymptotic behavior of the sequence (γε(r)nε), whose knowledge is necessary to obtain for instance **[16, (126)]**, would be undetermined.
4.2. Random case
Assuming 1<p≤2, fixing d>0, we set
[TABLE]
One can check that the finite metric dO:(ω,ω′)∈O2→inf{1,dH(ω,ω′)}, where dH stands for the Hausdorff distance, turns O into a complete metric space.
Denoting by B(O) the associated Borel σ-algebra on O, we consider a probability measure P on (O,B(O)) satisfying
[TABLE]
We introduce the σ-algebra F and the random variable n0 on O defined by
[TABLE]
We denote by EPFn0
the conditional expectation of n0 given F w.r.t. P.
We fix a positive real number γ(p) and set
[TABLE]
where εr is given by (4.3).
The following statement is proved in [10, Th. 2.4.2]:
Theorem 4**.**
Under the assumptions stated above, there exists a sequence (rk)k∈N converging to [math] and a P-negligible set N∈B(O) such that, for each ω∈O∖N, the sequence (nGrk(ω)) defined by (4.4) weakly⋆ converges in L∞(Ω′) to
the constant function EPFn0(ω).
The assumptions (4.3) and (4.8) ensure that Gr(ω) satisfies (4.2).
By combining Theorems 3 and 4, we obtain:
Corollary 1**.**
There exists a sequence (rk)k∈N converging to [math] and a P-negligible set N∈B(O) such that,
for all ω∈O∖N, the solution urk(ω)
to the problem
(Prk(Grk(ω))) considered in Section 4.1
weakly converges in Wb1,p(Ω;R3)
to the unique solution u(ω) to
(Phom,μ(ω)) defined by (4.5), where
μ(ω):=EPFn0(ω)L⌊Ω′2.
4.3. Large applied body forces
We consider the problems
[TABLE]
[TABLE]
where, given f∈Lp′(Ω;R3), g0∈C(Ω×S;R3) and a0,β0∈C(Ω×S),
fε is defined in terms of yε given by (3.11), by
[TABLE]
with
[TABLE]
We establish that the limit problems associated to
(4.9) and (4.10) are
The statements deduced from Theorems 1 and 2
by substituting
(4.9), (4.12) for
(1.1), (3.12) and (4.10), (4.13) for
(1.16), (3.26) hold.
Remark 4**.**
When κ>0 and Φfibers is given by (3.23),
if (u,vtuple) is a solution to (4.12) or (4.13), then
[TABLE]
5. Preliminary results and a priori estimates
This section is devoted to the study of the asymptotic behaviors of a sequence (uε) satisfying
[TABLE]
and of the auxiliary sequences (vε), (θε), (rεvε3), (rεθε)
defined by (3.10). Our main results, stated in Section 5.5,
will be deduced from a series of inequalities proved in Section 5.4.
5.1. Auxiliary sequences
The proof of Theorem 1 rests on
a choice of sequences of auxiliary fields
•
weakly differentiable w.r.t. x3,
•
weakly relatively compact in Lp,
•
locally characterizing a rigid motion
approximating to the velocity in the fibers.
The selection of surface integrals in (3.10) is required to ensure the first condition, and is a hindrance to the third.
To circumvent this difficulty, we introduce another couple of auxiliary sequences
constructed from volume instead of surface integrals,
satisfying the second and third conditions and having the same cluster points as vε, θε.
For convenience, we rewrite
(3.10) as follows
[TABLE]
where, for any Lipschitz domain A of R2 verifying
[TABLE]
\textswabvεA and θεA
are the linear operators defined on Lp(Ω;R3) by (see (3.3), (3.4), (3.11)):
[TABLE]
The ”volumic” couterparts \textswabvε(uε) and θε(uε) of (3.10) will be defined by splitting
each fiber Trεi into small cylinders Trεij of size rε given by
(see fig. 2)
[TABLE]
and setting for every φ∈Lp(Ω;R3),
[TABLE]
where φ is an extension of φ
(see (5.47))
bringing meaning to
∫−Trεijmφ
where
[TABLE]
We will prove that the piecewise rigid motion \textswabrε(uε)
approximates to the velocity in each set Trεij and asymptotically behaves
as \textswabrεS(uε).
To derive formula (1.13),
we will employ the auxiliary field
\textswabuε(uε)
defined by
[TABLE]
in terms of a sequence (Rε) satisfying (3.17). The assumption (3.17)
ensures that (\textswabuε(uε)) approximates to uε in Ω.
The selection of integrals over DRεi∖DRε/2i in (5.8)
will yield a crucial estimate at a technical step.
5.2. Two-scale convergence with respect to (mε)
To particularize the effective behavior of uε in the fibers,
we introduce variant of the two-scale convergence [1, 29, 45]
defined in terms of the measures
[TABLE]
which are supported on the fibers and weakly⋆ converging in M(Ω) to L⌊Ω3.
We say that a sequence
(φε) in Lp(Ω;Rk)
two-scale converges to φ0∈Lp(Ω×S;Rk) with respect to (mε) if
(see (3.11))
[TABLE]
The main properties of this convergence are stated in the following Lemma.
Lemma 1**.**
For every η∈C(Ω;L∞(S;Rk)),
[TABLE]
Any sequence (φε) in Lp(Ω;Rk)
such that
[TABLE]
two-scale converges w.r.t. (mε), up to a subsequence, to some
φ0.
Setting φ(x):=∣S∣1∫Sφ0(x,y)dy, the following implications holds:
[TABLE]
If φε⇀⇀\buildrelmεφ0, for any convex function j:Rk→R, we have
[TABLE]
Proof. The proof of (5.11) is similar to that of [1, Lemma 1.3]. Setting ∫η⋅dνε:=∫φε(x)⋅η(x,rεyε(x′))dmε(x)∀η∈C(Ω×S;Rk),
the
two-scale convergence w.r.t. (mε) of (φε) to φ0
is equivalent to
the weak⋆ convergence in M(Ω×S) of (νε) to ∣S∣1φ0(x,y)L⌊Ω×S5.
By (5.12) and Hölder’s inequality, we have
[TABLE]
thus
(νε), bounded
in M(Ω×S;Rk),
weakly⋆ converges up to a subsequence to some ν.
Passing to the limit in (5.16), taking
(5.11) into account, we obtain ∫η⋅dν≤C(∫Ω×S∣η(x,y)∣p′dxdy)p′1, hence, by
the Riesz representation theorem
ν=∣S∣1φ0(x,y)L⌊Ω×S5 for some φ0∈Lp(Ω×S;Rk).
The assertion (5.13) is straightforward.
By choosing in (5.10) test fields independent of y,
we obtain (5.14).
Denoting by j⋆ the Fenchel transform of j,
we deduce from Fenchel’s inequality, (5.10) and (5.11), that
[TABLE]
Noticing that by the convexity of j we have j⋆⋆=j, we infer
[TABLE]
the second line being justified in Remark 7.
The second inequality in (5.15) results from (5.14) and Jensen’s inequality. ∎
Remark 5**.**
Combining Lemma 1 with [11, Lemma 4.2], one can derive
a non-periodic version of Lemma 1 adapted for the proof of Theorem 3.
5.3. Properties of the auxiliary sequences
In this section, we compare diverse types of convergence for the auxiliary sequences (\textswabpε(uε)) introduced in Section 5.1. It turns out that, except for
\textswabpε∈{\textswabrε,\textswabrεS}, the weak limit of (\textswabpε(uε)) in Lp,
its two-scale limits w.r.t. to
(mε) and the weak* limits of (\textswabpε(uε)mε) in M(Ω) are the same, and that (\textswabvεS(uε)) and (\textswabvε(uε)) have
the same weak limits in Lp.
We will prove in Proposition 3 that the same holds
for (θεS(uε)e3) and (\textswabwε(uε)). This, combined with (5.23), shows that (\textswabrεS(uε)) and (\textswabrε(uε)) have the same two-scale limits w.r.t. (mε).
Proposition 1**.**
*Let \textswabpε be a linear combination of
\textswabwε, \textswabvε, θεe3, \textswabvεA, θεAe3, \textswabuε
defined by (5.4),
(5.6),
(5.8), and (φε) a sequence in Lp(Ω;R3) satisfying (5.12).
*
(i) For all φ∈W1,p(Ω;R3), ε>0, k∈{1,2,3}, we have
[TABLE]
[TABLE]
(ii) For every p∈Lp(Ω;R3), the following equivalences hold
[TABLE]
(iii) If (\textswabpε(φε)) two-scale converges to p0 w.r.t. (mε), then
p0(x,y)=p(x) a. e. in Ω×S for some p∈Lp(Ω;R3).
Furthermore, for every p∈Lp(Ω;R3),
[TABLE]
(iv) If \textswabpε is a linear combination of
\textswabwε, \textswabvε, θεe3, then
[TABLE]
where zε(x) is defined by (5.5).
In particular, the following holds:
[TABLE]
Proof. (i) Since \textswabpε(φ) is constant in each set Yεi×{x3}, by (3.3) and (5.9),
[TABLE]
We deduce from (5.3), (5.4), (5.9) and Jensen’s inequality, that
for k∈{1,2,3},
[TABLE]
The inequality ∫∣\textswabpε(φ)∣pdmε≤C∫Ω∣φ∣pdmε
is obtained in a similar way. ∎
As ψ is compactly supported in Ω, the following estimate holds:
[TABLE]
Taking (5.4), (5.6), (5.8), (5.9), and the constant nature of ψ(ε) in each set Yεi×Δεj into account, elementary computations yield
[TABLE]
By (5.12), (5.19) and
Lemma 1,
the sequences (φεmε), (\textswabvεS(φε)mε), and (\textswabvε(φε)mε) weakly⋆ converge in M(Ω;R3), up to a subsequence,
to p1, p2, and p3, respectively, for some pi∈Lp(Ω;R3).
By passing to the limit as ε→0 in (5.26), taking (5.25) into account,
we infer ∫Ωψ⋅p1dx=∫Ωψ⋅p2dx=∫Ωψ⋅p3dx and deduce from the arbitrariness of ψ that
p1=p2=p3. ∎
(iii) Assume that \textswabpε(φε)⇀⇀\buildrelmεp0 for some p0∈Lp(Ω×Y;R3). Since
\textswabpε(φε) is constant in each set Yεi×{x3}, for any
Ψ∈D(Ω×S;R3×R3) we have \displaystyle\int{{\textswab{\boldsymbol{p}}}}_{\varepsilon}({\boldsymbol{\varphi}}_{\varepsilon})\cdot{\rm\bf div}_{y}\boldsymbol{\Psi}\Big{(}x,\frac{y_{\varepsilon}(x^{\prime})}{r_{\varepsilon}}\Big{)}dm_{\varepsilon}=0. Passing to the limit as ε→0, we get
∣S∣1∫Ω×Sp0⋅divyΨdxdy=0
and deduce from the arbitrary choice of Ψ that p0(x,y)=p(x) a.e. in Ω×S, where p(x):=∫−Sp0(x,y)dy∈Lp(Ω;R3). It follows from (5.14) that
⇀pL3\textswabpε(φε)mε\buildrel⋆ weakly⋆ in M(Ω;R3).
On the other hand, by (5.12), (5.18), and (5.19), (\textswabpε(φε)) is bounded in Lp(Ω;R3) and weakly converges, up to a subsequence, to some p~.
By passing to the limit as ε→0 in (5.27), taking (5.25) into account, we obtain ∫Ωψ⋅pdx=∫Ωψ⋅p~dx and deduce p=p~. Conversely, if \textswabpε(φε)⇀p weakly in Lp(Ω;R3), then by (5.12), (5.19) and Lemma 1, (\textswabpε(φε)) two-scale converges w.r.t. (mε), up to a subsequence, to some
element of Lp(Ω;R3) which, by virtue of the last established implications, necessarily equals p.
The equivalences stated in (5.21) are proved and the implication is obtained
by choosing y∧η(x,y) as a test function for the two-scale convergence of (\textswabpε(φε)) to p w.r.t. (mε).
(iv)
Assume that \textswabpε is a linear combination of
\textswabwε, \textswabvε, θεe3 and \textswabpε(φε)⇀p weakly in Lp(Ω;R3). By (5.12) and (5.19), ∫\textswabpε(φε)∧rεzεpdmε≤C, hence, by Lemma 1, (\textswabpε(φε)∧rεzε)⇀⇀\buildrelmεq0,
up to a subsequence, for some q0.
We fix η∈D(Ω×S;R3) and set
[TABLE]
The fields \textswabpε(φε) and
η[ε](x,rεyε(x′)) are independent of x3 in each Yεi×Δεj and, by (3.11) and (5.5), ∫Δεjzε(x′,x3)dx3=∫Δεjyε(x′)dx3,
therefore
[TABLE]
Summing w.r.t. (i,j) over Iε×(Lε∖{jεm}), passing to the limit as ε→0, taking
the estimate ∣η−η[ε]∣L∞(Ω×S;R3)≤Cε into account and noticing that, by (5.21), \textswabpε(φε)∧rεyε(x′)⇀⇀\buildrelmεp∧y, we obtain
∣S∣1∫Ω×Sq0(x,y)⋅η(x,y)dxdy=∣S∣1∫Ω×S(p∧y)⋅η(x,y)dxdy
and deduce q0=p∧y. The assertion (5.22) is proved.
The assertion (5.23) results from (5.6), (5.21), and (5.22). ∎
5.4. Key inequalities
In this section, we establish a series of inequalities which, combined with Proposition 1,
will yield
a number of a priori estimates and convergences for a sequence satisfying (5.1) and its associated auxiliary sequences.
The following Korn’s inequalities, proved in [37, 41], will be employed at several occurrences:
Lemma 2**.**
(i) We have, for N∈{2,3},
[TABLE]
(ii) If U is a bounded Lipschitz domain of RN and V a subspace of W1,p(U;RN)
such that V∩R={0}, where R is the space of
rigid motions in U, then
[TABLE]
We set
[TABLE]
Proposition 2**.**
Let A be a bounded Lipschitz domain of R2 satisfying (5.3) and let α∈{1,2}.
The following inequalities hold for every φ∈W1,p(Ω;R3):
[TABLE]
The following inequalities hold for every φ∈Wb1,p(Ω;R3):
[TABLE]
[TABLE]
where
[TABLE]
Proof.Proof of (5.32), (5.33), (5.34).
By (5.5), we have (see fig. 2)
[TABLE]
We consider the linear
operators \textswabr,\textswabv,\textswabw defined on W1,p(T;R3) by
[TABLE]
We have \textswabr∘\textswabr=\textswabr and \textswabr(φ)=φ∀φ∈R, hence V:={φ−\textswabr(φ),φ∈W1,p(T;R3)} satisfies
V∩R={0}.
Applying Lemma 2, noticing that e(φ)=e(φ−\textswabr(φ)),
we infer
[TABLE]
By making appropriate changes of variables, taking (5.6) and (5.44) into account, we deduce that for every (i,j)∈Iε×Lε,
[TABLE]
If φ∈W1,p(Ω;R3), its extension
defined on Ω′×(L,2L) by
yielding (5.32).
Applying a similar argument to
W:={φ∈W1,p(T;R3),\textswabv(φ)=0},
we obtain
∫∣φ−\textswabvε(φ)∣pdmε≤Crεp∫∣∇(φ)∣pdmε and, taking (1.3), (5.9) and (5.31) into account, deduce
(5.33).
To prove (5.34), we start from the inequality
∣b∣p≤C∫T∣b∧z∣pdz∀b∈R3 (easily proved by contradiction).
By suitable changes of variables, we obtain
[TABLE]
Substituting for b the constant value taken by
\textswabwε(φ) in each set Trεij,
we infer
[TABLE]
which, combined with (5.32) and (5.33), yields (5.34). ∎
Proof of (5.35), (5.36). We put Yε:=⋃i∈IεYεi.
By (5.4) and (5.8), we have
[TABLE]
By (3.3), Ω′∖Yε⊂{x′∈Ω′,dist(x′,∂Ω′)≤2ε}, thus, since Ω′ is Lipschitz, L2(Ω′∖Yε)≤Cε.
By Hölder’s and Poincaré’s inequalities and the continuous embedding of W1,p(Ω;R3) into Lq(Ω;R3) (see (5.31),
[22, Corollary 9.14]), the following inequalities hold in Wb1,p(Ω;R3):
[TABLE]
Let E be a bounded Lipschitz domain of R2 such that ∂D⊂E.
We prove below the existence of C>0 such that, for every α∈(0,1] verifying αE⊂Y,
[TABLE]
Let us see how the claim follows from (5.52):
let (dε)⊂R be such that 0<dε≪ε and Edεi defined by (3.4).
An elementary change of variables yields
[TABLE]
Summing w.r.t. i over Iε and integrating w.r.t. x3 over (0,L),
choosing successively (E,dε):=(D∖D/2,Rε) and (E,dε):=(A,rε),
noticing that, by (1.3) and (5.52), εpl(εdε)≤γε(p)(dε)1, we infer
[TABLE]
Taking (5.50) and (5.51) into account, (5.35) and (5.36) are proved. We turn to the proof of (5.52).
A straightforward variant of [12, Lemma A4] yields
Applying (5.55)
to ψ:=φ′(.,x3), (5.56) to η:=φ3(.,z3), integrating w.r.t. z3 over (−1,1), noticing that ∣e(ψ)∣+∣∇η∣≤C∣ex′(φ)∣,
we infer
[TABLE]
and deduce, by suitable changes of variables,
[TABLE]
Noting that,
by (5.47),
∫Trεijεm∣ex′(φ)∣pdx≤C∫Trεi∣ex′(φ)∣pdx,
we infer (5.37).
∎
Proof of (5.38). Fixing φ∈W1,p(Ω;R3), applying (5.37) to
φ′,
observing that ∣ex′(φ′)∣≤∣e(φ)∣ and \textswabrεA(φ′)=\textswabrεA(φ)′ (see (5.4)), we obtain
[TABLE]
Combining this with the following inequality, deduced from (5.6) and (5.32),
By (5.4) and (5.6), θεA((\textswabrε−\textswabrεA)(φ)′)=(θε−θεA)(φ) and
\textswabvεA((\textswabrε−\textswabrεA)(φ)′)=(\textswabvε−\textswabvεA)(φ))′. This, along with the above inequalities, proves (5.38). ∎
Proof of (5.39).
Given f∈W1,p(T), we set g(t):=−∫Af(x′,t)dH2(x′) and TA:=A×(−1,1). For every x3∈(−1,1)2, we have
[TABLE]
This, combined with the following Poincaré-Wirtinger inequality
[TABLE]
implies
[TABLE]
By suitable changes of variables,
we infer (see (5.4), (5.6))
[TABLE]
Summing w.r.t. (i,j),
in view of (5.9), (5.31) and (5.47), we
obtain (5.39). ∎
Proof of (5.40). Let us fix φ∈W1,p(Ω;R3). By (5.4), for a. e. x3∈(0,L),
[TABLE]
By (3.11), ∂x2∂yε1=∂x1∂yε2=0 in Trεi
and
rεyε coincides on ∂Drεi×{x3} with the outward normal to ∂Drεi, hence
∫−Drεi−rεyε2∂x1∂φ3(.,x3)+rεyε1∂x2∂φ3(.,x3)dH2=0.
Summing this with (5.58), we obtain
[TABLE]
and deduce
∫∂x3∂θεD(φ)pdmε≤C∫∣e(φ)∣pdmε.
If φ∈Wb1,p(Ω;R3),
by (5.4) we have θεD(φ)∈Lp(Ω′;Wb1,p(0,L)) where Wb1,p(0,L):={η∈W1,p(0,L),η(0)=0}, therefore,
by (5.9) and Jensen’s inequality,
[TABLE]
Taking (5.38) into account, the assertion (5.40) is proved. ∎
Proof of (5.41).
The Banach space V:={ψ∈W1,p(S×(0,L);R3),ψ=0 on S×{0}}
satisfies V∩R={0}.
By applying (5.30)
to
ψ∈V defined
by ψα(z1,z2,z3):=φα(rε(z1−i1),rε(z2−i2),z3), α∈{1,2} and
ψ3(z):=rε1φ3(rε(z1−i1),rε(z2−i2),z3), and by making a suitable change of variables, we infer
∫Trεi∣rεφ1∣p+∣rεφ2∣p+∣φ3∣pdx≤C∫Trεi∣e(φ)∣pdx, yielding (5.41).
∎
In the next three propositions, we
establish a series of convergences
for a sequence (uε) satisfying (5.1), the sequence of its symmetrized gradients and the associated auxiliary sequences defined by (5.2), (5.8), (3.17), (5.43).
Proposition 3**.**
Let uε be a sequence in Wb1,p(Ω;R3)
satisfying (5.1).
The following convergences hold, up to a subsequence:
[TABLE]
where τ is defined by (3.10).
If in addition κ>0, then
Noticing that Wb1,p(Ω;R3)∩R={0}, we infer from (5.30) that
[TABLE]
Hence (uε) is bounded in Wb1,p(Ω;R3), thus weakly converges,
and strongly in Lp(Ω;R3), up to a subsequence, to some u. Since γ(p)>0, by (1.3) and (3.17), limε→0γε(p)(Rε)=+∞, therefore, by
(5.35),
(\textswabuε(uε)) strongly converges to u in Lp(Ω;R3).
Observing that, by (5.36),
(\textswabvεS(uε)) is bounded in Lp(Ω;R3), we infer from (5.18),
(5.33), (5.38), (5.39), and (5.63) that
By (5.21), up to a subsequence, each of the above sequences (\textswabpε(uε))
weakly converges in Lp(Ω;R3) to some p
and two-scale converges w.r.t. (mε) to the same p.
Taking (5.20), (5.23), (5.38) and (5.63) into account, we infer
[TABLE]
for some suitable v,w,θ. It follows from (5.32) and (5.63) that
and
deduce that w2=0. We likewise obtain w1=0, yielding
[TABLE]
By (5.13), (5.67) and (5.69), the sequence (−rεyε2uε1+rεyε1uε2) two-scale converges w.r.t. (mε) to
−y2v1+y1v2+diamS2∣y∣2θ, therefore, by (3.1) and (5.14),
⇀∫−S(−y2v1+y1v2+diamS2∣y∣2θ)dy=τθ(−rεyε2uε1+rεyε1uε2)mε\buildrel⋆ weakly⋆ in M(Ω).
The
assertion (5.59) is proved. ∎
hence by Lemma 1, (5.20), and (5.21),
the following convergences hold
[TABLE]
up to a subsequence, for some w, δ. By
(5.42) and (5.63), we have
[TABLE]
hence the sequence (rεsS(uε)mε) weakly⋆ converges in M(Ω), up to a subsequence, to some
g∈Lp(Ω).
We set
[TABLE]
One can check that
[TABLE]
We fix ζ∈C1(Ω), α,β∈{1,2}.
By (3.1) and (3.11), ∫Srεiyε(x′)dx′=0,
hence,
since vεβS(uε)ζε
is constant in each
Srεi×{x3}, ∫rεyεαvεβS(uε)ζεdmε=0. Taking (5.4)
and (5.43) into account, we deduce
[TABLE]
where τ is given by (3.10).
By passing to the limit as ε→0, thanks to (5.70) and (5.72), we obtain
∫Ωgηdx=τ∫Ωδηdx.
The
assertion (5.60) is proved. ∎
therefore, by (5.18), (5.66) and (5.68), θ=v3=0 if k=+∞,
w=0 and θ=0 if
p=2 or γ(p)=+∞. By (5.59), the sequence (\textswabuε(uε)−\textswabvεS(uε)) weakly converges in Lp(Ω;R3) to u−v, therefore, by
(5.35), (5.36) and (5.63),
[TABLE]
thus u=v if γ(p)=+∞.
If κ=+∞, by (5.40), (5.41), and (5.63), we have
[TABLE]
hence, by (5.59) and (5.70),
v=0 and δ=w=0.
Proposition 3 is proved.
∎
In the following proposition, we identify the two-scale limits w.r.t. (mε) of the sequence
(e(uε))
in terms of the functions v3 and θ given by (5.59).
Proposition 4**.**
Let (uε) be a sequence in Wb1,p(Ω;R3)
satisfying (5.1,5.59).
Then
[TABLE]
and, up to a subsequence,
[TABLE]
for some q∈Lp(Ω;W1,p(S;R3)), where L is defined by (3.15).
Proof. By (5.63) the sequence (e(uε)) satisfies (5.12), hence, by
Lemma 1,
[TABLE]
up to a subsequence for some Ξ0∈Lp(Ω×S;S3).
We fix ψ∈C∞(Ω;D(S)) such that ψ=0 on Ω′×{L}. By passing to the limit as ε→0 in the equation
By Hölder’s inequality, (5.32) and (5.63), we have
[TABLE]
Let us fix ξ∈R3.
Suitable changes of variables in the equations
0=−∫S(ξ1y1+ξ2y2)(∂y1∂φ1+∂y2∂φ2)(y)dy=∫Sξ⋅φ(y)dy, deduced from (5.77), yield
[TABLE]
Since
∂x3∂ηε\textswabvε(uε) is constant in each Srεi×{x3},
we infer
[TABLE]
By (3.11), (5.5), (5.6) and (5.77), the following holds: zε(x)=yε(x′)+zε3(x)e3, φ3=0 and wε3(uε)=θε(uε), therefore
[TABLE]
Since zε3, \textswabwε(uε), and ∂x3∂ηε are constant in
Srεi×{x3}, we deduce from (5.81) that
[TABLE]
By (5.21) and (5.66), we have θε(uε)mε⇀⇀\buildrelmεθe3∧y, hence, by (5.72),
Therefore ∫Ωθ∂x3∂ηdx≤C∣η∣Lp′(Ω)
and, by the arbitrariness of η,
[TABLE]
By (5.76) and (5.87), (5.73) is proved.
By integrating (5.86) by parts,
we find
[TABLE]
and, by the arbitrary choice of η∈D(Ω) and φ verifying (5.77), deduce that
[TABLE]
for some q3∈Lp(Ω;W1,p(S)).
We fix Ψ∈D(S;S3) such that Ψi3=0∀i∈{1,2,3} and divΨ=0. By (5.71), we have
∫e(uε):ηε(x)Ψ(rεyε(x′))dmε=0.
Passing to the limit as ε→0, in view of (5.72) and (5.75),
we obtain
∫Ω×SΞ0:η(x)Ψ(y)dxdy=0.
By the arbitrary choice of η and Ψ, we deduce from
a generalization of the Donati’s theorem [44, Th. 1] that
[TABLE]
for some q1,q2∈Lp(Ω;W1,p(S)).
Combining (3.15), (5.76), (5.88), and (5.89), we obtain Ξ0=\textgothL(q,∂x3∂v3,∂x3∂θ). Proposition 4 is proved. ∎
In the next proposition, we derive the two-scale limits w.r.t. (mε) of the sequences
(rεuε3) and
(rε1e(uε)) in terms of w,δ,v1,v2 given by (5.59) and (5.60), in the case κ>0.
Proposition 5**.**
Let (uε)⊂Wb1,p(Ω;R3)
satisfying (5.1), (5.59), (5.60). If κ>0,
[TABLE]
and, up to a subsequence,
[TABLE]
for some l∈Lp(Ω;W1,p(S;R3)), being I defined by (3.16).
Moreover,
[TABLE]
Proof. If κ>0, by (1.2), (5.41), and (5.63), we have ∫rεuε3p+rε1e(uε)pdmε≤C. Applying Lemma 1, taking (5.14) and (5.60) into account, we infer
[TABLE]
up to a subsequence,
for some w0∈Lp(Ω×S), Υ0∈Lp(Ω×S;S3).
We fix α∈{1;2} and ψ∈C∞(Ω;D(S)) such that ψ=0 on Ω′×{L}. By integration by parts,
[TABLE]
Passing to the limit as ε→0, taking
(5.93) into account and
noticing that, by (5.59) and (5.61), uε⇀⇀\buildrelmεv1e1+v2e2, we obtain
[TABLE]
We deduce from the arbitrariness of ψ that ∂x3∂vα∈Lp(Ω), ∂yα∂w0∈Lp(Ω×S), ∂yα∂w0(x,y)=−∂x3∂vα(x) a. e. in Ω×S and vα=0 a. e. in Ω′×{0}.
Thus ∇y(w0(x,y)+∂x3∂v1(x)y1+∂x3∂v2(x)y2)=0 and, by (3.1) and
(5.93),
[TABLE]
Passing to the limit as ε→0 in
∫rε1∂x3∂uε3ψdmε=−∫rεuε3∂x3∂ψdmε,
we obtain
∫Ω×SΥ033ψdxdy=−∫Ω×S(−∂x3∂v1y1−∂x3∂v2y2+w)∂x3∂ψdxdy, and infer
[TABLE]
Substituting (rεθε(uε)) for (θε(uε)), and rε1e(uε) for e(uε)
in the argument employed to establish (5.87), (5.88), (5.89),
we find
[TABLE]
for some l∈Lp(Ω;W1,p(S;R2)), and deduce Υ0=\textgothJ(l,∂x32∂2v1,∂x32∂2v2,∂x3∂w,∂x3∂δ) (see (3.16)).
The proofs of (5.90) and (5.91) are completed.
The assertion (5.92) follows from (3.19), (5.61), (5.73), and (5.90).
∎
6. Properties of capf
The proof of Theorem 1 relies in the apriori estimates and convergences established in Section 5 and an investigation into the properties of the capacity capf developed below.
In what follows, f:S3→R denotes a convex function, not necessarily strictly so, satisfying (3.5)
for some p∈(1,+∞) and some positive constants c,C.
For every nonempty bounded Lipschitz domain S⊂R2
verifying (3.1)
and every open set V⊂R2
such that S⊂V, we consider the mapping capf(.;S,V):(R3)2→[0,+∞) defined by (1.14).
The infimum problem Pf(a,α;S,V) may fail to be attained
when V is not bounded (see Remark 6).
We prove below that, if 1<p<2, capf is equivalently defined by a well posed minimum problem, namely
[TABLE]
where K0p(V;R3) is the closure of D(V;R3) in the reflexive Banach space
[TABLE]
By the Poincaré inequality, K0p(V;R3) is equal to W01,p(V;R3) when V is bounded in one direction. Otherwise, it
may be larger.
Lemma 3**.**
*(i) The functional φ→∫Vf(ey(φ))dy is
strongly continuous and weakly lower semi-continuous on W1,p(V;R3) and, if 1<p<2,
on Kp(V;R3).
(ii) Problem Pf(a,α;S,V) defined by (1.14) and, if 1<p<2, PKf(a,α;S,V) given by (6.1), have minimizing sequences in D(V;R3).
(iii)
If 1<p<2,
PKf(a,α;S,V) has a solution,
unique if f is strictly convex. Moreover, (6.1) holds and, for all φ∈K0p(V;R3),*
[TABLE]
*for some C(p)>0 independent of V.
(iv) If V is bounded in one direction and p∈(1,+∞), the same holds for Pf(a,α;S,V).
Moreover, for all φ∈W01,p(V;R3),*
[TABLE]
(v) There exists positive constants c(p,S,V), C(p,S,V) such that
[TABLE]
and, if 1<p<2 or V is bounded in one direction,
[TABLE]
Remark 6**.**
We prove below (see Remarks 8, 9) that, if p≥2, capf(a,0;S,R2)=0, therefore (6.6) doesn’t hold for V=R2 and the infimum
problem
Pf(a,0;S,R2) is not attained when a=0.
This is comparable to the Stokes’ paradox.
Proof. (i) By
(3.5),
the functional φ→∫Vf(ey(φ))dy is convex and bounded on the unit ball of W1,p(V;R3),
hence
strongly continuous, thus weakly lower semi-continuous.
The same holds with Kp(V;R3) in place of W1,p(V;R3) if 1<p<2.
(ii) Follows from (i) and a density argument as in [8, Lemma 1].
(iii) Given φ∈D(V;R3), applying (5.29)
for N=2 to φ′ defined by
(2.1) seen as an element of D(V;R2), noticing that ∣e(φ′)∣≤∣ey(φ)∣ and
∣∇φ3∣≤C∣ey(φ)∣,
we deduce
[TABLE]
By the Sobolev embedding theorem in R2
[22, Th. 9.9], we have
[TABLE]
The estimate (6.3) results from (6.7), (6.8), and the density of D(V,R3) in K0p(V;R3).
We fix a sequence (φn)n∈N⋆⊂Kp(a,α;S,V) satisfying
∫Vf(ey(φn))dy≤capf(a,α;S,V)+n1∀n∈N⋆.
The set Kp(a,α;S,V) is convex and strongly closed in Kp(V;R3), thus weakly closed.
By (3.5)
and (6.3),
(φn)n∈N⋆ is bounded
in
Kp(V;R3),
thus weakly converges, up to a subsequence,
to some
φ∈Kp(a,α;S,V).
Taking (i) into account, we deduce
that capf(a,α;S,V)≤∫Vf(ey(φ))dy≤liminfn→+∞∫Vf(ey(φn))dy=capf(a,α;S,V),
thus φ is a solution to
(6.1). Its uniqueness when f is strictly convex is straightforward.
Taking (ii) into account, we deduce that capf
satisfies (6.1).
(iv) The estimate (6.4) follows from (6.7) and the Poincaré inequality.
We conclude
by replacing (6.3) by (6.4)
and K0p(V;R3) by W01,p(V;R3) in the previous argument.
(v) Given η∈D(V) such that η=1 in S, the field
φ(y):=η(y)(a+diamS2α∧y) belongs to Wp(a,α;S,V), thus
[TABLE]
yielding (6.5).
Assume that 1<p<2 and
let
V1 be a bounded Lipschitz domain of R2 such that S⊂V1⊂V.
Let φ be a solution to (6.1).
By Hölder’s inequality, ∣φ∣Lp(V1;R3)≤∣V1∣21∣φ∣Lp⋆(V1;R3), hence
∣φ∣W1,p(V1;R3)≤(1+∣V1∣21)∣φ∣Kp(V;R3),
yielding
∣φ∣LH1p(∂S;R3)≤C(S,V)∣φ∣Kp(V;R3) by the continuity of the trace from W1,p(V1;R3) to LH1p(∂S;R3). Since φ⊂Kp(a,α;S,V),
taking (3.5) and (6.3) into account, we deduce
[TABLE]
Noticing that ∣a∣p+∣α∣p≤C(S)∫∂Sa+diamS2α∧ypdH1(y) for some C(S)>0,
(6.6) is proved for 1<p<2.
The proof of the other case is similar.
∎
We establish below some continuity properties for
(a,α)→capf(a,α;S,V).
Lemma 4**.**
The mapping
[TABLE]
is convex (resp. strictly convex if f is strictly convex and either 1<p<2 or V is bounded in one direction), hence continuous. The functional
[TABLE]
is convex (resp. strictly convex under the above stated additional assumptions) and strongly continuous, hence weakly lower-semicontinuous.
Proof. Let us fix λ∈(0,1), t>0, and ((a1,α1),(a2,α2))∈(R3×R3)2
such that (a1,α1)=(a2,α2).
By (1.14), there exists (φ1,φ2)∈Wp(a1,α1;S,V)×Wp(a2,α2;S,V) satisfying
[TABLE]
Observing that λφ1+(1−λ)φ2∈Wp(λa1+(1−λ)a2,λα2+(1−λ)α2;S,V),
we deduce from
the convexity of f that
[TABLE]
hence
the mapping (6.9)
is convex.
If f is strictly convex and 1<p<2, let
φk now denote a solution to
PKf(ak,αk;S,V) (see Lemma 3 (iii)). Then (6.11) and (6.12)
hold with t=0 and, since ey(φ1)=ey(φ2), the second inequality in (6.12) is strict, thus the mapping
(6.9) is strictly convex.
The same holds for the functional (6.10).
We obtain the same conclusion when V is bounded in one direction.
By (6.5),
the convex functional (6.10) is
bounded on the unit ball of (Lp(Ω;R3))2, thus strongly continuous, so weakly lower-semicontinuous. ∎
We state below
some monotonicity properties of capf(a,α;S,V) with respect to V,
S and f, whose proofs
are so easy that we omit them.
Lemma 5**.**
(i) Let V1, V2 be two open subsets of R2 such that S⊂V1⊂V2, and
S1, S2 two nonempty bounded Lipschitz domains of R2 such that S1⊂S2⊂V
and yS1=yS2=0. Then
[TABLE]
(ii) Let f1 and f2 be two convex functions on S3 such that 0≤f1(M)≤f2(M)∀M∈S3. Then,
for all (a,α,λ)∈(R3)2×(0,+∞),
[TABLE]
(iii) We have
[TABLE]
We now address
the asymptotic behavior of capf(a,α;.,.)
w.r.t. monotonous sequences of sets. Recall
that, for any convex function h:S3→R satisfying (3.5), the following holds [32, Prop. 2.32]:
[TABLE]
Lemma 6**.**
(i) Let (Vn)n∈N⋆ be an nondecreasing sequence of open subsets of R2 and (Sn)n∈N⋆ a nonincreasing sequence of bounded Lipschitz domains of R2 such that
[TABLE]
Then,
[TABLE]
Moreover, for all (a,α)∈(Lp(Ω;R3))2,
[TABLE]
(ii) Assume that either p<2 or V is bounded in one direction and let
(Sn′)n∈N⋆ be an nondecreasing sequence of bounded Lipschitz domains such that
[TABLE]
Then,
[TABLE]
Proof.
(i) We fix t>0, ε>0, set St:={x∈R2,dist(x,S)<t} and choose
ηt∈D(V;[0,1]) and φ∈Wp(a,α;S,V)∩D(V;R3)
such that
[TABLE]
By (6.17), there exists n0∈N⋆ such that
Sn⊂St and suppφ⊂Vn, ∀n≥n0. We set φt:=ηt(a+diamS1α∧y−φ)+φ.
Noticing that
[TABLE]
we deduce from (6.16) that ∫Vf(ey(φ))dy−∫Vf(ey(φt))dy≤C∣S2t∖S∣≤Ct. Taking (6.13) and (6.17) into account, we infer
[TABLE]
and deduce (6.18).
If (a,α)∈(Lp(Ω;R3))2,
by (6.5), (6.13) and (6.17) we have
[TABLE]
thus (6.19) results from (6.18) and the dominated convergence theorem. ∎
(ii) Assume that 1<p<2 and let φn be a
solution to
PKf(a,diamSdiamSn′α;Sn′,V).
By (6.13) and (6.20), we have
[TABLE]
thus, by (3.5) and (6.3), (φn)n∈N⋆ is bounded in K0p(V;R3), hence weakly converges in K0p(V;R3),
and a.e. converges in V,
up to a subsequence,
to some φ. Since φn∈Kp(a,diamSdiamSn′α;Sn′,V), we deduce from (6.20) that φ=a+diamS1α∧y a. e. in S, thus φ∈Kp(a,α;S,V) and, by Lemma 3 (i), (6.1) and (6.22),
[TABLE]
The proof of (6.21) when V is bounded in one direction is similar. ∎
The next lemma is crucial
for the proof of the upper bound.
We set:
[TABLE]
Lemma 7**.**
For every (a,α)∈(C∞(Ω;R3))2,
[TABLE]
Proof. The direct proof of this lemma is considerably shortened as follows:
we set
[TABLE]
We check below that
[TABLE]
therefore, by [17, Lemma 4.3]
(a corollary from [21, Th. 1]),
Verification of (6.25). For each i∈I, let ηi∈\textgothWbp(a,α;S,V) be such that gi(x)=∫Vf(ey(ηi(x,y)))dy. Then,
by (6.23), ∑i∈Iφi(x)ηi(x,y))∈\textgothWbp(a,α;S,V),
hence g(x):=∫Vf(ey(∑i∈Iφi(x)ηi(x,y)))dy∈H and,
by the convexity of f, g≤∑Iφigi.
∎
Remark 7**.**
If 0<k<+∞, by the same argument, one can check that
[TABLE]
where ghom is given by (3.15).
An analogous statement holds in the case 1<κ<+∞.
The first equation in (5.17) and the existence of φ0 satisfying (7.65) in the proof of Theorem 2
can be justified in the same manner.
We particularize below the case of a positively homogeneous function.
Lemma 8**.**
If f is positively homogeneous of degree p∈(1,+∞), then for every (a,α)∈(R3)2, the following holds:
[TABLE]
Proof. The first line of (6.26) is straightforward. Fix (t,λ)∈(0,+∞)2, (a,α)∈(R3)2, φ∈Wp(a,α;λS,λV) such that
capf(a,α;λS,λV)+t≥∫λVf(ey(φ))dy and set φ~(y):=φ(λy). Then φ~∈Wp(a,α;S,V), ey(φ~)(y)=λey(φ)(λy), and, by the
p-positive homogeneousness of f and the change of variables formula,
[TABLE]
∎
The main purpose of the following three propositions is to prove that cf is well defined by (1.13) and satisfies
the estimates (1.9) and (1.10) and the formula (1.15) synthetized in
(3.18) (see Corollary 2).
To that aim, we investigate the asymptotic behavior as r→0 of the sequence (capf(a,α;rS,RrD))r>0 when (Rr)r>0 satisfies
[TABLE]
Proposition 6**.**
Under
(3.8) and (6.27), the following holds
for p∈(1,+∞):
[TABLE]
Proof.
Since f∞,p is p-positively homogeneous, by (6.18) and (6.26) we have
[TABLE]
By (3.5) and (3.7),
f≤Cf∞,p
for some C>0, hence
capf≤Ccapf∞,p.
Denoting by φr a solution to
PKf∞,p(a,α;rS,RrD) if
capf(a,α;rS,RrD)≥capf∞,p(a,α;rS,RrD),
to PKf(a,α;rS,RrD) otherwise, we deduce
[TABLE]
and infer from (3.8), (6.27) and Hölder’s inequality
[TABLE]
Observing that, by (6.13) and (6.26), capf∞,p(a,α;rS,RrD)≤r2−pcapf∞,p(a,α;S,V)
as soon as (Rr/r)D⊃V, we deduce from (6.5) and (6.27) that
When p=2,
explicit computations
developed in [15, pp. 73-75] and [8, p. 147]
in the setting of linear isotropic elasticity, show that (see (3.20))
[TABLE]
where Capfλ0,μ0(rD,RrD)∈S4 is diagonal and satisfies, as r→0,
[TABLE]
Hence, by (6.13) and (6.14), if f is an arbitrary convex function
satisfying (3.5),
[TABLE]
for some c,C∈(0,+∞). In what follows, we denote by T the Fréchet topology on C(R3) of the uniform
convergence on the compact subsets of R3.
We set
[TABLE]
Proposition 7**.**
Assume that p=2. There exists a convex function c0f
independent of S, verifing
[TABLE]
and such that, for every sequence (Rr)r>0 satisfying (6.27),
[TABLE]
The demonstration of this result, postponed to Section 8.2,
is delicate: the convergence
(6.36)
is obtained
by capitalizing on its consequence, Theorem 1.
Proof. Let φ be
a solution to Pf(a,diamSdiamDζe3;rD,RrD).
By (3.1), (3.5), (6.7), and (6.13), we have
[TABLE]
Applying to each component of φ the estimate (see [12, Lemma A3])
[TABLE]
noticing that, by (1.14),
∫−∂RrDφdH1=0, ∫−∂rDφdH1=a,
we deduce
[TABLE]
Since p>2, by Hölder’s inequality we have φ∈W2(a,diamSdiamDζe3;rD,RrD),
hence, by (3.5), (6.13), (6.33), (6.39), and the quadratic nature of the mapping (a,ζ)→cap∣.∣2(a,ζ;rD,RrD) (see (6.31)), the following holds:
[TABLE]
We infer that capf(a,ζe3;rS,RrD)≥Rrp−2Cζ2−∣logr∣C2p which, combined with
(6.40), yields (6.38). ∎
Remark 9**.**
One can check that, if p=2 and 0<R1<R2,
[TABLE]
where s:=p−1p−2.
The field defined in polar coordinates by
φ(ρ,θ):=η(ρ)a, where η is the solution to the above problem,
belongs to Wp(a,0;R1D,R2D), therefore
cap∣.∣p(a,0;R1D,R2D)≤C∫R2D∣∇φ∣pdx≤(R2s−R1s)p−1C∣a∣p.
Setting R1=1, letting R2 tend to +∞, we deduce from (3.5), (6.14) and (6.18) that cap∣.∣p(a,0;D,R2)=capf(a,0;S,R2)=0.
The following corollary results from
(6.5), (6.6) and
Proposition 6 if p<2, from (6.33) and Proposition 7 if p=2, and from Proposition 8 if p>2.
Corollary 2**.**
For every sequence (Rε) verifying (3.17), the convergence (1.13) holds, where
cf is convex,
independent of (Rε), satisfies (3.18). If p=2, cf is independent of S and
Following the Γ-convergence method [33], we establish
a lower bound in Section 7.1.1, an upper bound
in Section 7.1.2,
and conclude
the proof of Theorem 1
in Section 7.1.3.
7.1.1. Lower bound
Proposition 9**.**
Let (uε) be a sequence in Wb1,p(Ω;R3) verifying (5.1) and
(uεk)k∈N a subsequence satisfying the convergences
stated in Propositions 3, 4 and 5.
Then
[TABLE]
Proof. To simplify, the subsequence will still be denoted by (uε).
By Lemma 10, there exists a sequence (Rε)
verifying (3.17) and
[TABLE]
We set
[TABLE]
Since ∣DRε∣≪1, by (5.59) the sequence (e(uε)\mathds1Ω∖(DRε×(0,L))) weakly converges to e(u) in Lp(Ω;S3). Therefore, by
the weak lower-semicontinuity on Lp(Ω;S3) of the functional
Ψ→∫Ωf(Ψ)dx, resulting from (3.5) and the convexity of f, we have
[TABLE]
Let us check that
[TABLE]
If (k,κ)∈{(+∞,0),(+∞,+∞)}, then by (3.15), (3.16), (3.19) and (5.92), we have
\textgothDvtuple=0, thus there is nothing to prove.
If 0<k<+∞, by
(5.74) we have e(uε)⇀⇀\buildrelmε\textgothL(q,∂x3∂v3,∂x3∂θ) for some q∈Lp(Ω;W1,p(S;R3)). Applying (5.15), taking (1.2) and (3.15) into account, we obtain
[TABLE]
yielding (7.5).
If 0<κ<+∞, by (1.2),
(3.8) and (5.63),
for some l∈Lp(Ω;W1,p(S;R3)).
Taking (3.9) into account, (7.5) is proved.
By (3.14), (5.92), (7.3), (7.4) and (7.5), the demonstration of Proposition 9 is achieved
provided we show that
[TABLE]
If γ(p)=0, by (3.18) we have ∫Ωcf(v−u,θ)dx=0, hence there is nothing to prove.
Likewise, if γ(p)=+∞,
by (5.61) there holds u=v and θ=0, hence, by (3.18), ∫Ωcf(v−u,θ)dx=0.
From now on, we assume that 0<γ(p)<+∞, thus p∈(1,2].
We fix c0∈(0,1) and a Lipschitz domain S′ such that
[TABLE]
By Lemmas 9 and 11
there exists
sequences
(\widetriangleuε), (\wideparenuε) in Lp(0,L;W1,p(Ω′;R3))
satisfying
[TABLE]
[TABLE]
[TABLE]
We set
[TABLE]
The following weak convergences in Lp(Ω) hold:
[TABLE]
We respectively denote by
aεi(x3), αεi(x3), and \textswabuε(\widetriangleuε)i(x3) the constant values taken by aε, αε, and \textswabuε(\widetriangleuε) in Yεi×{x3}.
For each i∈Iε and a. e. x3∈(0,L), \wideparenuε(.,x3) belongs to W1,p(DRεi;R3),
ex′(\wideparenuε(.,x3))=ex′(\wideparenuε−\textswabuε(\widetriangleuε)i)(.,x3) in DRεi
and,
by (1.14), (5.4), (7.10), (7.11), (7.12),
Substituting Sn′ for S′,
where
(Sn′)n∈N is an nondecreasing sequence of Lipschitz domains such that
⋃n∈N↑Sn′=S, Sn′⊂⊂S, ySn′=0, we infer
[TABLE]
By (3.18), (6.13), (6.21), the sequence (cn)n∈N is positive, nondecreasing, and pointwise converges to γ(p)capf∞,p(v−u,θe3;S,R2)=cf(v−u,θ).
Applying the
monotone convergence theorem, we obtain (7.7).
∙Case p=2, 0<γ(2)<+∞
By (1.3), (6.34) and (7.14), we have
[TABLE]
Summing w.r.t. i over Iε and integrating w.r.t. x3 over (0,L),
taking (7.8) and (8.27) into account, we infer
[TABLE]
After possibly extracting a subsequence, by (1.3) we can assume that
[TABLE]
For each N∈N, we set
[TABLE]
By (6.26) and (6.34), crε,Rε/c0f∞,2,S(tξ)=t2crε,Rε/c0f∞,2,S(ξ), hence
[TABLE]
Since ∣aεN∣L2(Ω;R3)≤∣aε∣L2(Ω;R3)≤C,
by the metrizability of the weak topology on bounded subsets of L2(Ω;R3),
there exists a subsequence of (aε)ε>0 (denoted the same) such that, for every N∈N,
(aεN)ε>0 weakly converges in L2(Ω;R3) to some aN.
As (aN)N∈N is bounded in L2(Ω;R3), a subsequence still denoted (aN)N∈N
weakly converges as N→+∞ to some a~.
We fix N∈N and α>0. By Proposition 7, (crε,Rε/c0f∞,2,S)ε>0 uniformly converges
on BR3(0,N) to c0f, thus
[TABLE]
We infer from the weak lower semicontinuity on L2(Ω;R3) of
ψ→∫Ωc0f(ψ)dx
[TABLE]
and then from (7.11), (7.16), (7.17), (7.19)
and the arbitrariness of α,
[TABLE]
By passing to the limit inferior as N→+∞, we obtain
By passing to the limit as ε→0, taking (7.13) into account, we infer
[TABLE]
The sequence (aε)ε>0 is bounded in L2(Ω;R3),
hence equiintegrable,
therefore (see [2, Prop. 1.27])
limN→+∞supε>0∫∣aε∣>N∣aε∣dx=0.
By passing to the limit as N→+∞, we infer
∫Ω(v−u−a~)⋅φdx=0
and
deduce
a~=v−u. By (7.20)
and the arbitrariness of S′, (7.7) holds.
Proposition 9 is proved. ∎
and let α>0.
There exists a sequence (φε) in Wb1,p(Ω;R3) such that
[TABLE]
and
[TABLE]
Remark 11**.**
Proposition 10 implies, by density and diagonalisation arguments, that for every u∈Wb1,p(Ω;R3), there exists a sequence (uε) weakly converging to u in Wb1,p(Ω;R3)
such that limsupε→0Fε(uε)≤F(u) (see (3.13)). By
Proposition 9, for every u∈Wb1,p(Ω;R3)
and every sequence (uε) weakly converging in Wb1,p(Ω;R3) to u,
liminfε→0Fε(uε)≥F(u).
Thus, since (Fε)ε>0 is equicoercive on Wb1,p(Ω;R3),
it Γ-converges to F in the weak topology of Wb1,p(Ω;R3) (see [33]).
Proof. Let us fix R>0, a sequence (Rε) satisfying (3.17), and a Lipschitz domain S′ such that
[TABLE]
The approximating sequence will be defined by
[TABLE]
where ξε is such that
[TABLE]
χε is defined either by (7.35), (7.37), or (7.40), depending on the magnitude of k and κ, and ηε is given by (7.42)
if γ(p)<+∞ and p<2,
by (7.47) if γ(p)<+∞ and p=2, and otherwise ηε=0.
These choices ensure that φε belongs to Wb1,p(Ω;R3) and satisfies (7.23).
The support of the function ηε is included in DRε×(0,L), where
It follows from (3.9), (5.11), (7.39) and the abobe inequality, that
[TABLE]
which, combined with (7.36), yields (7.31).
In the remaining cases, we set
[TABLE]
If k=+∞ and κ=0, by (3.19) and (7.21),
ζ=ψ3=0, ψtuple=0, hence ∣e(χε)∣≤Crε in Trε and, by (1.2),
Iε3≤Crεpkε→0=∫Ωghom(\textgothD(ψtuple))dx.
The case k=κ=+∞ is straightforward.
The assertion (7.31) is proved. ∎
Proofs of (7.32) and (7.33).
If γ(p)=+∞, these assertions straightforwardly results from ηε=0. Otherwise,
we distinguish two cases:
∙Case p<2 and γ(p)<+∞.
By (3.18) and (6.19), we can assume that R and S′ satisfy, besides (7.25),
The assertion (7.32) results from the estimates
\Big{|}\boldsymbol{e}\Big{(}\left(1-\xi_{\varepsilon}\right)\left({\boldsymbol{\psi}}+\frac{2}{{\rm diam}S}\zeta\boldsymbol{e}_{3}\wedge\frac{\boldsymbol{y}_{\varepsilon}}{r_{\varepsilon}}\Big{)}+\xi_{\varepsilon}{\boldsymbol{\chi}}_{\varepsilon}\right)\Big{|}≤C and ∣Srε′∖Srε×(0,L)∣≤Cε2rε2,
deduced from (7.26), (7.35), (7.37) and (7.40).
Applying (6.16) with
[TABLE]
integrating over DRrε∖Srε′×(0,L), noticing that ∣Mε−Mε′∣≤C,
∣Mε∣+∣Mε′∣≤rεC, and
∣DRrε∖Srε′×(0,L)∣≤Cε2rε2, we infer
By (7.41), ey(η)(x,rεyε(x′))=0 in Srε′×(0,L),
hence, by (7.44),
[TABLE]
where νε:=rε2∣RD∣ε2\mathds1(RD)rε×(0,L)L⌊Ω3. Notice that νε is deduced from mε by substituting RD for S in (5.9).
By applying
(5.11) to (νε,f∞,p(ey(η))) in place of (mε,ψ), taking the above estimates into account,
we obtain
[TABLE]
The assertion (7.33) is proved.
The proof of Proposition 10 in the case p<2, γ(p)<+∞ is achieved. ∎
∙Case p=2, γ(2)<+∞.
By Lemma 3, we can fix for each a∈R3, ε>0 a field ηεa∈W2(a,0;rεS′,RεD) such that
[TABLE]
We extend each ηεa by [math] to R2. We
split (0,L) into a suitable family of intervals (Jεk)k∈{1,..,nε} and set
where ηεϕεik is given by (7.45). The choice of the function ρε (see below)
ensures that ηε approximates to η~ε and
belongs to H1(Ω;R3).
By (7.45) and (7.46), we have
[TABLE]
The intervals Jεk will be defined by fixing
sequences (aε) and (bε) such that
[TABLE]
and setting ([s]: integer part of s)
[TABLE]
The function ρε introduced in (7.47) is chosen such that
The estimates
e((1−ξε)(ρε(x3)ϕε(x)+φ)+ξεχε)2≤C(1+∣ρ′(x3)∣2), deduced from (7.26), (7.35), (7.37) and (7.40),
and ∫0L∣ρε′∣2dx3≤Caεbε1,
infered from (7.50), yield Iε4≤Cε2rε2(1+aεbε1),
which by (7.49) implies (7.32).
By (1.3), (6.34) and (7.45) we have, for every (i,a)∈Iε×R3,
The solution uε to (1.1)
satisfies
Fε(uε)−∫Ωf.uεdx≤Fε(0)=0, thus, since Wb1,p(Ω;R3)∩R={0},
by (5.30) the following holds
[TABLE]
We deduce
(5.1), and then the existence of a subsequence
(uεk)k∈N satisfying the convergences stated in Propositions 3, 4 and 5 for some (u,vtuple)∈Wb1,p(Ω;R3)×D.
Applying
Proposition 9, we infer
[TABLE]
Let us fix α>0.
The functional Φ is convex and bounded on the unit ball of the Banach space Wb1,p(Ω;R3)×D, hence strongly continuous.
By density, there exists
a couple (φ,ψtuple) satisfying (7.21), (7.22), (7.23)
and
such that
[TABLE]
Applying Proposition 10, we fix (φε)⊂Wb1,p(Ω;R3) verifying (7.23) and (7.24).
Since uε is the solution to (1.1), we have
[TABLE]
By (7.60), (7.61), and the arbitrariness of α,
(u,vtuple) is a solution to (3.12), hence u is a solution to (3.13).
By the strict convexity of F, this solution is unique, hence the entire sequence (uε) weakly converges to u in W1,p(Ω;R3). Theorem 1 is proved.
∎
The proof follows the same pattern as that of Theorem 1.
We establish
the analog of Proposition 9,
assuming (3.25) and
[TABLE]
By (3.24), (3.25) and (7.62), whatever the operator \textswabpε∈{\textswabvε,...}
introduced in Section 5.1,
the sequence (\textswabpε(uε)) is bounded in Lp and
satisfies (5.12).
Repeating their proofs, we find that every assertion stated in Propositions 3, 4 and
5 holds, except the first line of (5.59) where the convergence of (uε) only holds in the weak topology of Lp(Ω;R3)
and that of (\textswabuε(uε)) is irrelevant.
Setting
[TABLE]
we infer (7.5) (same proof).
We deduce from
[13, Lemma 3.2] that, up to a subsequence, the following holds for some u0∈Lp(Ω;W♯1,p(Y;R3)):
[TABLE]
where "⇀⇀" denotes the ”usual” two-scale convergence, defined by substituting
L3 for mε in (5.10). In other words, uε⇀⇀u0 means that for all η∈C(Ω×Y;R3),
[TABLE]
On the other hand, by (5.10) and the third line of (5.59),
[TABLE]
By varying η in D(Ω×S;R3), taking (3.24) and (5.9) into account, we infer that u0(x,y)=v(x)+diamS2θ(x)e3∧y in Ω×S,
yielding, by (1.17) and (7.63),
[TABLE]
By (7.63) we have εe(uε)\mathds1Ω∖Tε⇀⇀ey(u0)\mathds1Y∖S
(see [7, Lemma 1]).
In view of (3.8) and
(5.15) applied to ”⇀⇀”, we infer
liminfε→0Iε1≥∫Ω×Y∖Sf∞,p(ey(u0))dxdy and deduce from (3.26) and (7.5) and (7.64) that
[TABLE]
To prove the couterpart of Proposition 10, we fix α>0 and (φ,ψ,ζ,a,b) satisfying (7.21), (7.22). Mimicking the proof of Lemma 7, we obtain (see (1.17))
[TABLE]
hence there exists φ0 such that φ0−φ∈\textgothW(ψ−φ,ζ) and
[TABLE]
Setting Tεε:={x∈Ω,dist(x,Tε)<ε2}, we fix
(ξε)⊂C∞(Ω) satisfying
[TABLE]
and put
[TABLE]
where χε is defined
by (7.35), (7.37) or (7.40)
(depending on k, κ).
Writing
[TABLE]
noticing that
e(φ0(x,rεyε(x′)))=ε1ey(φ0)(x,rεyε(x′))+O(1) and
∣e(φε)∣≤εC in Tεε∖Tε,
taking (3.8) and (5.11) applied to "⇀⇀"
into account,
we deduce that
thus, by (1.2), (5.41), (5.42) and the first inequality in (5.63), the following holds
[TABLE]
Combining (7.66) and (7.68), we infer (5.1), hence the convergences stated in Propositions 3, 4, 5 hold
for some subsequence (uεk). We deduce from (4.15), (5.59), (5.60), (5.61), (5.91), and (7.67), that
By density, there exists
(φ,ψtuple) satisfying (7.21), (7.22)
and
[TABLE]
and, by Proposition 10 a sequence (φε) verifying (7.23) and (7.24). By (4.15), (7.23) and (7.67),
limε→0∫Ωφε⋅fεdx=∫Ωφ⋅fdx+L(ψtuple), hence, by (7.24) ,
[TABLE]
This, along with (7.69), shows, by the arbitrariness of α, that (u,vtuple) is a solution to Problem 4.12.
The case κ=0 can be proved in a similar manner, by using the weak⋆ convergence in M(Ω),
stated in (5.59), of (−rεyε2uε1+rεyε1uε2)mε to τθ.
∎
Problem 4.10.
Noticing that the solution uε to (4.10) satisfies
(7.66), (7.68) with Fεsoft in place of Fε (same proof), hence, by
(3.25), is bounded in Lp(Ω;R3), we
adapt the proof of Theorem 2 in the same way as above.∎
8. Appendix
8.1. Technical lemmas related to the lower bound
In this section, given a sequence (uε) satisfying (5.1),
we establish the existence of
sequences (\widetriangleuε) and (\wideparenuε) verifying
(7.9), (7.10) and
(7.11).
Lemma 9**.**
Assume that p≤2, (Rε) verifies (3.17), and
(uε)
satisfies (5.1) and (5.59).
Let (Iε2) be defined by (7.3). There exists a sequence (\widetriangleuε) in Lp(0,L;W1,p(Ω′;R3)) such that
[TABLE]
Proof.
To construct \widetriangleuε, we fix two sequences (aε) and (bε) such that
[TABLE]
and choose
a family (lk,ε)ε>0,k∈{1,...,nε} verifying
[TABLE]
and such that, setting
[TABLE]
the following holds:
[TABLE]
To prove that this choice is possible,
we set, for every
m∈{0,...,[2bεaε]}
[TABLE]
Since Hεm∩Hεn=∅ if m=n, we have
[TABLE]
Accordingly, for each ε>0, there exists mˇε∈{0,...,[2bεaε]} such that
[TABLE]
Therefore
(lk,ε)k∈{1,...,nε} defined by lk,ε:=lk,εmˇε
satisfies (8.3) and
(8.5).
Let us fix (ρε) verifying (7.50).
To each sequence (ψε) in Lp(Ω;R3), we associate the sequence (\widetriangleψε) defined by
[TABLE]
Since 0≤ρε≤1, for every Borel subset E of Ω′, we have
Since ∫∣uε∣pdmε≤C (see (5.64)),
by (8.10) and Lemma 1
the following holds, up to a subsequence, for some (g,h,η)∈(Lp(Ω;R3))2×Lp(Ω):
[TABLE]
Taking (5.21) into account, the convergences stated in (8.1) are proved provided we establish that
(g,h,η)=(u,v,θ).
To that aim, we fix ξ∈D(Ω;R3) and set
ξε(x):=∑k=2nε(∫−Jεkξ(x′,s3)ds3)\mathds1Jεk(x3).
By (5.4) and (8.7), we have
[TABLE]
By (5.21) and (5.59), (\textswabvεS(uε)mε) weakly⋆ converges to v in M(Ω;R3) and,
by (5.65), (7.50), (8.2), (8.4), and (8.5),
[TABLE]
therefore (ρε\textswabvεS(uε)mε) weakly⋆ converges to v in M(Ω;R3).
Observing that ∣ξ−ξε∣L∞(Ω;R3)≤Caε≪1,
we deduce from (8.11) and (8.12) that
[TABLE]
and infer from the arbitrariness of ξ that h=v. We likewise obtain η=θ and g=u. The convergences stated in (8.1) are proved.
Let Eε be a Borel subset of DRε.
By (3.8), (5.63), and the estimate ∣DRε∣≤Cε2Rε2, we have
[TABLE]
By (7.50) and (8.4), ρε(x3)uε(x)=uε(x) in (Eε×(0,L))∖Hε and ∣ρε′∣<bεC,
hence by
(3.5), (5.63), (8.2), and (8.5),
[TABLE]
It follows from (5.63), Hölder’s inequality, (8.2), (8.4), and the continuous embedding of W1,p(Ω;R3) into L3−p3p(Ω;R3)
that
[TABLE]
Since ρε=0 on ⋃k=1nε{lk,ε} (see (7.50)), we have ∫−Jεk∂s3∂(ρε(s3)uε(x′,s3))ds3=0 a. e. x′∈Ω′,
∀k∈{2,...,nε}.
Hence, for a. e. x∈Ω′×Jεk, there holds (see (2.1), (8.7))
[TABLE]
Applying Jensen’s inequality, we infer
[TABLE]
which, combined with (8.13), (8.14), and (8.15), yields
[TABLE]
Choosing first Eε=DRε∖Srε, taking (7.3) into account,
we deduce
the lower bound stated in the second line of (8.1).
Next, substituting ∣.∣p for f in the above argument, we obtain
the third line of (8.1). ∎∎
We check below the existence of a sequence (Rε)
verifying (3.17) and (7.2).
Lemma 10**.**
Assume that
(uε) verifies (5.1).
For any sequence (Rε′) satisfying (3.17), there exists (Rε)
verifying (3.17) and (7.2) such that 0<Rε≤Rε′∀ε>0.
Proof. Let (Rε′) be satisfying (3.17).
There exists a sequence of positive integers (nε) such that limε→0nε=+∞ and (2−nεRε′) verifies (3.17).
By (5.1), we have
[TABLE]
hence there exists mε∈{0,..,nε−1} such that
[TABLE]
We set Rε:=2−mεRε′.
∎
Lemma 11**.**
Assume that p≤2,
(uε)
satisfies (5.1), (Rε) verifies (3.17) and (7.2), and
the convergences (5.59) hold. Let
(\widetriangleuε) given by Lemma 9 and
S′ a Lipschitz domain verifying (7.8).
Then, there exists (\wideparenuε)⊂Lp(0,L;W1,p(Ω′;R3)) satisfying
(7.10) and (7.11).
Proof of (8.18).
By (5.63) and the last line of (8.1), the following holds
[TABLE]
We deduce
[TABLE]
By (2.1) and (5.8) we have ex′(\textswabuε(\widetriangleuε))=0 in DRε×(0,L) and, by (8.9), (8.16) and (8.17),
\widetriangleuε−\wideparenuε=ψε(\widetriangleuε−\textswabuε(\widetriangleuε))=ψε(\widetriangleuε−\textswabuε(uε)) in
(DRε∖Srε)×(0,L).
Applying the last line of (8.1) to Eε:=DRε∖DRε/2, taking (7.2), (8.8),
(8.16) and (8.21) into account, we infer
[TABLE]
Let E be a bounded Lipschitz domain of R2.
One can check that
\int_{E}\big{|}\psi-{\int\hskip-8.5359pt-}_{E}\psi d{\mathcal{L}}^{2}\big{|}^{p}d{\mathcal{L}}^{2}\leq C\int_{E}\left|\nabla\psi\right|^{p}d{\mathcal{L}}^{2}∀ψ∈W1,p(E).
By making suitable changes of variables,
we infer
[TABLE]
Applying this to (E,ψ(.))=(D∖21D,uε(.,x3)), summing w.r.t. i over Iε and integrating w.r.t. x3 over (0,L),
taking (5.8) and (7.2) into account,
we deduce
[TABLE]
Combining this with (8.22) and (8.23) the proof of (8.18) is achieved. ∎Proof of (8.19). If 1<p<2, by (3.5), (8.16) and (8.17), we have
[TABLE]
Noticing that, by (5.4), ex′(\textswabrεS(\widetriangleuε))=0 in Trε,
we deduce from (5.9), (5.37), (5.63), and the last line of (8.1) applied to Eε=Srε, that
Proof of (8.20). If p=2, setting Grε:=(Srε∖Srε′)×(0,L), we infer from (5.4), (8.16), (8.17), and (8.24) (which also holds for p=2) that
[TABLE]
Noticing
that ex′(θεS(\widetriangleuε)e3∧rεyε(x′))=0 in Grε, θεS(\widetriangleuε) takes constant values in each set Srεi×{x3},
and ∣Trε∣∣Grε∣=∣S∣∣S∖S′∣, taking (5.9) and (8.16) into account, we deduce
[TABLE]
By
(5.34), (5.38),
(5.63), (8.8), and (8.9), we have, recalling that θε(uε)=wε3(uε)
and cε=ε2 if p=2,
[TABLE]
Combining the above estimates, we obtain (8.20). ∎
hence c0f=c0f∞,2 and it suffices to prove Proposition 7 when f is 2-positively homogeneous, which is assumed in the sequel.
We first establish the T-relative compactness of (cr,Rrf,S)r>0.
Lemma 12**.**
Under (6.27), the sequence (cr,Rrf,S)r>0 defined by (6.34) is uniformly equicontinuous on the compact subsets of R3,
hence T-relatively compact.
Its cluster points are
convex and satisfy (6.35).
Proof. By (6.33) and Lemma 4, any cluster point of (cr,Rrf,S) is convex and satisfies (6.35).
We prove below that
[TABLE]
hence (cr,Rrf,S)r>0 is uniformly equicontinuous on the compact subsets of R3, thus, by
Ascoli’s Theorem and Cantor’s diagonal process,
T-relatively compact. We turn to the proof of (8.26):
by
(6.13), (6.26) and (6.34), for all r,R,R′,t>0 and
all bounded Lipschitz domains S1, S2 of R2 such that R≤R′, S1⊂S2⊂RD,
yS1=yS2=0,
tS⊂RD, and S⊂RD, we have
[TABLE]
Let φr be the radial function defined on RrD in polar coordinates by
φr(ρ,θ):=ψr(ρ), where ψr is the solution to
[TABLE]
One can check that
[TABLE]
Let ψa,r be a solution to
Pf(a,0;rS,RrD).
Then ψa,r+hφr∈W2(a+h,0,rS,RrD) and, by (3.5), (6.16), (6.27), (6.33), and (8.28),
if ∣h∣≤C,
[TABLE]
Applying the same argument to (a′,h′):=(a+h,−h), we deduce (8.26).
∎
By Lemma 7, Proposition 7 is proved if we show that
(cr,Rrf,S) has a unique T-cluster point as r→0, independent of (Rr)r>0 and S.
We achieve this
task by making use of
the proof of Theorem 1. Setting
[TABLE]
we consider the functional Fε given by
(1.1). By (1.2), (1.3) and (8.29), we have
[TABLE]
Let If be the weak Γ-lower limit of (Fε)
in Wb1,2(Ω;R3), that is the functional defined by (see [33])
[TABLE]
We set
[TABLE]
We will prove that every cluster point of (cr,Rrf,S) as r→0 belongs to C. Hence, by Lemma 12, C=∅. The next implication, easily proved by contradiction,
[TABLE]
shows that C has a unique element
c0f. It follows that the whole sequence (cr,Rrf,S)T-converges
to c0f as r→0, for every sequence (Rr)r>0 verifying (6.27).
To prove that c0f is
independent of S, let us fix λ>1 such that D⊂S⊂λD.
By what precedes,
(cr,Rr/λf,D) and (cr,Rrf,D)
both T-converge to some c0f,D∈C(R3).
By passing to the limit as r→0 in the inequalities
cr,Rrf,D≤cr,Rrf,S≤cr,Rrf,λD=cr,Rr/λf,D, deduced from (8.27), we infer c0f=c0f,D,
thus c0f is
independent of S.
The proof of Proposition 7 is achieved provided we establish the following lemma:
Lemma 13**.**
Every T-cluster point of (cr,Rrf,S) as r→0
belongs to C.
Proof.
Let c be such a cluster point and (rk)k∈N a decreasing sequence converging to [math] such that r1<1 and
We prove below that for all φ∈D(Ω;R3), there exists
(φε)ε>0∈Cφ such that
[TABLE]
Since liminfε→0Fε(φε)≤limsupk→+∞Fεk(φεk), we deduce from (8.31) that
[TABLE]
Next, we establish that
[TABLE]
We infer from (8.37), (8.38), and
the arbitrariness of φ, u that c∈C.
It remains to exhibit (φε)∈Cφ satisfying (8.36)
and to prove(8.38).
∙ Given φ∈D(Ω;R3), the sequence (φε)∈Cφ verifying (8.36)
will be deduced from (7.26) (upperbound)
by specifying the choice of
Rε in (7.45), namely by setting
By (7.46), (φε) is bounded in L∞ and uniformly converges
on each compact subset of Ω to φ,
therefore, by (8.33),
[TABLE]
Taking (7.29) and (8.44) into account, the assertion (8.36) is proved.
Verification of (8.42).
By (6.27), rk≪Rrk, hence, by (8.40),
rε≪Rε(1).
By (8.29),
the mapping ε→rεε is decreasing on (0,2) hence
Rε(1)=rkrεRrk≤εkεRrk for all ε∈]εk+1,εk]
if εk<2.
Since
Rrk≪∣logrk∣1=εk,
we deduce that Rε(1)≪ε=∣logrε∣1 (see (8.29)), thus
(ε,rε,Rε(1)) verifies (3.17).
By (8.27), (8.35), (8.40), and the
decrease of ε∈(0,+∞)→∣logrε∣,
∙Proof of (8.38). We fix u∈D(Ω;R3),
(uε)ε>0∈Cu, and a subsequence (uεl)l∈N such that
[TABLE]
By (8.31), (uεl) verifies (5.1),
hence, up to a subsequence, every assertions stated in Propositions 3,4,5.
Let us fix S′ and c0 as in (7.8).
By (8.42) and Lemma 10, there exists (Rε) such that
up to a further subsequence,
for some c∈C(R3) satisfying, by (8.33) and (8.48), c≥c.
Using (8.49) in place of (6.36) in the argument of the lower bound,
noticing that by (3.19), (5.92) and (8.30), vtuple=0, we obtain
The reference list from the paper itself. Each links out to its DOI / PubMed record.
1[1] Allaire, G.: Homogenization and two-scale convergence. SIAM J. Math. Anal. 23, 1482–1518 (1992)
2[2] Ambrosio, L., Fusco, N., Pallara, D.: Functions of bounded variation and free discontinuity problems. Oxford Mathematical Monographs, Oxford University Press, New York, (2000)
3[3] Arbogast, T., Douglas, J., Hornung, U.: Derivation of the double porosity model of single phase flow via homogenization theory, SIAM J. Math. Anal., 21, 823–836 (1990)
4[4] Avila, A., Griso, G., Miara, B., Rohan, E.: Multiscale modeling of elastic waves: theoretical justification and numerical simulation of band gaps. Multiscale Model. Simul., 7, 1–21 (2008)
5[5] N.O. Babych, I.V. Kamotski, V.P. Smyshlyaev: Homogenization in periodic media with doubly high contrasts, Networks and heterogeneous Media, 3 (3), 413–436 (2008)
6[6] Bellieud, M.: Homogenization of evolution problems in a fiber reinforced structure, J. Convex Anal. Vol. 11, No 2, 363–385 (2004)
7[7] Bellieud, M.: Torsion effects in elastic composites with high contrast. SIAM J. Math. Anal. 41(6), 2514–2553 (2010)
8[8] Bellieud, M.: A notion of capacity related to linear elasticity. Applications to homogenization. Arch. Ration. Mech. Anal., vol. 203, no 1, 137–187 (2012)