François Murat
email: murat@ann.jussieu.fr

Ali Sili
email: ali.sili@univ-amu.fr

A REMARK ABOUT THE PERIODIC HOMOGENIZATION OF CERTAIN COMPOSITE FIBERED MEDIA

Keywords: 2010 Mathematics Subject Classification. Primary: 35B25, 35B27, 35B40; Secondary: 76M50, 74K10 Composite media, fibers, strange term, nonlocal problem, homogenization, dimension reduction

We explain in this paper the similarity arising in the homogenization process of some composite fibered media with the problem of the reduction of dimension 3d -1d. More precisely, we highlight the fact that when the homogenization process leads to a local reduction of dimension, studying the homogenization problem in the reference configuration domain of the composite amounts to the study of the corresponding reduction of dimension in the reference cell. We give two examples in the framework of the thermal conduction problem: the first one concerns the reduction of dimension in a thin parallelepiped of size ε containing another thinner parallelepiped of size rε ≪ ε playing a role of a «hole ». As in the homogenization, the one-dimensional limit problem involves a "strange" term. In addition both limit problems are of the same structure. In the second example, the geometry is similar but now we assume a high contrast between the order 1 of the conductivity in the small parallelepiped of size rε := rε, (r > 0 fixed ) and the conductivity in the big parallelepiped assumed to be of order ε 2 . We prove that the limit problem is a nonlocal problem and that it has the same structure as the corresponding periodic homogenized problem.

1. Introduction. During the last years, the study of the homogenization of composite heterogeneous media has given rise to an extensive literature and a significant part of that works was designed to the homogenization of media characterized by high heterogeneities, (see [START_REF] Arbogast | Derivation of the double porosity model of single phase flow via homogenization theory[END_REF], [START_REF] Bellieud | Homogenization of elliptic problems in a fiber reinforced structure. Nonlocal effects[END_REF]), [START_REF] Boughammoura | Homogenization and correctors for composite media with coated and highly anisotropic fibers[END_REF], [START_REF] Braides | Homogenization of non-uniformly bounded periodic diffusion energies in dimension two[END_REF], [START_REF] Braides | A variational approach to double-porosity problems[END_REF], [START_REF] Brillard | Asymptotic behaviour of a cylindrical elastic structure periodically reinforced along identical fibers[END_REF], [START_REF] Rauch | Potential and scattering theory on wildly perturbed domains[END_REF]). The pioneer work for problems of this kind was done in [START_REF] Arbogast | Derivation of the double porosity model of single phase flow via homogenization theory[END_REF] in the study of the double porosity model of single phase flow. Thereafter, the main idea of [START_REF] Arbogast | Derivation of the double porosity model of single phase flow via homogenization theory[END_REF] was taken up in [START_REF] Allaire | Homogenization and two-scale convergence[END_REF] and [START_REF] Nguetseng | A general convergence result for a functional related to the theory of homogenization[END_REF] to give rise to the two-scale convergence method. Fibered media is an example of composite heterogeneous media with high contrasting properties since usually the material constituting the fibers is very different from the material around it. For instance, in elasticity one can consider rigid fibers immersed in a soft matrix while in the framework of heat conduction one can consider fibers with high conductivity surrounded by a material with a poor conductivity. Among the first works devoted to the homogenization of such composite media one can quote [START_REF] Caillerie | A perturbation problem with two small parameters in the framework of the heat conduction of a fiber reinforced body[END_REF] where the homogenization process was performed using asymptotic expansions.

From the mathematical point of view, the contrast between the properties of the two materials leads to a degenerate problem in the sense that in general it leads to a lack of compactness. Indeed, the operators under consideration are in general not uniformly bounded (see [START_REF] Braides | Homogenization of non-uniformly bounded periodic diffusion energies in dimension two[END_REF]) or not uniformly coercive with respect to the small parameter (see [START_REF] Sili | Homogenization of a nonlinear monotone problem in an anisotropic medium[END_REF]).

In general, the configuration domain of such media may be described by a domain Ω of R 3 which is a periodic replication with a period of size ε of a set of Y ε i , i ∈ I ε . More precisely, Y ε i is assumed to be the union of a fiber F ε i with its complement M ε i ( M stands for matrix ) in Y ε i . It is well known that the homogenization process in such degenerate problems gives rise to homogenized problems with a different form from the equation at the microscopic level since memory effects, strange terms or nonlocal effects may appear at the limit (see [START_REF] Bellieud | Homogenization of elliptic problems in a fiber reinforced structure. Nonlocal effects[END_REF], [START_REF] Braides | Homogenization of non-uniformly bounded periodic diffusion energies in dimension two[END_REF], [START_REF] Braides | A variational approach to double-porosity problems[END_REF], [START_REF] Brillard | Asymptotic behaviour of a cylindrical elastic structure periodically reinforced along identical fibers[END_REF], [START_REF] Casado-Diaz | Two-scale convergence for nonlinear Dirichlet problems[END_REF], [START_REF] Cioranescu | Un terme étrange venu d'ailleurs[END_REF], [START_REF] Sili | Homogenization of a nonlinear monotone problem in an anisotropic medium[END_REF]).

The aim of the present work is to show that in the case of fibered media such effects at the limit are not due to the homogenization process itself but to the local structure of the composite media; more precisely, we show that the form of the homogenized problem is already determined by the study of the 3d -1d reduction of dimension which occurs locally. To illustrate that, we give here two examples in the framework of the thermal conduction. In the first one we show that one can obtain an extra term (or a strange term, see [START_REF] Cioranescu | Un terme étrange venu d'ailleurs[END_REF], [START_REF] Bellieud | Homogenization of elliptic problems in a fiber reinforced structure. Nonlocal effects[END_REF]) in the study of the reduction of dimension 3d -1d (see [START_REF] Murat | Problèmes monotones dans des cylindres de faible diamètre formés de matériaux hétérogènes[END_REF]). To prove that result, we consider in section 2 below a thin domain Ω ε := εY × (0, L) as a copy of the local cell arising in a periodic fibered medium for which the configuration domain Ω is the union of 1 ε 2 such cells, in other words, Ω := ∪ i∈Iε ( (εY + εi) × (0, L) ) . In fact, in this example the "fibers" play a role of vertical parallel holes and we prove that a strange term already appears at the limit in the 3d -1d problem under exactly the same assumptions on the size of the hole as the ones assumed in the corresponding homogenization problem. Furthermore, it is shown that the structure of the homogenized problem is very close to the one of the limit problem obtained in the reduction of dimension 3d -1d. The homogenization problem will be considered in section 3.

In section 4, we consider another example for which the reduction of dimension 3d -1d leads to a limit problem involving a nonlocal effect; the comparison with the corresponding homogenization problem is made in section 5 and once again the similarity between the two limit problems is pointed out.

One can explain the similarity between the homogenization problem and the reduction of dimension problem by the fact that the geometry of the fibered medium is such that the homogenization process implies a local reduction of dimension so that for such media, the homogenization may be viewed as a repeat of local reductions of dimension. More precisely, when we assume that the domain Ω is made from a single cell of size ε, Ω = Ω ε := Y ε := εY × (0, L) where for example Y is the square defined by

Y := ] -1 2 , 1 2
[ 2 and L > 0, one can denote the variable x in

Ω ε by x = (x ′ , x 3 ) = (εy, x 3 ), y ∈ Y, x 3 ∈ (0, L)
, in such a way that x ′ = εy. This classical dilation transforms the reduction of dimension problem posed in Ω ε into a singular perturbation problem posed on the fixed domain Ω (see [START_REF] Cioranescu | Un terme étrange venu d'ailleurs[END_REF]) and the limit problem (as ε goes to zero) is written in terms of (y, x 3 ). When we study the periodic homogenization with

Ω := ∪ i∈Iε ( (εY + εi) × (0, L) ) , the variable x in Ω writes as x = (x ′ , x 3 ) with x ′ = εy + εi, y ∈ Y, i ∈ I ε , x 3 ∈ (0, L).
It is well known (see [START_REF] Allaire | Homogenization and two-scale convergence[END_REF], [START_REF] Cherednichenko | Nonlocal limits for composite media with highly anisotropic periodic fibers[END_REF]) that the homogenized problem may be formulated in terms of the variables (y, x). Of course, it is not reasonable to expect a homogenized problem involving only x 3 as a macroscopic variable since we have many cells distributed in the plane (x 1 , x 2 ) so that several reductions of dimension occur. For that reason, the role of the part x ′ = (x 1 , x 2 ) of the macroscopic variable in the homogenized problem is restricted to that of a parameter so that the main variables are still (y, x 3 ) as in the 1d-model obtained after the reduction of dimension 3d -1d. Remarkably, dealing with one small cylinder or 1 ε 2 small cylinders does not affect the form of the limit problem.

For the sake of simplicity and brevity we consider in this work the case of the Laplacian but the results remain valid for other more general operators as we will show in coming works.

Strange term in the reduction of dimension.

We consider a thin structure described as follows. Let ε, r ε be two decreasing sequences of positive numbers such that lim ε→0 r ε ε = 0 and let Ωε , Fε , Mε be respectively the parallelepiped, the fiber and the matrix defined by

   Ωε = εY × (0, L), Y = ( ] -1 2 , 1 2 [ ) 2 , L > 0, Fε = r ε D × (0, L) ⊂ Ωε ,
where D(0, r) is the closed disk of radius 0 < r < 1 2 .

(

) 1
For the sake of brevity, in the sequel the interval (0, L) is sometimes denoted by I. We will assume Dirichlet condition on the part

∂ ΩD ε = {(x ′ , x 3 ) ∈ R 3 : x 3 = 0 or x 3 = L or x ′ ∈ r ε ∂D} (2)
of the boundary of Ωε and Newmann condition on the rest

∂ ΩN ε = ∂ Ωε \ ∂ ΩD ε (3)
of the boundary. Remark that ∂ ΩD ε is made of the upper and the lower faces of the parallelepiped together with the boundary of the fiber Fε which plays here the role of a hole. We consider the following equation:

             -∆ū ε = f in Ωε , ūε = 0 on ∂ ΩD ε , ∂ ūε ∂n = 0 on ∂ ΩN ε . ( 4
)
Introduce the change of variables x ′ = εy, u ε (y, x 3 ) := ūε (εy, x 3 ) and set

Ω := Y × (0, L), F ε = 1 ε Fε , Ω ε = Ω \ F ε , ( 5
)
H 1 D (Ω ε ) := {u ∈ H 1 (Ω), u(y, 0) = u(y, L) = 0, u = 0 on ∂F ε } (6)
. Finally, denote by ∇ ′ the gradient with respect to the two first variables y = (y 1 , y 2 ). Similarly, ∆ ′ denotes the Laplacian with respect to the same variables.

With these notations, one deduces easily that (4) is transformed into an equation satisfied by u ε with the following variational formulation

       u ε ∈ H 1 D (Ω ε ), ∫ Ωε ( 1
ε 2 ∇ ′ u ε ∇ ′ ϕ + ∂u ε ∂x 3 ∂ϕ ∂x 3 ) dydx 3 = ∫ Ωε f ϕ dydx 3 , ∀ ϕ ∈ H 1 D (Ω ε ). (7)
Assuming the source term f ∈ L 2 (Ω), problem ( 7) is well-posed; in addition, u ε may be extended by zero inside the hole F ε so that the new sequence still denoted by [START_REF] Caillerie | A perturbation problem with two small parameters in the framework of the heat conduction of a fiber reinforced body[END_REF], we obtain easily that there exist u ∈ H 1 0 (0, L) (subspace of functions in H 1 D (Ω) depending only on the variable x 3 ) and a subsequence (still denoted by ε) such that:

u ε belongs to H 1 D (Ω) := {u ∈ H 1 (Ω), u(y, 0) = u(y, L) = 0}. Furthermore, taking ϕ = u ε in
u ε ⇀ u weakly in H 1 D (Ω). (8)
The limit function u depends on the size of the hole and it is characterized by the following theorem. Indeed, defining f (x 3 ) :=

∫ Y f (y, x 3
) dy, we have the following:

Theorem 2.1.                                              If there exists k > 0 such that r ε ∼ e -k ε 2 , the function u
is the unique solution of the one-dimensional problem u ∈ H 1 0 (0, L), -

d 2 u dx 2 3 + 2πr k u = f in (0, L), if r ε ≪ e -k ε 2 , u is the unique solution of u ∈ H 1 0 (0, L), - d 2 u dx 2 3 = f in (0, L), if r ε ≫ e -k ε 2 , then u = 0. ( 9
)
Sketch of the Proof. Introduce the following function (see [START_REF] Cioranescu | Un terme étrange venu d'ailleurs[END_REF])

w ε (y) =            0 in rε ε D, ln(|y|)-ln( rε ε ) ln(r)-ln( rε ε ) in D \ rε ε D, 1 in Y \ D. ( 10
)
Then:

∂w ε ∂n = 1 ln( ε rε ) on ∂D, -∆ ′ w ε = 0 in D\ r ε ε D, w ε ⇀ 1 weakly in H 1 (Y ).
Use as a test function in [START_REF] Caillerie | A perturbation problem with two small parameters in the framework of the heat conduction of a fiber reinforced body[END_REF]:

ϕ(y, x 3 ) = w ε (y)ū(x 3 ), ū ∈ H 1 0 (I).Then ∫ L 0 ∫ D\ rε ε D 1 ε 2 ∇ ′ u ε ∇ ′ w ε ūdydx 3 + ∫ Ω ∂u ε ∂x 3 w ε dū dx 3 dydx 3 = ∫ Ω f w ε ū dydx 3 . ( 11
)
Integrating by parts with respect to y the first integral on the left, and bearing in mind the property -∆ ′ w ε = 0 in D \ rε ε D, we see that it reduces to the boundary integral:

∫ L 0 ∫ ∂D 1 ε 2 ∂w ε ∂n u ε dσūdx 3 = ∫ I ∫ ∂D 1 ε 2 1 ln( ε rε ) u ε dσūdx 3 .
Hence, we have the following conclusions in terms of the thickness of the hole:

if r ε = e -k ε 2 , the limit of the last integral is equal to

1 k |∂D| ∫ L 0 uū dx 3 = 2πr k ∫ L 0
uū dx 3 , so that the limit equation is the first equation of ( 9).

If r ε ≪ e -k ε 2 , the limit of that integral is equal to zero in such a way that the extra term disappears and the limit equation is now the second equation arising in [START_REF] Casado-Diaz | Two-scale convergence for nonlinear Dirichlet problems[END_REF].

Finally, let us prove that u = 0 under the assumption r ε ≫ e -k ε 2 . It is well known that for given ε > 0, the eigenvalues of the operator -

1 ε 2 ∆ ′ - ∂ 2 ∂x 2 3 with eigenvectors in H 1 D (Ω ε
) is an increasing positive sequence of real numbers and that for f ∈ L 2 (Ω), the solution of ( 7) is given by

u ε = ∞ ∑ i=1 1 λ ε i (f, v ε i )v ε i where (λ ε i , v ε i ) i≥1
is the sequence of the eigenelements under consideration and where (, ) denotes the scalar product of L 2 (Ω). As a consequence, we deduce the estimate

∥ u ε ∥ L 2 (Ω) ≤ 1 λ ε 1 ∥ f ∥ L 2 (Ω) . ( 12
)
On the other hand, the first eigenvalue is given by λ ε 1 = min

v∈H 1 D (Ωε)\{0} ∫ Ωε ( 1
ε 2 |∇ ′ v| 2 + | ∂v ∂x 3 | 2 ) dx ∫ Ωε |v| 2 dx.
. As an immediate consequence, we derive the following inequality

λ ε 1 ≥ min v∈H 1 D (Ωε)\{0} ∫ L 0 ∫ Y \Fε 1 ε 2 |∇ ′ v| 2 dydx 3 ∫ L 0 ∫ Y \Fε |v| 2 dydx 3 . ( 13
)
On the other hand, since for

x 3 ∈ (0, L) any function v ∈ H 1 D (Ω ε ) is such that v(., x 3 ) ∈ H 1 (Y ) and v(., x 3 ) = 0 on ∂( rε ε D)
, we can use the following estimate proved in [START_REF] Rauch | Potential and scattering theory on wildly perturbed domains[END_REF], pages 44-45:

min v∈H 1 (Y )\{0}, ∫ Y \Fε 1 ε 2 |∇ ′ v| 2 dy ∫ Y \Fε |v| 2 dy ≥ 1 ε 2 1 ∫ 1 2 rε ε t dt ∫ 1 2 rε ε 1 t dt . ( 14
)
An easy computation shows that the quantity in the right hand side of (32) goes to

+∞ if r ε ≫ e -k ε 2
and then the result follows from ( 13) and ( 12).

3. Strange term in the homogenization. In order to make easy the comparison with the results of the last section, we keep here analogous notations; in particular, Y × (0, L) and D(0, r) denote respectively the parallelepiped and the disk defined in the previous section. We assume now that Ω is the configuration domain of a set of parallelepipeds Y i ε := (εY +εi)×(0, L) distributed with a period εY in the horizontal x ′ -directions; each parallelepiped Y i ε contains another small parallelepiped playing the role of a hole

F i ε := (r ε D + εi) × (0, L).
The set of all the holes contained in Ω will be denoted by F ε . We then define

Ω ε = Ω \ F ε , F ε = ∪ i∈Iε F i ε , ( 15
)
where I ε denotes as usual

I ε := {i = (i 1 , i 2 ) ∈ Z 2 , Y i ε ⊂ Ω}.
Note that the analogous configuration domain in the reduction of the dimension problem is (εY \ r ε D) × (0, L) denoted there by Ωε (to be distinguished from Ω ε denoting there the domain obtained after scaling). In other words, the real domain (before scaling) of the 3d -1d problem is nothing but the representative cell of the present homogenization problem.

The equation we want to homogenize is the following

u ε ∈ H 1 D (Ω ε ), ∫ Ωε ∇u ε ∇ϕ dx = ∫ Ωε f ϕ dx, ∀ ϕ ∈ H 1 D (Ω ε ), (16)
where

H 1 D (Ω ε ) := {u ∈ H 1 (Ω), u(x ′ , 0) = u(x ′ , L) = 0, u = 0 on ∂F ε }.
To give the limit problem, we need to define

Γ D := {x = (x ′ , x 3 ) ∈ Ω, x 3 = 0 or x 3 = L}, Γ N := ∂Ω \ Γ D , ( 17
)
together with the space

H 1 D (Ω) := {u ∈ H 1 (Ω), u = 0 on Γ D }. ( 18
)
We consider a sequence δ ε of positive numbers such that

r ε ≪ δ ε ≪ ε, lim ε→0 ε 2 ln(δ ε ) = 0. (19)
For example if r ε = e -k ε 2 , k > 0, one can take δ ε = e -k ε . Denote by C δε the circle of radius δε ε centered at the origin and set C i δε = εC δε + εi. Finally, we introduce the sequence

ũε = ∑ i∈Iε 1 2πδ ε ∫ C i δε u ε dσ χ Y i ε (x ′ ). (20)
The following estimate, proved in [START_REF] Bellieud | Homogenization of elliptic problems in a fiber reinforced structure. Nonlocal effects[END_REF], will be helpful in the sequel.

Lemma 3.1. There exists a constant

C > 0 such that ∫ Ω |u ε -ũε | 2 dx ≤ Cε 2 ( 1 + | ln( δ ε √ 2 ε )| ) ∫ Ω |∇u ε | 2 dx. ( 21
)
We now state the main result of this section through the following theorem.

Theorem 3.2.

There exists u ∈ H 1 D (Ω) such that the sequence of solutions of ( 16) may be extended by zero to H 1 D (Ω) on a sequence still denoted by u ε such that

u ε ⇀ u weakly in H 1 D (Ω). ( 22
)
Furthermore, u is characterized by

                                                                             If there exists k > 0 such that r ε ∼ e -k ε 2 ,
then u is the unique solution of the homogenized problem

u ∈ H 1 D (Ω), -∆u + 2π k u = f in Ω, ∂u ∂n = 0 on Γ N . If r ε ≪ e -k ε 2 , then u is the unique solution of u ∈ H 1 D (Ω), -∆u = f in Ω, ∂u ∂n = 0 on Γ N . If r ε ≫ e -k ε 2 , then u = 0. ( 23
)
,

Proof. First of all, we remark that the sequence u ε of solutions of ( 16) extended by zero to the holes is bounded in H 1 D (Ω); this can be seen using ϕ = u ε as test function in [START_REF] Sili | Homogenization of a nonlinear monotone problem in an anisotropic medium[END_REF] and then applying the Cauchy-Schwarz inequality in the right hand side. Therefore, we can assume that there exists u ∈ H 1 D (Ω) such that ( up to a subsequence) u ε ⇀ u weakly in H 1 D (Ω). We want to find the equation satisfied by u.

Let ϕ ∈ C ∞ (Ω) ∩ H 1 D (Ω). Define the following sequence of piece-wise constant functions

ϕ ε = ∑ i∈Iε 1 πr 2 ε ∫ D i rε ϕ(x ′ , x 3 ) dx ′ χ Y i ε (x ′ ), (24)
where D i rε denotes the disk of center ε(i 1 , i 2 ) and of radius r ε . Using the regularity of ϕ, one can check easily by the use of the mean value Theorem that there exists a constant C > 0 such that |ϕ -

ϕ ε | ≤ Cε in Ω. Define the function d ε (x ′ ) := dist ( x ′ ; {εi, i = (i 1 , i 2 ) ∈ Z 2 } )
and denote by D i δε the disk with center εi and with radius δ ε , then we will use the following function defined for each i ∈ Z 2 such that Y i ε ⊂ Ω (i.e., i ∈ I ε ) by:

w ε (x ′ ) =            0 in D i rε , ln(dε(x ′ ))-ln(rε) ln(δε)-ln(rε) in D i δε \ D i rε , 1 in Y i ε \ D i δε . ( 25
)
In view of the definition of w ε , one can check easily by the use of cylindrical coordinates in the tube

U i ε := {(x ′ , x 3 ) ∈ Y i ε , r ε < d ε (x ′ ) < δ ε , x 3 ∈ (0, L)} the following estimate: ∫ Ω |∇ ′ w ε (x ′ )| 2 dx ≤ C ε 2 ln( δε rε ) . ( 26
)
Hence, the second equality arising in (19) allows us to conclude that w ε is bounded in H 1 (Ω) under the hypothesis r ε ≈ e -k ε 2 or the hypothesis r ε ≪ e -k ε 2 . Note also that the latest hypothesis which corresponds to "small holes" implies that the sequence w ε which is bounded in H 1 (Ω) strongly converges to 1 in H 1 (Ω). If we assume only the hypothesis r ε ≈ e -k ε 2 , then w ε only converges to 1 in the weak topology of H 1 (Ω), (see [START_REF] Casado-Diaz | Two-scale convergence for nonlinear Dirichlet problems[END_REF]).

We now choose a test function in the form ϕ(x)w ε (x ′ ) in equation ( 16) with ϕ ∈ H 1 D (Ω) and we get In view of the definition of w ε , that integral reduces to the integral on the tube

         ∫ Ω ∇u ε ∇ ( w ε ϕ ) dx = ∫ Ω ∇u ε ∇w ε (ϕ -ϕ ε ) dx + ∫ Ω w ε ∇u ε ∇ϕ dx + + ∫ Ω ϕ ε ∇u ε ∇w ε dx.
U i ε := {(x ′ , x 3 ) ∈ Y i ε , r ε < d ε (x ′ ) < δ ε , x 3 ∈ (0, L)}.
Using cylindrical coordinates (r, θ, x 3 ) ∈ (r ε , δ ε ) × (0, 2π) × (0, L) and the definition (25) of w ε , it is not difficult to check that one can write in terms of (r, θ, x 3 )

∇ ′ u ε ∇ ′ w ε = 1 r ln( δε rε ) ∂u ε ∂r (r, θ, x 3 ) in U i ε , ( 28
)
so that bearing in mind the formula (24 ) and the fact that u ε (r ε , θ, x 3 ) = 0, we get:

                                                                     ∫ Ω ϕ ε ∇u ε ∇w ε dx = ∑ i∈Iε ∫ U i ε ϕ ε ∇u ε ∇w ε dx = ∑ i∈Iε ∫ L 0 ∫ 2π 0 ∫ δε rε 1 r ln( δε rε ) ∂u ε ∂r (r, θ, x 3 )ϕ ε (x 3 ) r dr dθdx 3 = 1 ln( δε rε ) ∑ i∈Iε ∫ L 0 ∫ 2π 0 u ε (δ ε , θ, x 3 )ϕ ε (x 3 ) dθdx 3 = 2π ln( δε rε ) ∑ i∈Iε ∫ L 0 ∫ 2π 0 1 2π u ε (δ ε , θ, x 3 )ϕ ε (x 3 ) dθdx 3 = 2π ε 2 ln( δε rε ) ∑ i∈Iε ∫ L 0 ∫ Y i ε ũε (x ′ , x 3 )ϕ ε (x) dx ′ dx 3 = 2π ε 2 ln( δε rε ) ∫ Ω ũε (x ′ , x 3 )ϕ ε (x) dx ′ dx 3 , ( 29
)
where ũε is defined by (20). On the other hand, from Lemma 3.1 and the property lim ε→0 δ ε = 0 we deduce that ũε → u strongly in L 2 (Ω) since u ε converges strongly to u in L 2 (Ω). Finally using the property ϕ ε → ϕ strongly in L 2 (Ω), we obtain the first two items of the theorem by passing to the limit in the last term of (29). We now prove the last part of the theorem using the same idea as in the 3d -1d problem showing that the first eigenvalue λ 1 ε of the operator under consideration tends to infinity when ε → 0. We establish a lower bound of λ 1 ε as follows

λ ε 1 = min v∈H 1 D (Ωε)\{0} ∫ Ωε ( |∇ ′ v| 2 + | ∂v ∂x 3 | 2 ) dx ∫ Ωε |v| 2 dx ≥ min v∈H 1 D (Ωε)\{0} ∫ Ωε |∇ ′ v| 2 dx ∫ Ωε |v| 2 dx . (30) In each cell Y i ε ⊂ Ω, one can write for all v ∈ H 1 D (Ω ε ) by the use of cylindrical coordinates as above ∫ Y i ε |∇ ′ v| 2 dx ≥ ∫ L 0 ∫ 2π 0 ∫ ε rε | ∂v ∂r | 2 r dr dθ dx 3 . ( 31
)
As a consequence, we infer

∫ Y i ε |∇ ′ v| 2 dx ∫ Y i ε |v| 2 dx ≥ ∫ L 0 ∫ 2π 0 ∫ ε rε | ∂v ∂r | 2 r dr dθ dx 3 ∫ L 0 ∫ 2π 0 ∫ ε rε |v| 2 r dr dθ dx 3 . ( 32
)
Using the following inequality proved in [START_REF] Rauch | Potential and scattering theory on wildly perturbed domains[END_REF],

∫ ε rε | ∂v ∂r | 2 r dr ≥ 1 γ(ε) ∫ ε rε |v(r, θ, x 3 )| 2 r dr, a.e. (θ, x 3 ) ∈ (0, 2π) × (0, L), ( 33
)
where γ(ε) :=

∫ ε rε r dr ∫ ε rε 1
r dr, we get by integrating over (0, 2π) × (0, L) and summing up over i ∈ I ε ,

min v∈H 1 D (Ωε)\{0} ∫ Ωε |∇ ′ v| 2 dx ∫ Ωε |v| 2 dx ≥ 1 γ(ε) . ( 34
)
One can check that the constant γ(ε) goes to zero under the hypothesis

r ε ≫ e -k ε 2
and then the conclusion of the last part of the theorem is a consequence of the inequality ( 12) which still holds true in the homogenization problem.

We now give the second example regarding nonlocal effects in the limit problem. We begin with the reduction of dimension.

Nonlocal effects in the reduction of dimension.

In order to describe the heterogeneities of the medium, we need here to introduce some other notations and change slightly those of section 2. In this section D(0, r) still denotes the disk defined in the previous sections, Ω ε is defined by Ω ε = εY ×(0, L) = εY ×I, Y being defined by (1) and I := (0, L). In addition, we introduce the sets

M ε = Ω ε \ F ε where now F ε := εD × (0, L) in such a way that Ω ε = M ε ∪ F ε . Similarly, we put F := D × I, M := Ω \ F, in such a way that Ω = Y × (0, L) = M ∪ F. We still denote by H 1 D (Ω) the space H 1 D (Ω) := {u ∈ H 1 (Ω), u(y, 0) = u(y, L) = 0}. We then consider the problem                u ε ∈ H 1 D (Ω), ∫ Ω ( χ F + ε 2 χ M )( 1
ε 2 ∇ ′ u ε ∇ ′ ϕ + ∂u ε ∂x 3 ∂ϕ ∂x 3 ) dydx 3 = ∫ Ω f ϕ dydx 3 , ∀ ϕ ∈ H 1 D (Ω). (35)
Clearly, problem (35) is the variational formulation in the fixed domain Ω of the equation satisfied by u ε (y, x 3 ) := ūε (εy, x 3 ) and obtained from the equation posed in the variable domain Ω ε and satisfied by ūε . Loosely speaking, equation (35) takes into account the contrast between the diffusion in the fiber which is assumed to be of order 1 (before scaling) and the diffusion in the matrix (outside the fiber) assumed to be of order ε 2 . As regards of the asymptotic behaviour of the sequence u ε , the main result may be stated as follows.

Theorem 4.1. The sequence u ε of solutions of ( 35) is such that:

             u ε → u 0 (x) + v(x 3 ) strongly in L 2 (I; H 1 (Y )), 1 ε ∇ ′ u ε χ F → 0 strongly in ( L 2 (Ω) ) 2 , ∂uε ∂x3 χ F → dv dx3 χ F in L 2 (Ω), ∇ ′ u ε χ M → ∇ ′ u 0 χ M strongly in ( L 2 (Ω) ) 2 , ε ∂uε ∂x3 χ M → 0 in L 2 (Ω), (36)
where the pair (u 0 , v) is the unique solution of

               (u 0 , v) ∈ {ϕ ∈ L 2 (I; H 1 (Y )), ϕ(., x 3 ) = 0 in D} × H 1 0 (I), ∫ Ω ( dv dx 3 dv dx 3 χ F + ∇ ′ u 0 ∇ ′ ūχ M ) dydx 3 = ∫ Ω f (ū + v) dydx 3 , ∀ (ū, v) ∈ {ϕ ∈ L 2 (I; H 1 (Y )), ϕ = 0 in D} × H 1 0 (I). ( 37
)
Furthermore, the convergence of the energies

E ε := ∫ Ω ( χ F + ε 2 χ M )( 1
ε 2 ∇ ′ u ε ∇ ′ u ε + ∂u ε ∂x 3 ∂u ε ∂x 3 ) dydx 3 towards the limit energy E 0 := ∫ Ω ( dv dx 3 dv dx 3 χ F +∇ ′ u 0 ∇ ′ u 0 χ M ) dydx 3 holds true.
Proof. Taking ϕ = u ε in equation ( 35), we get

E ε = ∫ Ω f u ε dydx 3 ≤ C ∥ u ε ∥ L 2 (Ω) . ( 38
)
By the assumption u ε (y, L) = u ε (y, 0) = 0 one can apply the one dimensional Poincaré inequality in (0, L) to the function u ε (y, .) for a given y ∈ Y to obtain

∫ L 0 |u ε (y, x 3 )| 2 dx 3 ≤ C ∫ L 0 | ∂u ε ∂x 3 (y, x 3 )| 2 dx 3 , a.e. in Y. ( 39
)
Integrating (39) with respect to y ∈ D and bearing in mind that F = D × (0, L), we get

∥ u ε ∥ 2 L 2 (F ) ≤ C ∥ ∂u ε ∂x 3 χ F ∥ 2 L 2 (Ω) ≤ CE ε . ( 40
)
On the other hand, for a given x 3 ∈ (0, L) the Poincaré-Wirtinger inequality gives the estimate

∥ u ε (., x 3 ) - 1 |D| ∫ D u ε (., x 3 ) dy ∥ L 2 (Y ) ≤ C ∥ ∇ ′ u ε (., x 3 ) ∥ L 2 (Y ) . ( 41
)
Integrating (41) with respect to x 3 ∈ (0, L), we derive with the help of (40) and the fact that

∥ ∇ ′ u ε ∥ 2 L 2 (Ω) ≤ E ε for sufficiently small ε, the estimate ∥ u ε ∥ 2 L 2 (Ω ≤ CE ε . ( 42
)
Turning back to equation ( 35) and applying the Young inequality in the right hand side, we deduce that E ε ≤ C which in turn implies by virtue of (42) that

∥ u ε ∥ L 2 (Ω ≤ C. ( 43
)
From estimate (43) and the boundedness of the sequence E ε , we deduce, up a subsequence still denoted by ε, the existence of

u 0 ∈ L 2 (I; H 1 (Y )) such that u ε ⇀ u 00 (y, x 3 ) weakly in L 2 (I; H 1 (Y )), (44)
which together with the estimates (which are consequences of

E ε ≤ C)      ∥ 1 ε ∇ ′ u ε χ F ∥ (L 2 (Ω)) 2 ≤ C, ∥ u ε ∥ H 1 (F ) ≤ C, ( 45
)
easily implies the existence of v ∈ H 1 0 (I) such that u 00 (y,

x 3 ) = v(x 3 ) in F. ( 46
)
Therefore we obtain all the convergences stated in (36) except the convergence 1 ε ∇ ′ u ε χ F ⇀ 0 replacing the strong convergence by the weak convergence, if we define u 0 by u 0 (y, x 3 ) = u 00 (y, x 3 ) -v(x 3 ) in Ω.

(47) To prove that the weak convergences are actually strong convergences, we first identify the limit problem and then we will prove the convergence of the sequence E ε to E 0 . The strong convergence of u ε in L 2 (Ω) will be derived from the convergence of the sequence E ε with the help of a kind of Poincaré-Wirtinger inequality as it will be seen below.

The first estimate in (45) amounts to say that the sequence defined by

w ε := 1 ε (u ε - 1 |D| ∫ D u ε dy) is bounded in L 2 (I; H 1 m (D))
where H 1 m (D) is the subspace of functions in H 1 (D) with zero average. Hence, one can assume that for a subsequence at least, it converges weakly in L 2 (I; H 1 m (D)) to some w. Taking a test function ϕ in (35) in the form ϕ = ū + v + ε w with ū ∈ D(Ω) such that ū = 0 in F and v ∈ H 1 0 (I), w ∈ D(Ω), we can pass to the limit to find the equation

∫ Ω ( (∇ ′ w∇ ′ w + dv dx 3 dv dx 3 )χ F + ∇ ′ u 0 ∇ ′ ūχ M ) dydx 3 = ∫ Ω f (ū + v) dydx 3 . ( 48
)
By a density argument we can choose ū = v = 0 and w = w in (48) so that we get w = 0 and equation ( 37) is obtained. Note that by the same, the convergence

1 ε ∇ ′ u ε χ F ⇀ 0 is proved.
Introduce now the sequence

X ε := ∫ Ω ( ( 1 ε 2 |∇ ′ u ε | 2 + ( ∂u ε ∂x 3 - dv dx 3 ) 2 )χ F + (|∇ ′ u ε -∇ ′ u 0 | 2 + ε 2 | ∂u ε ∂x 3 | 2 )χ M ) . (49)
Choosing ϕ = u ε in (35) and (ū, v) = (u 0 , v) in (37) and thanks to the previous weak convergences, we show that the limit of X ε is zero so that the stated strong convergences take place. To prove the strong convergence in L 2 (Ω) of the sequence u ε , we will use the following Poincaré-Wirtinger type inequality: there exists a positive constant C such that

     ∥ u ∥ L 2 (Ω) ≤ C ( ∥ ∇ ′ u ∥ L 2 (Ω) + ∥ ∂u ∂x 3 ∥ L 2 (F ) ) , ∀ u ∈ L 2 (I; H 1 (Y )) ∩ L 2 (D; H 1 0 (I)). (50)
The proof may be done arguing by contradiction. Assume that there exists a sequence u n in that space such that ∥ u n ∥ L 2 (Ω) = 1 for all n while the sequence

∥ ∇ ′ u n ∥ L 2 (Ω) + ∥ ∂un ∂x3 ∥ L 2 (
F ) goes to zero. Then the bidimensional Poincar é-Wirtinger inequality applied to the function u n (., x 3 ) for x 3 ∈ I together with an integration over I of such inequality allows one to get the estimate

∥ u n - 1 |D| ∫ D u n dy ∥ L 2 (Ω) ≤ C ∥ ∇ ′ u n ∥ L 2 (Ω) , ∀n. (51)
On the other hand, the one dimensional Poincar é Inequality applied to the function u n (y, .) for y ∈ D and then an integration with respect to y ∈ D of that inequality lead to the estimate

∥ 1 |D| ∫ D u n dy ∥ L 2 (Ω) ≤ C ∥ ∂u n ∂x 3 ∥ L 2 (Ω) , ∀n. (52)
From ( 51) and (52), we deduce that ∥ u n ∥ L 2 (Ω) goes to zero and this is a contradiction with the assumption. Applying (50) to the sequence u ε -

(u 0 + v), (recall that u 0 = 0 in F ) we get            ∥ u ε -(u 0 + v) ∥ L 2 (Ω) ≤ C ( ∥ ∇ ′ u ε -∇ ′ u 0 ∥ L 2 (M ) + + ∥ ∇ ′ u ε ∥ L 2 (F ) + ∥ ∂uε ∂x3 -dv dx3 ∥ L 2 (F )
) .

(

) 53
Applying the previous strong convergences, the right hand side of (53) tends to zero so that the proof of the theorem is now complete.

One can highlight the nonlocal character of the previous equation. Indeed, define

u(x 3 ) := ∫ Y ( u 0 (x) + v(x 3 ) ) dy = ∫ Y \D u 0 (y, x 3 ) dy + v(x 3 ). Let û be the unique solution of              -∆ ′ û = 1 in Y \ D, û = 0 on ∂D, ∂ û ∂n = 0 on ∂Y. ( 54
) Define m := ∫ Y \D
û dy > 0. Then, we have the following result.

Theorem 4.2. Assume that f (y, x 3 ) = f (x 3 ) does not depend on the variable y.

Then u 0 given in (36) may be written as u 0 (x) = f (x 3 )û(y) and the sequence ūε defined in the variable domain Ω ε = εY × I (recall that u ε (y, x 3 ) = ūε (εy, x 3 ) for (y,

x 3 ) ∈ Ω) is such that                              1 |εY | ∫ εY ūε (x ′ , x 3 ) dx ′ → u strongly in L 2 (I), 1 |εM | ∫ εM ε∇ ′ ūε (x ′ , x 3 ) dx ′ → 1 |M | ∫ M ∇ ′ u 0 (y, x 3 ) dy strongly in L 2 (I), 1 |εD| ∫ εD ∂ ūε ∂x 3 (x) dx ′ → dv dx 3 strongly in L 2 (I). (55)
The pair (u, v) is the unique solution of the nonlocal one-dimensional problem

       (u, v) ∈ L 2 (I) × H 1 0 (I), u(x 3 ) -v(x 3 ) = mf (x 3 ), -|D| d 2 v dx 2 3 = f (x 3 ) in I. ( 56
)
Proof. Making the change of variable x ′ = εy, the strong convergences (55) become immediate consequences of the convergences (36).

The second equation of ( 56) is obtained from the equation (37) by choosing ū = 0 and taking into account the fact that F = D × I.

On the other hand, one can check that the function f (x 3 )û(y) where û is the solution of (54), solves the same equation as u 0 , that is the equation obtained from (37) by choosing v = 0. By virtue of the uniqueness of u 0 , we conclude that u 0 (y, x 3 ) = f (x 3 )û(y) and then the first equation of ( 56) is nothing but the equality u(x 3 ) := ∫ Y \D u 0 (y, x 3 ) dy + v(x 3 ) which defines u.

Nonlocal effects in the homogenization.

To describe the geometry of the medium, we need further notations.

Let D(0, r) be the disk defined in the previous sections and let ω be the square ω :=] -1, 1[ 2 . Assume that Ω := ω × (0, L) = ω × I is now the configuration domain of a set F ε of cylindrical parallel fibers periodically distributed with a period

εY = ε ( ] -1 2 , 1 2 [
) 2 in the x ′ -horizontal directions which are surrounded by a poor conductor occupying the matrix M ε in such a way that

Ω = F ε ∪ M ε , F ε = ∪ i∈Iε F i ε , F i ε = (ε D(0, r) + εi) × I (57) M ε = Ω \ F ε . (58)
Hence, the medium is now a periodic replication of the one arising in the previous section. The equation we want to homogenize is the following

u ε ∈ H 1 D (Ω), ∫ Ω ( χ Fε + ε 2 χ Mε ) ∇u ε ∇ϕ dx = ∫ Ω f (x)ϕ(x) dx, ∀ ϕ ∈ H 1 D (Ω), (59)
where H 1 D (Ω) is still the space defined in the beginning of the previous section. For the sake of brevity, we consider only the case of a source term not depending on the microscopic variable but one can handle also that general case as pointed out in the Remark 5.1 below.

Before stating the main result, we recall the definition of two-scale convergence (see [START_REF] Nguetseng | A general convergence result for a functional related to the theory of homogenization[END_REF], [START_REF] Allaire | Homogenization and two-scale convergence[END_REF]) a well adapted tool for periodic homogenization. A sequence

t ε ∈ L 2 (Ω) two scale converges to a function t ∈ L 2 (Ω × Y ) if ∫ Ω t ε (x)ϕ(x, x ′ ε ) dx → ∫ Ω ∫ Y t(x, y)ϕ(x, y) dxdy, ∀ ϕ ∈ L 2 (Ω; C # (Y ))
where C # (Y ) denotes the space of functions which are continuous and Y -periodic. It is known that every bounded sequence in L 2 (Ω) admits a two-scale converging subsequence.

In the sequel, the notation ⇀⇀ will stand for the two-scale convergence. The main result may be stated as follows.

Theorem 5.1. The sequence u ε of solutions of ( 59) is such that :

u ε ⇀⇀ u 0 (x, y) + v(x), ( 60
)
1 ε ∇ ′ u ε χ Fε ⇀ 0 weakly in L 2 (Ω), (61)
∂u ε ∂x 3 χ Fε ⇀ |D| ∂v ∂x 3 weakly in L 2 (Ω), ( 62
)
ε∇ ′ u ε χ Mε ⇀⇀ ∇ ′ y u 0 χ Y \D , ε ∂u ε ∂x 3 χ Mε ⇀⇀ 0, ( 63
)
where the pair

(u 0 , v) ∈ {ϕ ∈ L 2 (Ω; H 1 # (Y )), ϕ(x, .) = 0 in D} × L 2 (ω; H 1 0 (I)) is the unique solution of        ∫ Ω×Y ( ∂v ∂x 3 ∂v ∂x 3 χ D (y) + ∇ ′ u 0 ∇ ′ ūχ ω\D ) dxdy = ∫ Ω×Y f (x)(ū + v) dxdy, ∀ (ū, v) ∈ {ϕ ∈ L 2 (Ω; H 1 # (Y )), ϕ(x, .) = 0 in D} × L 2 (ω; H 1 0 (I)). (64)
Furthermore, by the same approach already used in the 3d -1d problem, one can eliminate the microscopic variable y to derive a formulation of the limit problem involving only the macroscopic variable x. Indeed, let û be the unique solution of

           -∆ ′ y û = 1 in Y \ D, û = 0 on ∂D, û is Y -periodic. ( 65
) Define m := ∫ Y \D û dy > 0.
Theorem 5.2. The function u 0 given in (60) may be written as u 0 (x, y) = f (x)û(y) and the sequence u ε is such that

u ε (x) ⇀ u := ∫ Y \D u 0 (x, y) dy + v(x) in L 2 (Ω), (66)
the pair (u, v) is the unique solution of the nonlocal homogenized problem

       (u, v) ∈ L 2 (Ω) × L 2 (ω; H 1 0 (I)), u(x) -v(x) = mf (x), -|D| ∂ 2 v ∂x 2 3 = f (x) in Ω. (67)
In the homogenization setting we cannot reasonably expect strong convergences similar to those obtained in the reduction of dimension due to the oscillations induced by the microscopic variable. However, we have the following result. Theorem 5.3. Let u ε be the sequence of solutions of (59) and let (u 0 , v) be the pair defined by (64) , then the sequence

X ε = ∫ Ω ( (| ∂v ∂x 3 - ∂u ε ∂x 3 | 2 +| 1 ε ∇ ′ u ε | 2 )χ D ( x ′ ε )+(|∇ ′ y u 0 (x, x ′ ε )-∇ ′ u ε | 2 +|ε ∂u ε ∂x 3 | 2 )χ Y \D ( x ′ ε ) (68) converges to zero.
Remark that this result is a corrector result meaning that u ε may be approximated by v(x) + u 0 (x, x ′ ε ). Due to the close similarities between the reduction of dimension problem and the homogenization problem and for the sake of brevity, we will give only the outline of the proof of Theorem 5.1.

Proof. As above we take ϕ = u ε in the equation ( 59) and then we prove in the same way that the L 2 -norm of u ε is dominated up to a positive constant by the energy which is now defined by

E ε := ∫ Ω ( χ Fε + ε 2 χ Mε )( ∇ ′ u ε ∇ ′ u ε + ∂u ε ∂x 3 ∂u ε ∂x 3 ) dx ′ dx 3 .
Once again, we use the Poincaré-Wirtinger inequality in the reference cell Y to obtain the inequality

∫ Y |u - 1 |D| ∫ D u dy| 2 dy ≤ C ∫ Y |∇u| 2 dy, ∀u ∈ H 1 (Y ). (69)
For given ε, i ∈ I ε and x 3 ∈ (0, L), we choose u(y) := u ε (εy + εi, x 3 ) in ( 69) and then we make the change of variables x ′ = εy + εi in such a way we derive the inequality ∫

Y i ε |u ε - 1 |D i ε | ∫ D i ε u ε dx ′ | 2 dx ′ ≤ C ∫ Y i ε ε 2 |∇ ′ u ε | 2 dx ′ . ( 70
)
On the other hand, similarly to the proof in the reduction of dimension problem and due to the Dirichlet boundary condition u ε (x ′ , 0) = u ε (x ′ , L) = 0, one can apply the one dimensional Poincaré inequality to the sequence u ε (x ′ , .) for almost all x ′ ∈ ω to get after an integration with respect to x ′ and after summing up over i

∈ I ε , ∫ Ω ∑ i∈Iε | 1 |D i ε | ∫ D i ε u ε dx ′ | 2 dx ≤ C ∑ i∈Iε ∫ L 0 ∫ D i ε | ∂u ε ∂x 3 | 2 dx = C ∫ Ω | ∂u ε ∂x 3 | 2 χ Fε dx.
(71) Summing up over i ∈ I ε and integrating over (0, L) the inequality (70), we get with the helph of ( 71)

∫ Ω |u ε | 2 dx ≤ C ∫ Ω ( | ∂u ε ∂x 3 | 2 χ Fε + ε 2 |∇ ′ u ε | 2 ) dx ≤ CE ε . ( 72
)
Using (72) in equation ( 59) in which we take ϕ = u ε , we derive easily the estimate E ε ≤ C. In particular the last inequality implies that the sequence ∂u ε ∂x 3 χ Fε is bounded in L 2 (Ω) so that u ε χ Fε is bounded in L 2 (ω; H 1 0 (0, L)). Therefore one can extract a subsequence still denoted ε and find functions u 00 ∈ L 2 (Ω; H 1 # (Y )), K ∈ (L 2 (Ω × Y )) 2 , v ∈ L 2 (ω; H 1 0 (0, L)) such that the following convergences hold true:

u ε ⇀⇀ u 00 , ε∇ ′ u ε ⇀⇀ ∇ ′ y u 00 , 1 ε ∇ ′ u ε χ Fε ⇀⇀ Kχ D (y), (73)
Since ∇ ′ u ε χ Fε strongly converges to zero in L 2 (Ω), we deduce from (73) that u 00 (x, y) = v(x) in Ω × D with some v ∈ L 2 (Ω).

In addition, there exists v ∈ L 2 (ω; H 1 0 (0, L)) such that u ε χ Fε ⇀ v weakly in L 2 (ω; H 1 0 (0, L))),

∂u ε ∂x 3 χ Fε ⇀ ∂v ∂x 3 weakly in L 2 (Ω). (74)
From the first convergence of (73) we deduce v(x) := ∫ Y u 00 (x, y)χ D (y) dy = |D|v(x). Hence, v ∈ L 2 (ω; H 1 0 (0, L)). On the other hand, the Poincaré-Wirtinger inequality applied in H 1 (D) and the same change of variables as in the proof of (70) leads to

∫ D i ε |u ε - 1 |D i ε | ∫ D i ε u ε dx ′ | 2 dx ′ ≤ C ∫ D i ε ε 2 |∇ ′ u ε | 2 dx ′ . ( 75
)
Using the boundedness in L 2 (Ω) of the sequence 1 ε ∇ ′ u ε χ Fε , we deduce from (75) that the sequence w ε :=

∑ i∈Iε 1 ε (u ε - 1 |D i ε | ∫ D i ε u ε dx ′ )χ D i ε is bounded in L 2 (0, L; H 1 m (D))
and one can assume possibly by extracting a subsequence that it converges weakly in that space to some w. We then prove that the last limit K arising in (73) is equal to ∇ ′ y w. The function defined by u 0 (x, y) := u 00 (x, y) -v(x) in Ω × Y satisfies the convergence (60). Taking in (59) a test function in the form ϕ(x) = ū0 (x, x ′ ε ) + v(x) + w(x, x ′ ε ) with regular ū0 , w, v and passing to the limit, we obtain the limit equation(64) by a density argument. Choosing ū0 = v = 0 and w = w in that equation, we conclude that w = 0 completing the proof of Theorem 5.1.

Remark 1. In order to emphasize the nonlocal effect at the limit, we have assumed in Theorem 4.2 and in Theorem 5.2 that the function f does not depend on the variable y. One can handle the general case of a source term f (y, x) depending also on the variable y; in this case, setting f (x) := ∫ Y f (y, x) dy, one can check easily that the solution u 0 (x, y) of (60) takes the form u 0 (x, y) = f (x)û(y)+ û0 (x, y) where û is the solution of (65) and where û0 (x, y) is the solution of (60) but with the right hand side given by f (y, x)-∫ Y f (y, x) dy. This is due to the simple remark that any right hand side f (y, x) of (60) may be written as f (y, x) = f (y, x) -∫ Y f (y, x) dy + ∫ Y f (y, x) dy so that the solution u 0 is the superposition of the two solutions. Of course the same remark holds true in the problem of reduction of dimension. In both cases, the nonlocal effect is due to the term f (x)û(y).

  inequality sup x∈Ω |ϕ -ϕ ε | ≤ Cε pointed out above and due to the boundedness in H 1 (Ω) of the sequences u ε and w ε , the first integral in the right hand side of (27) clearly tends to zero, while the second one converges to ∫ Ω ∇u∇ϕ dx since one can assume the strong convergence to 1 in L 2 (Ω) of the sequence w ε by the Rellich Theorem. Hence, it only remains to compute the limit of the third integral.