Durham Research Online

Transcription

1 Durham Research Online Deposited in DRO: 23 April 2018 Version of attached le: Accepted Version Peer-review status of attached le: Peer-reviewed Citation for published item: Cooper, S. (2018 'Quasi-periodic two-scale homogenisation and eective spatial dispersion in high-contrast media.', Calculus of variations and partial dierential equations., 57 (3. p. 76. Further information on publisher's website: Publisher's copyright statement: c The Author(s 2018 Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License ( which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s and the source, provide a link to the Creative Commons license, and indicate if changes were made. Additional information: Use policy The full-text may be used and/or reproduced, and given to third parties in any format or medium, without prior permission or charge, for personal research or study, educational, or not-for-prot purposes provided that: a full bibliographic reference is made to the original source a link is made to the metadata record in DRO the full-text is not changed in any way The full-text must not be sold in any format or medium without the formal permission of the copyright holders. Please consult the full DRO policy for further details. Durham University Library, Stockton Road, Durham DH1 3LY, United Kingdom Tel : +44 ( Fax : +44 (

2 Quasi-periodic two-scale homogenisation and effective spatial dispersion in high-contrast media Shane Cooper April 23, 2018 Abstract The convergence of spectra via two-scale convergence for double-porosity models is well known. A crucial assumption in these works is that the stiff component of the body forms a connected set. We show that under a relaxation of this assumption the (periodic two-scale limit of the operator is insufficient to capture the full asymptotic spectral properties of high-contrast periodic media. Asymptotically, waves of all periods (or quasi-momenta are shown to persist and an appropriate extension of the notion of two-scale convergence is introduced. As a result, homogenised limit equations with none trivial quasimomentum dependence are found as resolvent limits of the original operator family. This results in asymptotic spectral behaviour with a rich dependence on quasimomenta. 1 Introduction The model problem to study time-harmonic waves, with frequency ω, in media with microstructure is div ( a ε ( x ε u = ω 2 u in where the wave u represents the information being propagated, such as pressure in acoustics, deformation in elasticity or electromagnetic fields in electromagnetism 1. The microstructured nature of the media is characterised by periodic coefficients a 2 ε : { a1ε (y, y Q a ε (y = 1, a 0ε (y, y Q 0, where a 0ε, a 1ε are (the square-root of the wave speeds of the individual constitutive material components, see figure 1. The parameter ε represents the ratio between the size of the microstructure and the observable length scale, and is typically taken to be small. From the point of view of applications, it is important to study the asymptotic behaviour of these waves in the limit of vanishing ε. A classical approximation, provided by the homogenisation theorem 3, states that for fixed frequency ω the microstructured media can be approximated by an effective homogeneous media whose wave speed a hom is constant and determined directly from the local periodic behaviour of the problem. The intuition behind why the homogenisation theorem holds is that the wavelength of u is long with respect to the microstructure: variations in u appear over much longer distances than the media s period. Mathematically, this is ensured by assuming that a ε are taken to be uniformly bounded and elliptic with respect to ε, for example a 1ε = a 1, a 0ε = a 0, for bounded elliptic a 0, a 1. 1 In elasticity and electromagnetism the wave equation describes certain polarised waves: e.g. Shear polarised wave in elasticity or Transverse Electric and Transverse Magnetic polarised waves for the Maxwell system. 2 The implied non-trivial dependence of a on ε is deliberate and, as we shall see, important. 3 Also called the long-wavelength or quasi-static approximation depending on the community. 1

3 Figure 1: A typical schematic of some three-dimensional composite media with period microstructure. It has been known for some time now that interesting effects appear when the above elliptic conditions are not uniform. This happens for example in so-called high-contrast media. In the context of waves, high-contrast media of particular interest are the so-called double porosity models which admit the critical scaling: a 1ε = a 1, a 0ε = ε 2 a 0. Physically, this critical scaling corresponds to the wavelength of u remaining long within the media Q 1 but in media Q 0 the wavelength is at the resonant scale, i.e. of the same order as the size of the microstructure. Thus violating the underlying intuition for the long-wavelength approximation. The mathematical analysis of high-contrast problems has given rise to rigorous descriptions of various scale-interaction phenomena such as memory effects and other non-local effects (e.g. [5, 7 9, 14, 17, 28]. Within the context of wave propagation, an important feature of high-contrast problems is that they contain spectral gaps (cf. [18, 30, 31]: frequencies at which no wave can propagate through the underlying medium. Such gaps are important from the point of view of wave-guiding applications such as photonic crystal fibres. An important initial work in the study of the spectrum of high-contrast elliptic operators was undertaken by V. V. Zhikov [30, 31]. Therein, the homogenisation theory for double porosity-type problems was developed within the framework of the so-called two-scale convergence of G. Nguetseng-G. Allaire [1, 27]. Using this theory, Zhikov derived two-scale limit spectral equations that contain a non-trivial coupling between microand macro-scales. Such a coupling leads to an eigenvalue problem with a highly non-linear spectral dependence, described by a function β. The convergence of spectra (in the appropriate sense was proved and, by doing so, demonstrates that this β function provides an explicit description of the asymptotic structure of the spectrum. Such an explicit description of the limit spectral behaviour via two-scale homogenisation has made way for mathematical studies of high-contrast media as wave-guides: in [20] using multi-scale asymptotics and supplemented with analysis based on two-scale convergence in [10]. Moreover, the Zhikov β function was later independently discovered by G. Bouchitté and D. Felbacq [6] in the specific context of TM-polarised electromagnetic wave propagation in a dilute dielectric two-dimensional photonic crystal fibre; therein the authors made the interesting interpretation of the β-function playing the role of effective negative magnetism. Later, in the context of elasticity, a matrix analogue of the β function is derived and plays the role of frequency-dependent effective density [3, 4, 32]. Such works demonstrate that the unusual phenomena observed in high-contrast media can be described by non-standard constitutive laws provided via two-scale homogenisation. The idea that high-contrast media can result in the appearance of non-standard constitutive laws and give rise to composite media with complex wave phenomena near micro-resonances has prompted a recent energetic pursuit of such laws in the contexts of elasticity (e.g. [24] and electromagnetism (e.g. [12, 22]. 2

4 Applications can be found in areas such as cloaking (e.g. [25, 26]. It was shown in the work of V. P. Smyshlyaev [29], building on related ideas in [14], that the two-scale homogenised limit of various anisotropic elastic media contain not only the temporal non-locality (as described by the Zhikov β function but also exhibit spatial non-locality. The presence of which leads to the phenomena of directional localisation: the number of admissible propagating wave modes depends not only on the frequency but on the direction of propagation. Such a feature is important for cloaking applications. These motivating works have led to recent systematic study containing rigorous asymptotic and spectral analysis of general mathematical constructions containing high-contrasts [21]. Analysis based on the work [21] has led to the demonstration that the two-scale convergence is insufficient to fully study the spectrum of general high-contrast problems, see [15, Chapter 5], [11]. The reason for this inconsistency is due to the presence of quasi-periodic microoscillations that persist at leading-order in general high-contrast media. In this work we appropriately develop homogenisation theory to study quasi-periodic micro-oscillations. This is achieved by extending the two-scale convergence framework to admit such oscillations. We explain the (lack of the role these micro-oscillations in the numerous previous works on high-contrast problems. Then, we apply this theory in the spectral analysis of a novel class of high-contrast media. In particular, we shall show that by relaxing the geometric assumptions in the double-porosity model leads to multi-scale homogenised models that contain a new feature: the effective wave speed depends on the quasi-momenta in a highly discontinuous fashion. Specifically, the non-standard constitutive equations for such high-contrast media exhibit spatial dispersion. The presence of this novel feature is related to the contribution of the quasi-periodic waves on the microscale. Notations We end the introduction with some words on the notation used in this article. Vectors and vector-valued functions are represented by lower-case boldface symbols with the exception of the co-ordinate points. {e 1, e 2, e 3 } denotes the Euclidean basis in R 3. For a vector u R 3, we denote by u i its component with respect to e i, and write u = (u 1, u 2, u 3 = 3 u i e i. Points in R 3 will be denoted by the symbol x and points in the unit cell := [0, 1 3 will be denoted by y. The notation i will be used to denote partial differentiation with respect to the i-th coordinate variable, and we shall replace the suffix i with x i or y i when we wish to emphasis the macroscopic or microscopic variable. Similarly, the notion div x, div y, x, and y, are used for the divergence or gradient of a function in terms of x or y. Throughout is a domain in R d, d 1, := (0, 1 d and [0, 2π d. All of the functions, even if real-valued, are considered to take values in the complex field. The space C# ( denotes the usual space of smooth -periodic functions. Whereas, C ( shall denote the space of smooth functions ϕ(y whose functions and derivatives are -quasi-periodic with respect to y: ϕ(y + e j = exp(i j ϕ(y for each y and each Euclidean basis vector e j, j = 1,..., d. Equivalently, i=1 C ( = {φ φ = e i y ψ, ψ C # ( }. Note that C0 ( = C# ( and use the latter to avoid confusion with the notation for the space compactly supported smooth functions. The Sobolev space H# 1 ( is the usual Sobolev space of H1 -periodic functions. Whereas H 1 ( is defined as the closure of C ( with respect to the H1 norm, or equivalently as H 1 ( := {e i y u # u # H 1 #( }. (1.1 Also, we note H0( 1 = H# 1 ( and, in this situation, we use the latter to avoid confusion with the Sobolev space of zero trace H 1 functions. 3

5 For subsets {S ε } ε and S of R d we say that S ε converges to S in the Hausdorff sense if the following conditions hold: 1. For every λ ε S ε such that λ ε converges to some λ 0, then λ 0 S. 2. For every λ 0 S there existsλ ε S ε such that lim ε λ ε = λ 0. We shall use the notation lim ε S ε = S when a sequence of sets S ε Hausdorff converges to S. The Einstein summation convention will not be used in this article, that is we do not sum with respect to repeated indices. 2 Quasi-periodic two-scale convergence In this section we introduce an appropriate notion of convergence that will account for the presence of microscopic oscillations that are quasi-periodic in nature. This convergence will turn out to be a natural extension of the standard (periodic two-scale convergence introduced by G. Nguetseng [27]-G. Allaire [1]. In particular, we aim to use this extended notion of two-scale convergence to study the spectral properties of operator families in homogenisation theory in a similar vein to that first introduced by V. V. Zhikov [30, 31]. 2.1 Motivation We shall motivate the notion of quasi-periodic two-scale convergence here. This motivation is based on the principle goal of characterising the spectral asymptotics of high-contrast elliptic operators. Let ε be a sequence of positive real numbers with limit zero. Consider the differential operator A ε : D(A ε L 2 (R d L 2 (R d whose action is given by A ε u := div ( a ε ( ε u and domain D(A ε consists of u for which A ε u L 2 (R d. Here a ε are -periodic measurable functions, that are bounded and elliptic: ν 1, ν 2 > 0 such that ν 1 I a ε ν 2 I. In this article, we focus on a ε that are uniformly bounded, i.e. ν 2 is independent of ε, but a ε may degenerate in the sense that ν 1 = ν 1 (ε with lim ε ν 1 (ε 0. We are interested in analysing the structure of the spectrum σ(a ε of A ε in the limit of ε. The strategy of the study is to establish the existence of some operator A 0 such that σ(a ε Hausdorff converges to σ(a 0 ; i.e. the following conditions hold: 1. For every λ ε σ(a ε such that λ ε converges to some λ 0 we deduce that λ 0 σ(a For every λ 0 σ(a 0 we find λ ε σ(a ε such that lim ε λ ε = λ 0. A crucial question is how to determine the operator A 0. For example, in classical and semi-classical highcontrast problems, A 0 turns out the be the strong two-scale resolvent homogenised limit of A ε, cf. [16,30,31]. To develop some intuition on what to expect in the general case, let us recall an important result from the spectral theory of elliptic operators with ε -periodic coefficients: the Floquet-Bloch decomposition (see for example [23] for more details. This result states that the following characterisation of σ(a ε holds: σ(a ε = σ(a ε (Θ [ Θ 0, 2π d ε where A ε (Θ : D(A ε (Θ L 2 (ε L 2 (ε, describe a family of densely defined self-adjoint operators with compact resolvent given by the action that A ε (Θu = f L 2 ( if u H# 1 (ε solves div ( a ε ( x ε eiθ x u = e iθ x f(x, x ε. 4

6 Taking the above into consideration we see that λ ε σ(a ε if, and only if, there exists Θ ε [ 0, 2π d ε and non-trivial u ε H# 1 (ε such that div ( a ε ( x ε eiθ x u ε = λε e iθ x u ε (x, x ε. By a change of variables y = x/ε and = εθ, we see that w ε (y := e i y u ε (εy solves div ( ε 2 a ε (y w ε = λε w ε (y, y, (2.1 and w ε belongs to the space of H 1 ( functions that satisfy the condition w ε (y + z = e i z w ε (y, y, z Z d, for some [0, 2π d. This condition is typically referred to as the Bloch or quasi-periodic condition and is known as the quasi-momentum. Note that = 0 is the usual periodicity condition. The Sobolev space of H 1 ( -quasi-periodic functions coincides with H 1 ( which, we recall, is be defined as H 1 ( := {e i y u # u # H 1 #( }. The general principle to observe here is that if we wish to study the asymptotic behaviour of the spectrum σ(a ε we need to keep track of the eigenfunctions that are -quasi-periodic on micro-scale y := x/ε, for all [0, 2π d. The notion of quasi-periodic two-scale convergence, introduced in Section 2.2 below, performs such a task. We note here that in the case of the whole space, as discussed above, one need not refer to a notion of quasi-periodic two-scale convergence to study the asymptotics of the spectrum; one may study the normresolvent limits of the operators A ε ( to study spectral asymptotics, cf. [18] where the point-wise (in limits or [13] where the uniform limits were considered in the double-porosity setting. That being said, for boundary-value problems, the Bloch decomposition does not hold; nevertheless, the whole space (or Bloch spectrum is expected to contribute asymptotically to the bounded domain spectrum and the precluding discussion is still relevant. It is this setting that the quasi-periodic two-scale convergence will be particularly useful. Finally, we comment that the above discussion leads to the natural question: why in previous cases considered was it sufficient to consider the standard (periodic two-scale limit of A ε to ensure spectral convergence? Or put another way, when in the asymptotic limit of ε do we need consider all quasi-periodicity and not just = 0. This shall be explained in Section Definition and basic properties This section is dedicated to the introduction of the notion quasi-periodic two-scale convergence and an exposition of results that are appropriate to homogenisation theory. Recall C# ( denotes the usual space of smooth -periodic functions and C ( = {φ φ = e i y ψ, ψ C # ( }. The following mean-value property will be important: For every ϕ C ( and every φ C 0 ( the following convergence ( lim φ(xϕ x 2 ε 0 ε dx = φ(xϕ(y 2, (2.2 holds. This fact follows by noting that the assertion holds for elements in C# (, see for example [1, Lemma 1.3], and observing that multiplication by exp( i y defines an isomorphism between C ( and C # ( that preserves absolute value. Indeed, ϕ belongs to C ( if, and only if, exp( i yϕ belongs to C # ( and ϕ = exp( i yϕ in. 5

7 We remark here that because C ( is isomorphic to C # ( = C 0 (, with isomorphism exp ( i y, then the results presented in this section 4 are immediately established for each [0, 2π 3 if proved for = 0. We shall demonstrate this with the first result of the section and omit the remaining proofs which follow in a similar manner. Definition 1. Let u ε L 2 ( be a bounded sequence and u L 2 (. Then, we say u ε (weakly 2 -quasi-periodic two-scale converges to u, denoted by u ε u, if the following convergence lim u ε (xφ(xϕ ( x ε 0 ε dx = u(x, yφ(xϕ(y dydx, φ C0 (, ϕ C ( (2.3 holds. Remark 1. Notice that for = 0, this is the standard notion of two-scale convergence. The next important result states that bounded sequences in L 2 ( are relatively compact with respect to quasi-periodic two-scale convergence. Proposition 1. If u ε is bounded in L 2 ( then, up to a subsequence, u ε weakly -quasi-periodic two-scale converges to some u L 2 (. Proof. The result has been established previously for the case = 0, see for example [1, 27, 30]. Let us consider 0. Note that the function ũ ε = exp ( i y u ε is bounded in L 2 and therefore by the assertion for = 0, up to a discarded subsequence, ũ ε (0-quasi-periodically two-scale converges to some ũ L 2 ( Q. Now, the result follows from this fact and noting that for fixed ϕ C ( one has u ε (xφ(xϕ( x ε dx = exp ( i y ũ(x, yφ(xϕ(y dydx = ũ ε (xφ(xexp ( i x ε ϕ( x ε dx, Q ũ(x, yφ(xexp ( i y ϕ(y dydx, and that exp ( i y is a smooth periodic function. Hence u ε 2 exp ( i y ũ. An important result from the point of view of homogenisation theory is that the test functions ϕ in (2.3 can be taken to be quasi-periodic elements of L 2 (, i.e. the following result holds. Proposition 2. If u ε L 2 ( -quasi-periodic two-scale converges to u L 2 (, then the following convergence lim u ε (xφ(xψ ( x ε 0 ε dx = u(x, yφ(xψ(y dydx holds for all φ C 0 (, and for all ψ L 2 ( such that ψ(y + e j = exp(i j ψ(y for almost every y and j = 1,..., d. Remark 2. If is a bounded domain, as in this article, then additionally the test functions φ can be taken to be elements of C(. The following results are of interest. Proposition For u ε L 2 ( -quasi-periodic two-scale converging to u L 2 ( one has that exp( i x ε u ε(x exp( i yu(x, y dy weakly in L 2 (. 4 The results in this section can be established by first principles making no reference to such an isomorphism. 6

8 2. For u ε L 2 ( -quasi-periodic two-scale converging to u L 2 ( then lim inf u ε (x 2 dx u(x, y 2 dydx. ε 0 A result of particular interest in high-contrast homogenisation problems is the following. Proposition 4. Let u ε H 1 ( satisfy sup u ε L 2 ( <, ε sup ε u ε L 2 ( <. ε Then, there exists u L 2 (; H 1 ( such that, up to a subsequence, the following convergences hold: u ε 2 u, ε u ε 2 y u. Recall here that H 1 ( is given by (1.1. Proof. Let φ and Ψ denote respectively fixed arbitrary elements of C 0 ( and C ( ; Cd. By Proposition 1, there exists u L 2 ( and χ L 2 ( ; C d such that, up to a discarded subsequence, the following convergences hold: u ε 2 u, ε u ε 2 χ. (2.4 Note that, since u ε is bounded in L 2 (, then εu ε strongly converges to zero in L 2 ( and from Proposition 3 part 2. we conclude that εu ε 2 0. (2.5 Let us prove that u L 2 (; H 1 (Q. Using the convergences (2.4 and (2.5 we pass to the limit in the identity ε u ε (x φ(xψ( x ε dx = u ε (xεdiv ( φ(xψ( x ε dx = u ε (xε x φ(x Ψ( x ε dx u ε (xφ(xdiv y Ψ( x ε dx, to deduce that χ(x, y φ(xψ(y dydx = u(x, yφ(xdiv y Ψ(y dydx. Therefore, for almost every x, the functions χ(x, and u(x, are related by the identity χ(x, y Ψ(y dy = u(x, ydiv y Ψ(y dy, Ψ C ( ; C d. ( It is clear that C 0 (0, 1 d C ( and so u H1 ( with y u = χ. It remains to show u belongs to H 1 (. This follows from noting that after performing integration by parts in the above identity we arrive at u(x, yψ(y ν ds(y = 0, Ψ C ( ; C d. Setting Ψ = e i y Ψ # above, for arbitrary smooth -periodic Ψ #, demonstrates that u # (x, := e i y u(x, is an element of H 1 ( that satisfies periodic boundary conditions with respect to y. That is, u # (x, H 1 # ( and so (see definition (1.1 u(x, H1 (. 7

9 We end this section with a result that is illuminating when it comes to studying the convergence of spectra for parameter-dependent operator families. It readily provides a one-sided justification for the Hausdorff convergence of the high-contrast spectra to the spectrum associated to quasi-periodic two-scale limits. This result is based on the following definition, which extends the notion of strong resolvent two-scale convergence first introduced by V. V. Zhikov in [30, 31]. Definition 2. Fix [0, 2π d, and let A ε and A be non-negative self-adjoint operators respectively defined in L 2 ( and H a closed subset of L 2 (. We say that A ε strong resolvent -quasi-periodic two-scale converges to A if for every f ε (x L 2 ( that -quasi-periodic two-scale converges to f(x, y L 2 (, the following convergence u ε = (A ε + I 1 f ε 2 u = (A + I 1 P f, as ε 0 holds. Here, P is the orthogonal projection onto H in L 2 (. Here, we state an important consequence of such resolvent convergence. The proof, omitted here, follows standard spectral theoretic arguments, see for example [30]. Proposition 5. If A ε strong resolvent -quasi-periodic two-scale converges to A then the spectrum σ(a of A is related to the spectrum σ(a ε of A ε in the following sense: For every λ σ(a there exists λ ε σ(a ε such that λ ε converges to λ as ε tends to zero. 2.3 On the relevance of quasi-periodic two-scale convergence in spectral asymptotics Proposition 5 informs us that, in principle, one should consider all strong quasi-periodic two-scale limits of an operator A ε to fully characterise its limit spectrum (in the Hausdorff sense. Yet, clearly this is not always the case: such a notion of convergence has not appeared previously, nor was it needed, to study the spectral asymptotics of classical and particular double-porosity operators. The reason for this shall be elucidated here. Moreover, at the end of this section we shall argue when quasi-periodic convergence is necessary via a model problem that we later study in detail in this article Classical homogenisation Consider the resolvent problem: For fixed f L 2 ( find u ε H 1 0 ( such that div ( a( x ε u ε + uε = f, (2.6 where the symmetric matrix-valued function a is -periodic, elliptic and bounded: ν > 0 such that The following homogenisation theorem is classical. ν ξ 2 a(yξ ξ ν 1 ξ 2, ξ C d, a.e. y. Theorem 1 (Classical homogenisation theorem. Let ε be a sequence with limit 0, and f ε L 2 ( a sequence such that f ε weakly converges in L 2 ( to some f 0 as ε tends to zero. Then u ε H 1 0 ( the solution to (2.6, for f = f ε, converges weakly in H 1 0 ( (and strongly in L 2 ( to u 0 H 1 0 ( the solution to div ( a hom u 0 + u0 = f 0. Here a hom is the constant symmetric homogenised matrix determined by a: a hom ξ ξ := a ( N + ξ ( N + ξ, ξ R d. min N H 1 # ( 8

10 It is well-known that the homogenisation theorem implies the Hausdorff convergence of spectra (cf. [2, Section 2]: lim ε σ(a ε = σ(a hom. Let us study the quasi-periodic two-scale limits of u ε. Proposition 6. Fix (0, 2π d and consider f ε L 2 2 ( such that f ε f 0. Then, u ε H0 1 ( -quasiperiodic two-scale converges to zero; that is u ε 2 0. Remark Proposition 6 informs us that for the classical resolvent problem (2.6, the non-zero quasiperiodic micro-oscillations at leading order do not contribute to the spectral asymptotics. So one need only study the = 0 quasi-periodic oscillations, i.e. the standard two-scale limit. It is well-known that the (periodic two-scale limit coincides with the classical limit provided by Theorem 1, see for example [1, 30]. 2. The part of the spectrum corresponding the -quasi-periodic micro-oscillations, for 0, actually resides in an ε 2 neighbourhood of infinity; this can be formally seen from the considerations of Section 2.1: for a ε independent of ε, the eigenvalues λ ε in (2.1 are clearly of the order ε 2. To study such high-frequency spectrum one can consider the re-scaled operator ε 2 A ε, that is consider coefficients of the form a ε = ε 2 a. The precise study of such high-frequency spectra was performed in [2] for a broader class of moderately contrasting locally periodic coefficients. Therein, the authors provide a rigorous description of the high-frequency spectral asymptotics in terms of non-trivial quasi-momenta. This was done by introducing an appropriate notion of Bloch wave homogenisation. For the reduced setting of (globally periodic coefficients, the Bloch-wave operator-limits determined therin can readily be shown to be equivalent to the -quasi-periodic two-scale limits. Proof of Proposition 6. The sequence f ε weakly converges, cf. Proposition 3 part 1., and so is bounded. Multiplying (2.6 (for f = f ε and integrating over, and using the ellipticity of a, produces the a-priori bound u ε 2 L 2 ( + ν u ε 2 L 2 ( f ε 2 L 2 ( C <. Applying Proposition 4, we deduce that there exists u L 2 (; H 1 ( such that, up to a subsequence, the following convergences hold: u ε 2 u, ε u ε 2 y u. Let us show u = 0: u ε is a bounded sequence and so ε u ε strongly converges to zero in L 2. Therefore, by Proposition 3 part 2., we deduce y u = 0. As is connected it follows that u is constant. Yet u ε H 1 ( and there are no non-trivial constant -quasi-periodic functions for 0, see (1.1. Hence, u = Double-porosity model Consider the resolvent problem: For fixed f L 2 ( find u ε H 1 0 ( such that Here div ( a ε ( x ε u ε + uε = f. (2.7 a ε (y = { a1 (y, y Q 1, ε 2 a 0 (y, y Q 0, where Q 0 is a smooth compactly contained subset of such that, for Q 1 := \Q 0, the periodic extension F 1 := z Z d (Q 1 + z forms a connected set in R d. The functions a i, i = 0, 1 are taken to be real-valued, elliptic and bounded on Q i. The following homogenisation theorem is established in [30, Theorem 5.1]. 9

11 Theorem 2. Suppose f = f ε in the right-hand side of (2.7 two-scale converges to some f 0, that is f ε 2 f 0 for = 0. Then, the sequence of solutions u ε two-scale converges to u 0 (x, y = u(x + v(x, y, where (u, v belongs to V 0 = H 1 0 ( L 2 (; H 1 0 (Q 0 and uniquely solves a hom dp x u(x x φ(x dx + a 0 (y y v(x, y y ϕ(x, y dydx + (u(x + v(x, y (φ(x + ϕ(x, y dydx = f 0 (x, y (φ(x + ϕ(x, y dydx, φ H 1 0 (, ϕ L 2 (; H 1 0 (Q 0. (2.8 Here, a hom dp by a 1 : is the constant symmetric and positive homogenised matrix for perforated domains determined a hom ( ( dp ξ ξ := min a 1 N + ξ N + ξ, ξ R d. N H# 1 (Q1 Q 1 This result informs us that A ε strongly two-scale converges to the operator A 0, defined in L 2 (, associated to the above two-scale limit resolvent problem. Therefore, by Proposition 5 for = 0, the lowersemicontinuity of the spectral convergence is ensured. In fact, Zhikov proved in [30, Theorem 8.1], under the condition that F 1 is connected in R d, the stronger result lim ε σ(a ε = σ(a 0. Let us determine the strong resolvent quasi-periodic two-scale limits of A ε. Proposition 7. Fix (0, 2π d 2. Suppose f = f ε in the right-hand side of (2.7 such that f ε f 0 to some f 0 L 2 (. Then, the sequence of solutions u ε -quasi-periodically two-scale converges to v 0 (x, y L 2 (; H0 1 (Q 0 the solution to a 0 (y y v 0 (x, y y ϕ(x, y dydx + v(x, y ϕ(x, y dydx = f 0 (x, y ϕ(x, y dydx, ϕ L 2 (; H 1 0 (Q 0. (2.9 Remark Note that A = A is independent of for 0, and its spectrum is the point spectrum given by the operator whose action is u div(a 0 u with domain {u H0 1 (Q 0 div(a 0 u L 2 (Q 0 }. 2. It is easy to see A A 0 (by noting that setting u = φ = 0 in (2.8 gives (2.9 and so σ(a = σ(a σ(a 0. The set 0 σ(a, 0 is the limit spectrum arriving from quasi-periodic micro-oscillations. 3. The restriction of the limit spectrum σ(a 0 to 0 σ(a is achieved by considering the purely macroscopic component u(x (of eigenfunctions to be zero. For this reason, we coin this spectrum to be pure Bloch spectrum. In the simplified setting of double-porosity the pure Bloch spectrum is point spectrum (due to the fact A = A is independent of for 0. In general, we expect this spectrum to have band-gap structure, 10

12 and the gaps have only contracted to points here due to the geometric constraint that F 1 is connected in R d. This expectation is verified in Section Even though the strong resolvent -quasi-periodic limit of A ε exists, it has trivial dependence on, 0 and more importantly is a restriction of the two-scale limit A 0. Hence, one need only consider A 0, and this explains why in this setting one is to expect Zhikov s result lim ε σ(a ε = σ(a 0. In general, the limit A 0 will not be sufficient to capture the full spectral asymptotics. Proof of Proposition 7. Let us consider a 1 (respect. a 0 to be extended by zero into Q 0 (respect Q 1, and consider ν > 0 to be the constant such that The solution u ε solves a 1 + a 0 ν. (a 1 ( x ε + ε2 a 0 ( x ε u ε φ + u ε φ = f ε φ, φ H 1 0 (. (2.10 Setting φ = u ε in the above variational problem and using the fact that a 1 + a 0 ν, we deduce the a-priori bound (for ε 1 u ε 2 L 2 ( + ν ε u ε 2 L 2 ( f ε 2 L 2 ( C <. Additionally, we have the bound Indeed, a 1 0 and a 1 ( x ε u ε u ε a 1 ( x ε u ε 2 L 2 ( f ε 2 L 2 ( C <. (a 1 ( x ε + ε2 a 0 ( x ε u ε u ε + By Proposition 4 it follows that, up to a discarded subsequence, for some v 0 L 2 (; H 1 (Q. Let us show that v 0 L 2 (; H 1 0 (Q 0. By Proposition 2 it follows that u ε 2 = f ε u ε. u ε 2 v 0, and ε u ε 2 y v 0, (2.11 a1 ( x ε ε u ε 2 a 1 (y y v 0. Yet a 1 ( x ε ε u ε strongly converges to zero in L 2 (. Therefore a 1 (y y v 0 = 0, which is equivalent to y v 0 = 0 on Q 1 (recall a 1 is positive on Q 1 and zero on Q 0. As Q 1 is connected then v 0 is constant in Q 1. Now, since the periodic extension F 1 = z Z d(q 1 + z forms a connected set then v 0 is constant in F 1. Yet, v 0 H 1 ( for 0 and consequently this constant is zero, cf (1.1. It remains to prove v 0 solves (2.9. This can easily be deduced by passing the the -quasi-periodic limit in (2.10 for test functions φ(x = ψ(xϕ( x ε, ψ H1 0 (, ϕ H0 1 (Q 0 and using convergences ( An example with non-trivial quasi-periodic limits Let us provide an example which demonstrates that in general the family A, of strong resolvent -quasiperiodic limits to A ε, are not restrictions of A 0. Suppose, we consider (2.7 for coefficients a ε = a 1 + ε 2 a 0, 11

13 and a i are real-valued -functions such that a 1 0 and a 1 + a 0 ν > 0. Let f ε be a bounded sequence and u ε solve (2.7 for f = f ε. Arguing as in the proof of Proposition 7, we see that u ε to solution to (2.7 will satisfy the a-priori bounds u ε 2 L 2 ( + ν ε u ε 2 L 2 ( + a 1 ( x ε u ε 2 L 2 ( f ε 2 L 2 ( C <. In particular, Proposition 4 informs us that up to a subsequence u ε 2 u 0, and ε u ε 2 y u 0, for some u 0 L 2 (; H 1 (Q. Moreover, by an application of Proposition 2 we deduce that a1 y u 0 = 0. Denoting by V the closed linear subspace of H 1 ( given by5 V = {v H 1 ( a 1 y v = 0}. Suppose we show an example where V is not a subset of V 0 for some non-trivial V, 0. Then, for such examples we should not expect that A 0 is an extension of A, nor should we expect σ(a σ(a 0. Let us provide such an example. The conjectures (stated immediately above based on this example will be proved rigorously in the remainder of the article. Suppose Q 1 is the cylinderical domain Q 1 := [0, 1 [ 1 4, 3 4 ]2, and a 1 (y := { 1, y Q1, 0, y Q 0. Notice that F 1 := z Z 3(Q 1 + z the periodic extension of Q 1 into R d consists of infinitely many mutually disjoint cylinders C l := R [ l 1, l 2] 2, l Z 2. That is, the assumptions of Subsection 2.3.2, and in particular [30, 31], do not hold. Now v V if, and only if, v H 1 ( with v = c for some c C on Q 1. This implies, cf (1.1, that v(1, y 2, y 3 = e i1 v(0, y 2, y 3, (y 1, y 2 (0, 1 2, j {1, 2, 3}, (in the sense of trace. Then, for (y 2, y 3 ( 1 4, we arrive at the condition c = e i1 c. Therefore, if 1 0, then the above condition only holds if c = 0. That is v must necessarily be zero on Q 1. On the other hand, if 1 = 0 then any H 1 ( function that is constant on Q 1 belongs to V. In particular we see that V does not belong to V 0 for all (0, 2π 3. Remark 5. Note that if Q 1 contains a connected subset which joins two opposite faces of the square then the space V non-trivially depends on. Consequently, non-trivial limit Bloch spectrum is expected for aligned orthogonally to these faces. The remainder of the article is dedicated to determining the strong quasi-periodic two-scale limits of A ε for such fibre-like inclusions. Moreover, we demonstrate that indeed A are not restrictions of A 0 and that σ(a form non-trivial subsets of lim ε σ(a ε. 5 The space V 0 was first introduced in [29] and coined the space of microscopic oscillations. We have appropriately extended this notion to -quasi-periodic oscillations here, 0. 12

14 3 Problem formulation and Homogenisation In this article we are concerned with the asymptotic analysis of the resolvent problem { Find uε H 1 0 ( such that div ( a ε ( x ε u ε + uε = f ε in (3.1 where ε < 1 is a small parameter, is a smooth open bounded star-shaped domain 6 and f ε L 2 ( known. The coefficient a ε is given by a ε (y = { a1 (y, y Q 1, ε 2 a 0 (y, y Q 0, 0 a i, a 1 i L (Q i, a i = 0 on Q 1 i, i = 0, 1, (3.2 and the regions Q 0 and Q 1 are described as follows (cf. Figure 2. Geometric assumptions. For j = 1, 2, 3 we consider smooth domains S j compactly contained in (0, 1 2 that have mutually disjoint closures. We denote by C j the cylinder aligned to the j-th co-ordinate axis with cross-section S j, i.e. C 1 := {y (0, 1 3 y (0, 1 S1 }, C 2 := {y (0, 1 3 y = (z2, z 3, z 1, z (0, 1 S 2 } and C 3 := {y (0, 1 3 y = (z 3, z 1, z 2, z (0, 1 S 3 }. Then, for a given non-empty subset I of {1, 2, 3}, we consider Q 1 = i I C i. We denote by Γ i = C i \ Q and Γ = i I Γ i = Q 1 \ Q. Figure 2: a An example stiff component Q 1 consisting of one cylinder C 1 extending in the x 1 direction, i.e. I = {1}. b The stiff component Q 1 consists of two disjoint cylinders extending in the co-ordinate directions x 1, and x 3, i.e. I = {1, 3}. c Q 1 consists of mutually disjoint cylinders extending in all co-ordinate directions, i.e. I = {1, 2, 3}. Under this geometric assumption we determine for each [0, 2π 3 the strong resolvent -quasi-periodic two-scale limit, cf. Section 2, of the self-adjoint operator A ε associated to resolvent problem (3.1. That is, for a fixed [0, 2π 3 and a given bounded sequence f ε L 2 2 ( such that f ε f, i.e. lim ε 0 f ε (xφ(xϕ ( x ε dx = f(x, yφ(xϕ(y dydx, φ C 0 (, ϕ C ( we aim to determine the -quasi-periodic two-scale limit behaviour of the solution u ε H0 1 ( to ( a1 ( x ε + ε2 a 0 ( x ε u ε (x φ(x dx + u ε (xφ(x dx = f ε (xφ(x dx, φ C0 (. (3.3 As f ε is bounded in L 2 (, upon setting φ = u ε in (3.3 we deduce that the sequences a 1 ( ε u ε L2 (;C 3, ε u ε L2 (;C 3, and u ε L2 (, (3.4 6 All the results and proofs follow through in an identical manner for the case where is the whole space. 13

15 are bounded. Let us describe the -quasi-periodic two-scale limit, referring to Section 4 for the details. The limit of u ε (x will be a function u(x, y, of two variables x, y Q, that is -quasi-periodic with respect to the second variable y, cf Proposition 4. Furthermore, due to the fact that in each cylinder C i, i I, the gradient of u ε is bounded, the limit u necessarily belongs to the (Bochner space L 2 (; V where V := {v H 1 (Q v is constant in C i for each i I}. (3.5 It follows from this (see (1.1 that u is non-zero in cylinder C i if and only if the i-th component i of is zero. If i = 0, then we determine that u i is not only non-trivial but it is actually more regular in the x i -th coordinate direction: xi u i L 2 (. More precisely, for I the subset of indexes I {1, 2, 3} given by I := {i I i = 0}, we denote by C the closed subspace of C 3 spanned by {e i } i I, 7 and show that the function u belongs to the set U = { u L 2 (; H 1 (Q u = ui on C i, for some u L 2 (; C with i u i L 2 ( and u i ν i = 0 on }, (3.6 which is clearly a Hilbert space when endowed with the inner product (u, v U := i u i (x i v i (x dx + y u(x, y y v(x, y dydx i I Q 0 + u(x, yv(x, y, dydx. Here, ν is the outer unit normal to. For each fibre C i there corresponds an effective constant material parameter a hom i > 0 given by a hom i = a 1 (y[ yi N (i (y + 1] dy, (3.7 C i where N (i H# 1 i (C i := {u H 1 (C i u is 1-periodic in the variable y i } is the unique solution to the cell problem a 1 (y [ N (i ] (y + e i φ(y dy = 0, φ H 1 #i (C i, C i (3.8 N (i = 0. C i Then, for each [0, 2π 3, the -quasi-periodic two-scale limit problem is formulated as follows: For f L 2 ( Q find u U such that a hom i xi u i (x xi φ i (x dx + a 0 (y y u(x, y y φ(x, y dydx i I Q 0 (3.9 + u(x, yφ(x, y dydx = f(x, yφ(x, y dydx, φ U. As a hom i Q Q are positive numbers and a 1 0 L (Q 0 it follows that the left-hand side of the above problem defines an equivalent inner product on the space U, and consequently the existence and uniqueness of solutions u to (3.9 are ensured by the Riesz representation theorem. Setting φ = 0 on Q 1 in (3.9 gives the equation div y ( a0 (y y u(x, y + u(x, y = f(x, y, x, y Q 0, 7 Note that C is either the whole space, a plane or a line in C 3. Q 14

16 and a subsequent integration by parts in (3.9 leads to the variational formula a hom i xi u i (x xi φ i (x dx + ( a 0 (y y u(x, y n(y dy φ i (xdx i I i I Γ i + C i u i (xφ i (xdx = ( f(x, y dy φ i (xdx, i I i I C i for all φ L 2 (, C, such that i φ i L 2 (. For each fixed j I we set φ j = φ, φ C0 ( and φ i = 0 for i j, above. This leads to the -quasi-periodic two-scale homogenised system of equations. a hom j x 2 j u j (x + T j (u + C j u j = f j (x, x, for each j I ( : div y a0 (y y u(x, y + u(x, y = f(x, y, x, y Q 0, (3.10 u 0 = u j on Γ j, u j ν j = 0 on, Here and u 0 = 0 on Γ i, for i I\I. T j (u(x = a 0 (y y u(x, y n(y ds(y, Γ j f j (x = f(x, y dy, C j for n the outer unit normal of Γ j = C j \ Q. We now state the main result of the article. Theorem 3. Consider f ε L 2 (, f L 2 ( Q such that f ε 2 f, and u ε the solution to (3.3. Then u ε converges, up to some subsequence, in the -quasi-periodic sense to u U the unique solution to (3.9, equivalently (3.10. An immediate consequence of Theorem 3 is that for each [0, 2π 3, the operator A ε strong resolvent -quasi-periodically two-scale converges to the operator A hom associated to problem (3.10, see Definition 2 in Section 2. Consequently, Proposition 5 informs us that the lower semi-continuity of the spectra in the Hausdorff sense is ensured: for every λ σ(a hom λ ε σ(a ε such that lim λ ε = λ. ε 0 [0,2π 3 The structure of the limit spectrum 10. σ(a hom [0,2π 3 is analysed in Section 6 and described in Proposition Remark 6. A seperate issue, not explored here, is the so-called spectral completeness statement, i.e. the question of whether or not the remaining criterion for Hausdorff convergence of spectra is satisfied: does it follow that for every λ ε such that lim λ ε = λ, then λ? ε 0 [0,2π 3 σ(a hom In general this will not be true due to the presence of the boundary, and the fact that Q 0 intersects the boundary. This leads to the expectation that there exists non-trivial spectrum due to surface waves asymptotically localised near the boundary, cf. [2] for analoguous results in the context of classical locally periodic media. For the case of being the torus or the whole space the above assertion is expected to hold and will be explored in future works. 4 Proof of the homogenisation theorem This section is dedicated to the proof of Theorem 3. To do this, we shall develop an appropriate quasi-periodic two-scale variation of a powerful method first introduced in [21] in the context of standard (periodic twoscale convergence, i.e. -quasi-periodic two-scale convergence for = 0. In what follows φ, ϕ and Φ will denote respectively fixed arbitrary elements of C 0 (, C (Q and C (Q; C3. 15

17 4.1 Technical preliminaries The following results will be of importance in the proof of the homogenisation theorem. Lemma 1. Let B be the closure of a smooth domain and let B 1 be a smooth bounded domain such that B B 1 and A = B 1 \B is a connected set. Then every u H 1( (0, 1 A can be extended to (0, 1 B 1 as a function ũ H 1( (0, 1 B 1 such that ũ 2 c u 2, (0,1 B 1 (0,1 A (4.1 ũ 2 c u 2, (0,1 B 1 where c does not depend on u H 1 ((0, 1 A. (0,1 A Proof. Suppose u H 1( (0, 1 A. Then, by Fubini s theorem, for almost every x 1 (0, 1 the function u(x 1, belongs to H 1 (A and let Eu(x 1, be the Sobolev extension of u(x 1, into B 1 given in [19, Lemma 3.2, pg. 88]. In particular, one has Eu(x 1, 2 c u(x 1, 2, B 1 A (4.2 Eu(x 1, 2 c u(x 1, 2, B 1 A where c does not depend on u nor x 1. Here, denotes the gradient vector (0, x2, x3. Consider ũ given by u(x 1, := Eu(x 1,, a.e. x 1 (0, 1. (4.3 Then u(x1, = Eu(x 1, and from assertion (4.2 it follows that ũ 2 c u 2, (0,1 B 1 (0,1 A ũ 2 c u 2. (0,1 B 1 (0,1 A To prove (4.1, it remains to demonstrate that x1 ũ L 2( (0, 1 B 1 and 1 0 B 1 x1 ũ 2 (x 1, x dx dx 1 c For each t R, the difference quotient is given by 1 D t ũ(x 1, x := ũ(x 1 + t, x ũ(x 1, x, t 0 A x1 u 2 (x 1, x dx dx 1. (4.4 where we have extended ũ trivially by zero into R\(0, 1. Notice that D t ũ = ED t u, i.e. the extension into B 1 of the function u(x1+t, u(x1, t = D t u(x 1 + t,, and consequently 1 1 D t ũ 2 (x 1, x dx dx 1 c D t u 2 (x 1, x dx dx 1. 0 B 1 0 A Since D t u converges strongly in L 2( (0, 1 A to x1 u, it follows that D t ũ is a Cauchy sequence ( in L 2( (0, 1 B 1 and this limit can be identified, using the fact that (Dt ũ, φ L 2 = (ũ, D t φ L 2 for φ C0 (0, 1 B1, as x1 ũ. Furthermore, passing to the limit in the above inequality yields (

18 Proposition 8. Fix [0, 2π 3. There exists a constant C = C > 0 such that inf v V u v H 1 (Q C a 1 u L2 (Q, u H 1 (Q. Here V = ker a 1 = {u H 1 (Q a 1 u 0}. Proof. For each i I, let S i be the cross-section of the cylinder C i. Since S i are compactly contained in (0, 1 2 and have mutually disjoint closures then there exists open A i such that S i A i (0, 1 2 and A i are mutually disjoint. Let χ i C0 (A i be smooth cut-off functions that are identity on S i, we extend χ i by zero to (0, 1 2. Now using Lemma 1, let ũ i be the extension of u Ci H 1 (C i to H 1 (D i, where D i is the cylinder whose axis is parallel to y i with cross-section A i. Note that since u is i -quasi-periodic in the variable y i, then the extension will be also, see (4.3. By Lemma 1 it follows that ũ i 2 c D i u i 2. C i (4.5 For i 0, the following Poincaré inequality ũ i 2 i 2 D i holds 8. For i = 0, one has D i ũ i 2 (4.6 ũ i ũ i 2 C D i ũ i 2 D i (4.7 for some C > 0. Here ũ i := 1 D i D i ũ i. Recalling, I = {i I i = 0}, we set ũ = i I χ i(ũ i ũ i + i I\I χ iũ i, here χ i are taken to be constant in the variable y i and as above the complementary directions. It follows that ũ H 1 (Q and u ũ V. Note that, by construction and inequalities (4.6, and (4.7, one has ũ 2 H 1 (Q c ũ i 2 L 2 (D. i Now, the positivity of a 1 on Q 1 and (4.5 imply that the element v := u ũ of V is such that u v 2 H 1 (Q = ũ 2 H 1 (Q C a 1 u 2, Q 1 and the result follows. i I 4.2 Proof of Theorem 3 Consider the sequence u ε of solutions to (3.3, i.e. ( a1 ( x ε + ε2 a 0 ( x ε u ε (x φ(x dx + u ε (xφ(x dx = 2 for f ε f, and recall, cf. (3.4, that sup ε f ε (xφ(x dx, φ C 0 (. (4.8 ( a1 ( ε u ε L 2 (;C 3 + ε u ε L 2 (;C 3 + u ε L 2 ( <. (4.9 8 This follows from noting the lower bound on the spectrum of the laplacian on the space of H 1 (D i functions that are i -quasi-periodic in direction y i. 17

19 Consequently, Proposition 4 informs us that a subsequence of u ε -quasi-periodic two-scale converges to some u L 2 (; H 1(Q, and moreover ε u 2 ε y u. Let us study the structure of this limit u further. We begin by introducing the densely defined unbounded linear operator a 1 : H 1(Q L2 (Q L 2 (Q; C 3 which is given by the action w a 1 w, for w H 1 (Q. We now argue that a generalised Weyl s decomposition holds, which was first introduced and proved for the case = 0 in [21]. Lemma 2. Let ( a 1 denote the adjoint of a 1. Then, the orthogonal decomposition holds. L 2 (Q; C 3 = ker ( ( a 1 Ran( a 1 Remark 7. Lemma 2 is a generalisation of the well-known fact that (periodic divergence-free vector fields are mutually orthogonal to gradients of (periodic potentials in L 2. In fact, this classical result can be deduced from the above lemma by (formally 9 setting a 1 = I on. Proof of Lemma 2. By the Banach closed ranged theorem, this result will follow if we demonstrate that the range of a 1 is closed, and this fact is implied by Proposition 8. Indeed, suppose u n Ran( a 1 converges strongly in L 2 (Q; C 3 to some u as n, i.e. there exists w n H 1(Q such that a 1 w n converges strongly in L 2 (Q; C 3 to u. By Proposition 8, the sequence wn, where wn denotes the orthogonal projection of w n onto the orthogonal complement V of V in H 1 (Q, is a Cauchy sequence in H 1(Q and therefore converges, up to some subsequence, to w H1 (Q. In particular, a1 w n = a 1 wn converges strongly in L 2 (Q to a 1 w and, consequently u = a 1 w. Hence, the range of a 1 is closed. Let us now describe u in detail. Lemma 3. The function u belongs to the Bochner space L 2 (; V. Proof. Recall that V = {v H 1 (Q v = 0 in Q1 } = ker( a 1, and so we aim to show that a 1 u = 0. On the one hand we deduce from (4.9 and (2.2 that lim ε a 1 ( ε 0 x ε u ε φ(xφ( x ε dx = 0. 2 Yet, on the other hand, Proposition 2 and the assertion ε u ε y u imply lim ε a 1 ( ε 0 x ε u ε φ(xφ( x ε dx = lim ε u ε φ(xa 1 ( x ε 0 ε Φ( x ε dx = a 1 (y y u(x, y φ(xφ(y dx. Therefore, as finite sums of φ(xφ(y are dense in L 2 ( Q; C 3 it follows that a 1 y u = 0 and since a1 1 L (Q we find that u L 2 (; V. The following result is of fundamental importance in characterising the (-quasi-periodic limit of the flux a 1 ( ε u ε in terms of the limit u of the function u ε. Put another way, this identity is crucial for determining the homogenised coefficients. 9 In fact, as expected the proof of this statement for I is much easier as I is positive where as a 1 is non-negative. Q 18

20 Lemma 4. There exists ξ L 2 ( Q; C 3 such that, up to a subsequence, a 1 ( ε u 2 ε ξ. Moreover, ξ belongs to the Bochner space L 2( ; ker ( ( a 1 and the pair (u, ξ satisfies the identity ξ(x, y φ(xψ(y dydx = a1 (yu(x, y x φ(x Ψ(y dydx, Q Q φ C (, Ψ ker ( ( a 1 (4.10. Proof. By Proposition 1 and (4.9 there exists ξ L 2 ( Q; C 3 such that, up to a subsequence that we discard, one has a1 ( ε u 2 ε ξ. (4.11 To prove ξ L 2( ; ker ( ( a 1, we take in (4.8 test functions of the form εφ(xϕ( x ε, φ C 0 (, ϕ C (Q, and use (4.9, (4.11 to pass to the limit in ε and deduce that a1 (yξ(x, y φ(x y ϕ(y dydx = 0. Q Therefore, for almost every x one has a1 (yξ(x, y y ϕ(y dy = 0, Q ϕ H 1 (Q, and, hence by Lemma 2 it follows that ξ(x, y L 2( ; ker ( ( a 1. Let us now prove assertion (4.10. Henceforth, we consider Ψ ker ( ( a 1 to be -quasi-periodically extended to R 3. We shall prove below the following integration by parts formula: a1 ( x ε u ε(x φ(xψ( x ε dx = a1 ( x ε u ε(x x φ(x Ψ( x ε dx, φ C (, Ψ ker ( ( a 1 ( Using Proposition 2, (4.11 and the convergence u ε u, we pass to the limit in the above formula to readily arrive at (4.10. To prove (4.12, it is sufficient to prove the following: for every w H 1 (R 3 one has a1 ( x ε w(x Ψ( x ε dx = 0. (4.13 R 3 Indeed, (4.12 follows from utilising (4.13 and the following facts: for φ C ( then u ε φ belongs to H 1 0 (, as u ε H 1 0 (, and can be trivially extended to H 1 (R 3, and that a1 ( x ε u ε(x φ(xψ( x ε dx = a1 ( x ε (u εφ(x Ψ( x ε dx a1 ( x ε u ε(x x φ(x Ψ( x ε dx. Let us now prove (4.13. For Q (z ε = 3 i=1 ε(z i, z i + 1, z Z 3, it follows that a1 ( x ε w(x Ψ( x ε dx a1 ( x ε w(x Ψ( x ε dx = R 3 z Z 3 = ε 3 z Z 3 Q (z ε Q a1 (y w(εy + εz exp(i zψ(y dy, 19

21 where the last equality comes from the change of variables x = ε(y + z and recalling that a 1 (y is periodic and Ψ is -quasi-periodic. By noting, for w H 1 (R 3, that w ε (y := z Z 3 w(εy + εzexp( i z, y Q, is an element of H 1 (Q, and that w ε (y := ε z Z 3 w(εy + εzexp( i z, y Q, the identity (4.13 follows. We are now ready to describe the properties of the macroscopic part of u and express the flux ξ in terms of u. Lemma 5. Let (u, ξ, u L 2 (; V and ξ L 2 (; ker ( ( a 1, be a pair which satisfies the identity (4.10. Then, u U, see (3.6. That is, for every for i I = {i I i = 0}, one has u = u i on C i, where i u i L 2 ( with u i ν i = 0 on, for ν the outer unit normal to. Furthermore, ξ(x, y = a 1 (y i I xi u i (x1 Ci (y[ y N (i (y + e i ], x, y Q 1. (4.14 Here, N (i solve (3.8. The following result immediately follows from the above lemma. Proposition 9. For every for i I, one has C i a1 (yξ(x, y dy = a hom i xi u i (xe i, for almost every x. Here, a hom are given by (3.7, i.e. a hom i = C i a 1 (y[ yi N (i (y + 1] dy > 0, for N (i H# 1 i (C i = {u H 1 (C i u is 1-periodic in the variable y i } is the unique solution to the cell problem a 1 (y [ N (i ] (y + e i φ(y dy = 0, φ H 1 #i (C i, C i (4.15 N (i = 0. C i Proof. Equation (4.14 implies a1 (yξ(x, y dy = xi u i (x a 1 (y[ y N (i (y + e i ] dy. C i C i For each j {1, 2, 3}\{i}, we set φ = y j in (4.15 to determine that a 1 (y [ y N (i ] (y + e i yj (y dy = a 1 (y [ y N (i ] (y + e i dy ej = 0. C i C i Hence, it follows that a1 (yξ(x, y dy = xi u i (x a 1 (y[ y N (i (y e i + 1]e i dy C i C i = a hom i xi u i (xe i, 20