Homogenization of a locally periodic oscillating boundary

TL;DR

This paper studies the homogenization process for Laplace equations in domains with oscillating boundaries, establishing weak and strong convergence results under natural regularity assumptions.

Contribution

It provides new homogenization results for mixed boundary value problems with locally periodic oscillating boundaries, including weak and strong convergence analyses.

Findings

Weak $L^2$ convergence of solutions and flows established

Strong $L^2$ convergence of extensions analyzed

Homogenization holds under natural regularity conditions

Abstract

This paper deals with the homogenization of a mixed boundary value problem for the Laplace operator in a domain with locally periodic oscillating boundary. The Neumann condition is prescribed on the oscillating part of the boundary, and the Dirichlet condition on a separate part. It is shown that the homogenization result holds in the sense of weak $L^2$ convergence of the solutions and their flows, under natural hypothesis on the regularity of the domain. The strong $L^2$ convergence of average preserving extensions of the solutions and their flows is also considered.

Figures

Tables

Table 1. Table 1. The numerically computed errors for the approximation in Theorem  6.1 , rounded to three significant digits. The number of degrees of freedom, rounded to two significant digits, is abbreviated to dof. The row ε = 0 𝜀 0 \varepsilon=0 is for the numerical solution u s 0 subscript superscript 𝑢 0 𝑠 u^{0}_{s} to the homogenized problem.

$\epsilon$	${\Vert u_s^{\epsilon }-u_s^0\Vert }_{L^2({\Omega }^{\epsilon },h)}$	${\Vert \nabla u_s^{\epsilon }-A^0\nabla u_s^0\Vert }_{L^2({\Omega }^{\epsilon },h)}$	dof
$0$	–	–	$78\cdot {10}^2$
$1/1$	$\mathrm{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}0.0488}$	$\mathrm{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}0.141}$	$39\cdot {10}^3$
$1/2$	$\mathrm{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}0.0502}$	$\mathrm{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}0.140}$	$46\cdot {10}^3$
$1/3$	$\mathrm{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}0.0336}$	$\mathrm{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}0.101}$	$33\cdot {10}^3$
$1/4$	$\mathrm{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}0.0226}$	$0.0918$	$27\cdot {10}^3$
$1/5$	$\mathrm{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}0.0149}$	$0.0830$	$59\cdot {10}^3$
$1/6$	$\mathrm{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}0.0104}$	$0.0736$	$69\cdot {10}^3$
$1/7$	$\mathrm{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}0.0120}$	$0.0781$	$77\cdot {10}^3$
$1/8$	$\mathrm{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}0.0115}$	$0.0784$	$11\cdot {10}^4$
$1/9$	$\mathrm{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}0.0101}$	$0.0738$	$13\cdot {10}^4$
$1/10$	$\mathrm{\hspace{0.17em}0.00941}$	$0.0689$	$15\cdot {10}^4$
$1/11$	$0.00940$	$0.0698$	$20\cdot {10}^4$
$1/12$	$0.00882$	$0.0704$	$23\cdot {10}^4$
$1/13$	$0.00801$	$0.0680$	$24\cdot {10}^4$
$1/14$	$0.00740$	$0.0649$	$35\cdot {10}^4$
$1/15$	$0.00716$	$0.0638$	$35\cdot {10}^4$
$1/16$	$0.00688$	$0.0642$	$37\cdot {10}^4$

Equations

$$\displaystyle\Omega^{\varepsilon}=\big{\{}x\in\mathbb{R}^{2}:0<x_{1}<1,\,\,0<x_{2}<\eta\big{(}x_{1},\frac{x_{1}}{\varepsilon}\big{)}\big{\}},$$

$${\Omega }^{\epsilon }$$

$$-\Delta u$$

$$u$$

$$\nabla u\cdot \nu$$

$$-\mathrm{div}(h{A}^0\nabla u)$$

$$u$$

$$h{A}^0\nabla u\cdot \nu$$

$${\eta }_{-}(x)$$

$$\Omega$$

$${\Omega }_{-}$$

$${\Omega }_{+}$$

$$A^0$$

$$Y(x)$$

$$h(x)$$

$${\Gamma }_a=\left\{(x_1,x_2)\text{leavevmode}|\text{leavevmode}{x}_1\in (0,1),\text{leavevmode}{x}_2={\eta }_{+}(x_1),\text{leavevmode}h(x_1,x_2)=0\right\}.$$

$$W(\Omega ,\Gamma )$$

$$\widetilde{u^{\epsilon }}\rightharpoonup h{u}^0\text{ weakly in }{L}^2(\Omega ),$$

$$\widetilde{\nabla u^{\epsilon }}\rightharpoonup h{A}^0\nabla u^0\text{ weakly in }{L}^2(\Omega ),$$

$$\Vert u^{\epsilon }-u^0{\Vert }_{L^2({\Omega }^{\epsilon },h)}\to 0,$$

$$\Vert \nabla u^{\epsilon }-A^0\nabla u^0{\Vert }_{L^2({\Omega }^{\epsilon },h)}\to 0,$$

$${\widetilde{u^{\epsilon }}}^m\to u^0\text{ strongly in }{L}^2(\Omega ,h),$$

$${\widetilde{\nabla u^{\epsilon }}}^m\to A^0\nabla u^0\text{ strongly in }{L}^2(\Omega ,h),$$

$$\widetilde{u^{\epsilon }}$$

$$\widetilde{\nabla u^{\epsilon }}$$

$$\displaystyle-{\mathop{\mathrm{div}}}\big{(}A\big{(}x,\frac{x_{1}}{\varepsilon}\big{)}\nabla u\big{)}$$

$$u$$

$$\displaystyle A\big{(}x,\frac{x_{1}}{\varepsilon}\big{)}\nabla u\cdot\nu$$

$$A(x,y)\xi \cdot \xi \ge C|\xi {|}^2,x\in \Omega ,y\in \mathbb{T},\xi \in {\mathbb{R}}^2.$$

$$-\mathrm{div}(A^0\nabla u)$$

$$u$$

$$A^0\nabla u\cdot \nu$$

$$A_{11}^0$$

$$A_{12}^0$$

$$A_{21}^0$$

$$A_{22}^0$$

$$-\Delta u^{\epsilon }$$

$$\displaystyle u^{\varepsilon}(x)\sim\big{(}u^{0}_{+}(x_{2},y)+\varepsilon u^{1}_{+}(x_{2},y)+\varepsilon^{2}u^{2}_{+}(x_{2},y)\big{)}\Big{|}_{\displaystyle y=\frac{x_{1}}{\varepsilon}},$$

$$\{(x,y):0<x_1<1,0<x_2<\eta (x_1,y)\}.$$

Full text

1. Introduction

This paper is concerned with the homogenization of a boundary value problem for the Laplace operator on a domain in $\mathbb{R}^{2}$ with locally periodic oscillating boundary. Specifically, the domain is given by

[TABLE]

where $\eta$ is a positive Lipschitz continuous function which is periodic in the second variable, and $\varepsilon$ is a small positive parameter (see Figure 1(a)). Certain further requirements are imposed on $\eta$ to ensure a particularly simple structure of the homogenized problem. On the oscillating part of the boundary the Neumann condition is prescribed, and on a separate part the Dirichlet condition, and the data are assumed to be $L^{2}$. We are interested in the asymptotic behavior of the solutions $u^{\varepsilon}$ to the boundary value problem, and their flows, as $\varepsilon$ tends to zero.

Domains with oscillating boundaries have attracted particular interest in the case where the domain is thin, that is to say when homogenization and dimension reduction may take place. This is natural in the mathematical physics program of the derivation of lower-dimensional theories from three-dimensional (see e.g. [15]). There is a rich literature on thin heterogeneous domains (see [47, 49, 7] and the references therein). Asymptotic analysis in thin domains with locally periodic oscillating boundary was conducted in for example [22, 45, 19, 6, 3, 7, 51, 8].

There are many works on homogenization in periodically oscillating domains with pillar type oscillations of fixed amplitude where the cross-section of each pillar is constant in the vertical direction. For the literature on pillar type oscillations, we refer to [29, 46, 32] and the references therein. Oscillating boundary domains with non-uniform cylindrical pillars, that is when the cross-sections of the pillars are varying in the vertical direction, have been considered in [30, 42, 1, 39]. In the mentioned works, the top boundary of the pillars have been assumed to be flat, that is the measure of the cross-section of each pillar at the maximum height is assumed to be positive, and also the base of each pillar assumed to be flat. There are few works on non-flat top boundaries, and [21, 2] stand out. In [21], the authors restrict the boundary graph functions to be smooth, periodic, and to have a unique maximum in each period. In [31], the authors consider an oscillating domain without explicit periodicity assumption and the base of each uniform pillar is allowed to be non-flat. Locally periodic flat pillar type domains were considered in [26], with respect to width and height. The works of Mel’nyk and his collaborators (c.f. [44, 41, 27, 28, 43, 33]) appear to have had a strong influence on later developments, after the initial works of Brizzi and Chalot [20, 21].

In this paper, the oscillating domain $\Omega^{\varepsilon}$ is a bounded region partially bounded by the graph of a locally periodic Lipschitz function $\eta(x_{1},x_{1}/\varepsilon)$, where $\eta:[0,1]\times\mathbb{T}\to\mathbb{R}$, $\mathbb{T}$ is the 1-torus, and $\varepsilon$ is a small parameter. These assumptions on $\eta$ ensure that the domain is connected and Lipschitz. In particular, the domain is not thin in the direction normal to oscillation, and not of pillar-type, and no assumptions are made on the flatness. Under these assumptions, the domain naturally becomes asymptotically disconnected (in the $x_{1}$ direction) between two curves, one that appears as a part of the limiting boundary and one that appears as an interior interface, as $\varepsilon$ tends to zero. The assumptions on $\eta$ guarantee that these curves are graphs of Lipschitz functions.

The analysis simplifies considerably with the Brizzi-Chalot condition of a single bump in each period. It appears to be worthwhile to remark that it is not generally true that as soon as there is more than one bump in each period, even if one restricts to the smallest possible periodicity cell, the asymptotic limit will not be decoupled (fast and slow variables). An example that shows that the connectedness of the sections is not necessary is given.

The expected influence of the domain oscillations on the asymptotic behavior of the solutions is that since the Laplace operator is local and the periodicity of the domain makes a region asymptotically become non-connected in the $x_{1}$ direction, $\partial/\partial x_{1}$ cannot be present in the homogenized equation in that region. This will show in the homogenized boundary value problem.

The result of this paper is the homogenization of the domain for the particular boundary value problem in which the heterogeneity is only in the domain. Three properties of homogenization are shown. Namely, (i) the weak $L^{2}$ convergence of the zero-extended solutions and their flows (Theorem 6.1), (ii) that the error of the zeroth approximation of the solutions and their flows converge strongly to zero in $L^{2}$ restricted to the oscillating domain (Theorem 7.1), and (iii) the strong $L^{2}$ convergence of average preserving extensions of the solutions and their flows (Theorem 8.1). In regard to (iii), in [21] reflection extensions were constructed and used, while the extension we use is the one used in [38, 26].

The analysis methods we use are standard techniques of asymptotic analysis and homogenization in particular. The method of homogenization is outlined in [35]. The method of periodic unfolding is described in [23, 24, 25], which is closely related to the notion of two-scale convergence [48, 4, 55], a generalization of weak convergence. Some works in which the unfolding method was used extensively in problems with oscillating boundary are [17, 18, 26, 1]. The homogenized problem is of degenerate elliptic type, as will be described below (c.f. [21]). The classical theory of Sobolev spaces for such is outlined in [36, 37]. For the method of asymptotic expansions we refer to [10, 12, 11, 50, 14, 52, 13, 53, 16, 9, 40].

The rest of this paper is organized as follows. The problem statement and the results are presented in Section 2. The homogenized problem is derived using formal asymptotic expansions of the solutions in Section 3. In Section 4, an example is provided that shows that the first term in the asymptotic expansion to $u^{\varepsilon}$ may be globally regular even in a degenerating case. In Section 5, the mean value property for the periodic unfolding is presented, adjusted to the present problem. The homogenization result of weak convergence of the solutions and their flows is established in Section 6. In Section 7, the convergence of energy is shown. In Section 8, the strong convergence of the extended solutions and their flows using average preserving extensions is shown. A numerical example, illustrating the rate of convergence is presented in Section 9. In Section 10, some information about a case of non-connected sections is provided.

2. Problem statement and results

In this section we state the boundary value problem and present the results.

With strictly positive Lipschitz $\eta:[0,1]\times\mathbb{T}\to\mathbb{R}$, the sequence of Lipschitz domains with periodically oscillating boundaries is for each $\varepsilon=1/k$, $k=1,2,\ldots$, defined by

[TABLE]

Here, $\mathbb{T}$ denotes the one-dimensional torus realized as $(0,1)$. See Figure 1(a) for an illustration of $\Omega^{\varepsilon}$. Additional restrictions on $\eta$ will be added below in order to ensure homogenization.

Let $u^{\varepsilon}\in H^{1}(\Omega^{\varepsilon},\Gamma)$ be the sequence of solutions to the following mixed boundary value problem:

[TABLE]

with $f$ in $L^{2}(\Omega)$, where $\Omega^{\varepsilon}\subset\Omega$. Here $\nu$ denotes the outward unit normal to the domain, and $H^{1}(\Omega^{\varepsilon},\Gamma)$ the functions in $H^{1}(\Omega^{\varepsilon})$ with zero trace on $\Gamma=(0,1)\times\{0\}$. Our goal is to describe the asymptotic behavior of the solutions $u^{\varepsilon}$ to (2) as $\varepsilon$ tends to zero.

The solutions $u^{\varepsilon}$ will be approximated in terms of the solution $u^{0}$ to the homogenized problem:

[TABLE]

where the coefficients $A^{0}$ and the domain $\Omega$ are defined as follows. In terms of the Lipschitz functions

[TABLE]

the domain

[TABLE]

is separated into the regions

[TABLE]

with interior interface $\Gamma_{-}=\partial\Omega_{-}\cap\partial\Omega_{+}$. The effective matrix is

[TABLE]

An illustration of $\Omega$, with the regions $\Omega_{+}$ and $\Omega_{-}$, and the interface $\Gamma_{-}$, indicated is shown in Figure 1(b), corresponding to the domain $\Omega^{\varepsilon}$ in Figure 1(a).

Let $h$ denote what we call the density of $\Omega^{\varepsilon}$ in $\Omega$:

[TABLE]

In $\Omega_{-}$ the density is $h=1$. A part of the upper boundary, denoted by $\Gamma_{a}$ is given by

[TABLE]

In order to ensure that homogenization takes place, the following hypotheses will be used:

1. (H1)

$Y(x)$ is connected, $x\in\Omega$. 2. (H2)

$Y(x)=k/N+Y(x)$, $x\in\Omega$, $k\in\mathbb{Z}$, for some natural number $N\geq 1$, and $Y(x)=y_{0}-Y(x)$ for some $y_{0}\in[0,1/N)$, and $Y_{0}(x)=\{y\in Y(x):y_{0}\leq y\leq y_{0}+1/(2N)\}$ is connected, $x\in\Omega$.

The hypothesis (H1) means that there is only one so-called pillar or bump in each period, and (H2) means restricting to a fundamental symmetry cell with respect to some translations and mirror symmetry. The hypothesis (H1) is stronger than (H2), and it gives a sufficient condition for homogenization. The hypothesis (H2) is included here to illustrate that it is not necessary there is only one ’bump’ in each period even if one uses the smallest periodicity cell. The graph of a function $\eta$ that satisfies (H2) but not (H1) is illustrated in Figure 3.

The assumption that $\eta$ is strictly positive ensures that the segment $\Gamma=(0,1)\times\{0\}$ is separated from the graph of $\eta(x_{1},x_{1}/\varepsilon)$, so $\Omega_{-}$ is a nonempty connected Lipschitz domain. The subdomains $\Omega_{+}$ and $\Omega_{-}$ have been chosen such that $\Omega_{+}$ covers the periodic region of $\Omega^{\varepsilon}$, and $\Omega_{+}$ is of positive measure if $\eta(x,y)$ is non-constant in $y$ for at least one $x$.

Denote by $L^{2}(\Omega,h)$ the Lebesgue space $\{v:\int_{\Omega}v^{2}h\,dx<\infty\}$, and $W(\Omega,\Gamma)$ the Sobolev space

[TABLE]

where $A^{0}$ is defined in (3) and $h$ is defined in (4)-(5). The homogenized problem (2) has a unique solution $u^{0}\in W(\Omega,\Gamma)$.

We will first establish the homogenization of (2), that is the convergence of the solutions $u^{\varepsilon}$ and their flows $\nabla u^{\varepsilon}$ to the solution $u^{0}$ to the homogenized problem (2) and its flow $hA^{0}\nabla u^{0}$, under hypothesis that (H1) holds.

Let $u^{\varepsilon}\in H^{1}(\Omega^{\varepsilon},\Gamma)$ be the solutions to (2), and let $u^{0}\in W(\Omega,\Gamma)$ be the solution to (2). In Theorem 6.1, it is shown under hypothesis (H1) that

[TABLE]

as $\varepsilon$ tends to zero. Here tilde denotes extension by zero.

At this point we know that there will be oscillations in $\widetilde{u^{\varepsilon}}$ in the upper part $\Omega_{+}$ due to the periodicity of $\Omega_{\varepsilon}$, while no oscillations in the lower part $\Omega_{-}$ due to the compact embedding of $H^{1}(\Omega_{-})$ into $L^{2}(\Omega_{-})$. Moreover, we cannot at this point exclude the possibility of oscillations in the solutions $\widetilde{u^{\varepsilon}}$ due to something else, as the above weak convergences may be expressed as weak unfolding or two-scale convergence. A next step is to check whether the strong unfolding convergence holds for the solutions $u^{\varepsilon}$ and their flows.

The weak convergence of the zero extensions in the upper part $\Omega_{+}$ cannot be strong, unless the limit is zero. To this end we first describe the error measured in the oscillating domain. Not only are the unfoldings of the solutions $u^{\varepsilon}$ and their flows strongly converging, there are no oscillations.

For $u^{\varepsilon}$ the solutions to (2) and $u^{0}$ the solution to (2), in Theorem 7.1 it is shown under hypothesis (H1) that

[TABLE]

as $\varepsilon$ tends to zero.

Turning back to the question of homogenization, there is an extension of the solutions $u^{\varepsilon}$ and their flows for which strong convergence holds in $L^{2}(\Omega)$.

For $u^{\varepsilon}$ the solutions to (2) and $u^{0}$ the solution to (2), in Theorem 8.1 it is shown that when the functions are extended in a way preserving their average, under hypothesis (H1) (c.f. [38, 26]),

[TABLE]

as $\varepsilon$ tends to zero, when $\sim{\!m}$ denotes the particular extension ((25), (26) in Section 8). It is remarked that the above mentioned convergences also hold under hypothesis (H2).

Without the hypotheses (H1), (H2), under a slightly milder restriction on the domain, in Theorem 10.1 it is shown that the solutions $u^{\varepsilon}$ to (2) converge to the solution $u^{0}$ to the limit problem (10), with $A^{0}$ given by (29), in the sense that

[TABLE]

as $\varepsilon$ tends to zero. One also has strong unfolding convergence. We do not analyze further the convergence in the case of non-connected sections.

In the case of the homogeneous Dirichlet condition on the oscillating part of the boundary, the limit is trivial in the oscillating part. This case was studied in [5].

Remark 2.1** (The effect of anisotropy and oscillating coefficients).**

To illustrate the effect of anisotropy and oscillating coefficients on the asymptotic behavior of the solutions $u^{\varepsilon}$ to (2), one can consider the following elliptic model problem:

[TABLE]

under the assumption that the not necessarily symmetric matrix $A\in C^{1}(\overline{\Omega},L^{\infty}(\mathbb{T}))$ satisfies the ellipticity condition

[TABLE]

Suppose further that hypothesis (H2) is satisfied, and that the coefficient matrix satisfies the corresponding symmetry conditions: $A(x,y)=A(x,k/N+y)$, $k\in\mathbb{Z}$, and $A(x,y)=A(x,y_{0}-y)$, for $x\in\Omega_{+}$, where $N$ and $y_{0}$ are the same as for $Y(x)$ in (H2). Then the weak limits of $\widetilde{u^{\varepsilon}}$ and $\widetilde{\nabla u^{\varepsilon}}$ are $hu^{0}$ and $A^{0}u^{0}$, respectively, in $L^{2}(\Omega)$, as $\varepsilon$ tends to zero, where $u^{0}$ solves the limit problem

[TABLE]

where the entries $A^{0}_{ij}$ of the effective matrix $A^{0}$ are given by, here including $h$,

[TABLE]

where $A_{ij}$ denote the entries of the local matrix $A$. In the fixed region $\Omega_{-}$, $A^{0}$ is the classical effective matrix for layered materials because there $Y(x)=\mathbb{T}$, while in the oscillating region $\Omega_{+}$ all but the last term in $A_{22}^{0}$ vanish due to the insulation in the $x_{1}$ direction. In particular, for $-\Delta$ the effective matrix reduces to the matrix given in (3). A similar statement can be made for the case of non-connected sections analogous to Theorem 10.1.

3. Asymptotic expansions for the solutions

To derive the homogenized equation (2) for the solutions $u^{\varepsilon}$ to the equation

[TABLE]

asymptotic expansions may be used, in the form of a formal power series in the small parameter $\varepsilon$. In the problem (7), in the periodic region $\Omega^{\varepsilon}_{+}$ of the domain $\Omega$, $\varepsilon$ represents both the periodicity in the $x_{1}$ direction, as well as the order of magnitude of the widths of the pillars with homogeneous Neumann condition on their sides. In the fixed region $\Omega_{-}$ of the domain, $\varepsilon$ represents the periodicity on the interface $\Gamma_{-}$.

The above reasoning leads us to consider the following Bakhvalov ansatz for inner expansion in the periodic region $\Omega^{\varepsilon}_{+}$:

[TABLE]

where $u^{i}_{+}$ are assumed to be periodic in $y$.

Let $\chi(x,y)$ denote the characteristic function of the set

[TABLE]

Then the characteristic function of $\Omega^{\varepsilon}$ is

[TABLE]

We write (7) in the homogeneous domain $\Omega$ as follows:

[TABLE]

With the ansatz (8),

[TABLE]

A substitution of (8) into (9) and collecting similar powers of $\varepsilon$ result in the following equations for the initial powers of $\varepsilon$:

[TABLE]

The equation for power $\varepsilon^{-2}$ suggests that $u^{0}_{+}$ is independent of $y$, under the assumption that $Y(x)$ is connected, or some discrete symmetry in the problem such as (H2). Viewing the above equations as definitions of $u^{i}_{+}$ in $\mathbb{T}$ with $x$ as a parameter, the compatibility condition for $u^{2}_{+}$ may be read off from the $\varepsilon^{0}$ equation:

[TABLE]

Because

[TABLE]

where

[TABLE]

the compatibility condition is

[TABLE]

The equation (10) is the homogenized equation for $u^{\varepsilon}$ in the periodic part $\Omega_{+}$ of $\Omega$.

In the fixed region $\Omega_{-}$ of the domain, the equation (7) reads

[TABLE]

Consider the following ansatz for inner expansion in the fixed region $\Omega_{-}$:

[TABLE]

where $u^{i}_{-}$ are assumed to be periodic in $y$. With the ansatz (12),

[TABLE]

A substitution of (12) into (11) and collecting similar powers of $\varepsilon$ result in the following equations for the initial powers of $\varepsilon$:

[TABLE]

The equation for power $\varepsilon^{-2}$ suggests that $u^{0}_{-}$ is independent of $y$. By the periodicity of $u^{1}_{-}$, the compatibility condition for $u^{2}_{-}$ reads

[TABLE]

In order $u^{0}=\chi_{\Omega_{-}}u^{0}_{-}+\chi_{\Omega_{+}}u^{0}_{+}$ to be globally defined, an interface condition is needed on $\Gamma_{-}=\partial\Omega_{-}\cap\partial\Omega_{+}$. If one requires continuity of $u^{0}$ and its flow on $\Gamma_{-}$, one must have $u^{0}_{+}=u^{0}_{-}$ on $\Gamma_{-}$, as well as

[TABLE]

where $\nu$ is one of the unit normals on $\Gamma_{-}$, and the upper left entry of the matrix in the first term is arbitrary. This condition becomes explicitly,

[TABLE]

which may be expressed as

[TABLE]

where $[\cdot]$ denotes the jump on $\Gamma_{-}$, and $A^{0}$ is given by (3).

4. Example of behavior of the solutions in a degenerating case

The solution $u^{0}$ to the homogenized problem (2) may belong to $H^{1}(\Omega)$ even if $h(x)$ tends to zero as $x$ approaches $M=\{x:h(x)=0\}\subset\partial\Omega$, as the following example illustrates. Consider the case $f=1$ with $\eta(x,y)=\hat{\eta}(y)$ the piecewise linear function given by $\hat{\eta}(0)=1$ and $\hat{\eta}(1/2)=2$, that is $\hat{\eta}(y)=2-2|y-1/2|$. The homogeneous domains are $\Omega=(0,1)\times(0,2)$, $\Omega_{+}=(0,1)\times(1,2)$, $\Omega_{-}=(0,1)\times(0,1)$, and the interface is $\Gamma_{-}=(0,1)\times\{1\}$. In Figure 2, illustrations of the domains $\Omega^{\varepsilon}$ and $\Omega$ are shown. In this case,

[TABLE]

and

[TABLE]

The solution to the homogenized problem (2) is

[TABLE]

In particular, $u^{0}\in H^{1}(\Omega)$ and it has continuous gradient over the interface $\Gamma_{-}$.

5. Periodic unfolding

The only apparent possible cause of oscillations in the solutions to (2) and their flows is the periodicity in the domain, in the $x_{1}$ direction. For the study of these oscillations we will use periodic unfolding.

The periodic unfolding at rate $\varepsilon$ of a function $v:\mathbb{R}^{2}\to\mathbb{R}$ along $x_{1}$ is

[TABLE]

where $[\cdot]$ denotes the integer part, and $v$ is extended by zero when necessary. Using this change of variables for the periodic coefficient $\chi_{\Omega^{\varepsilon}_{+}}$ in (2), the characteristic function of the region of $\Omega^{\varepsilon}$ where coefficients are periodic,

[TABLE]

gives $\chi_{\Omega^{\varepsilon}_{u}}$, the characteristic function of the domain

[TABLE]

There holds

[TABLE]

where

[TABLE]

The property (13) is the strong unfolding convergence of the sequence and it will be used in the passage from periodic domain to a fixed domain in integrals. It expresses that $\chi_{\Omega_{+}^{\varepsilon}}$ converges weakly in $L^{p}(\mathbb{R}^{2})$ while not strongly, and that the oscillation spectrum of the sequence belongs to the integers if not empty. To obtain (13) one uses the almost everywhere pointwise convergence of $\chi_{\Omega^{\varepsilon}_{u}}$ to $\chi_{\Omega_{u}}$ and the Lebesgue dominated convergence theorem, or views it as a consequence of Lemma 5.1 below.

The cost of replacing in integrals the $\varepsilon$ depending unfolded domain $\Omega^{\varepsilon}_{u}$ with the fixed domain $\Omega_{u}$ is described by the following lemma.

Lemma 5.1**.**

Let $\Omega$ contain $\Omega^{\varepsilon}_{+}$. Suppose that $\|v^{\varepsilon}\|_{L^{p}(\Omega)}\leq C$ and $p>1$. Then

[TABLE]

as $\varepsilon$ tends to zero.

Proof.

Because $T^{\varepsilon}\chi_{\Omega_{+}^{\varepsilon}}=\chi_{\Omega_{u}^{\varepsilon}}$, the discrepancy can be computed as follows:

[TABLE]

The differences $\Omega_{u}^{\varepsilon}\setminus\Omega_{u}$ and $\Omega_{u}\setminus\Omega_{u}^{\varepsilon}$ are contained in some strip of measure $O(\varepsilon)$:

[TABLE]

where $C$ may be chosen independent of $\varepsilon$ by the Lipschitz continuity of $\eta$ and $\eta_{-}$,

By the Hölder inequality,

[TABLE]

which gives the desired estimate. ∎

6. Homogenization

In this section we establish the homogenization of problem (2) to (2) in the sense of weak convergence of the solutions and their flows. The method we use is the unfolding Lemma 5.1 to pass to the fixed domain $\Omega_{u}$, and the weak compactness in $L^{2}(\Omega_{u})$ to characterize the asymptotic behavior of $u^{\varepsilon}$.

Throughout this section (H1) is assumed to hold.

Theorem 6.1**.**

Suppose that (H1) holds. Let $u^{\varepsilon}\in H^{1}(\Omega^{\varepsilon},\Gamma)$ be the solutions to (2), and let $u^{0}\in W(\Omega,\Gamma)$ be the solution to (2). Then

[TABLE]

as $\varepsilon$ tends to zero, where $\sim$ denotes extension by zero.

The convergence of the flows in Theorem 6.1(ii) means that $u^{\varepsilon}$ converges weakly to $u^{0}$ in $H^{1}(\Omega_{-},\Gamma)$, and strongly in $L^{2}(\Omega_{-})$ by the Relich theorem.

Lemma 6.1**.**

For any $\varepsilon=1/k$, $k=1,2,\ldots$, there exists a unique solution $u^{\varepsilon}\in H^{1}(\Omega^{\varepsilon},\Gamma)$ to (2). For the solutions $u^{\varepsilon}$, the following a priori estimate holds:

[TABLE]

where $C$ is independent of $\varepsilon$.

Proof.

The variational form of (2) is: Find $u^{\varepsilon}\in H^{1}(\Omega^{\varepsilon},\Gamma)$ such that

[TABLE]

for all $\psi\in H^{1}(\Omega^{\varepsilon},\Gamma)$. Using the Poincaré inequality

[TABLE]

one verifies that left hand side in (14) is an inner product on $H^{1}(\Omega^{\varepsilon},\Gamma)$. The right hand side in (14) is a bounded linear functional on $H^{1}(\Omega^{\varepsilon},\Gamma)$. The Riesz theorem guarantees the existence of a unique solution $u^{\varepsilon}\in H^{1}(\Omega^{\varepsilon},\Gamma)$. Using $u^{\varepsilon}$ as a test function in (14) gives

[TABLE]

uniformly in $\varepsilon$ by the uniform Poincaré constant, from which the desired a priori estimate is obtained. ∎

Lemma 6.2**.**

For the sequence of solutions $u^{\varepsilon}$ to (2), the following a priori estimates hold for the unfolded sequences:

[TABLE]

where $C$ is independent of $\varepsilon$.

Proof.

By the definition of unfolding,

[TABLE]

because $T^{\varepsilon}u^{\varepsilon}$ is zero outside $\Omega_{u}^{\varepsilon}$. By using the estimate for the solutions $u^{\varepsilon}$ in Lemma 6.1, the estimate $\|T^{\varepsilon}u^{\varepsilon}\|_{L^{2}(\Omega_{u})}\leq C$ is obtained. The same computation with $|\nabla u^{\varepsilon}|$ in place of $u^{\varepsilon}$ gives the estimate $\|T^{\varepsilon}\nabla u^{\varepsilon}\|_{L^{2}(\Omega_{u})}\leq C$. ∎

Regarding the function space, $W(\Omega,\Gamma)$ is associated to the homogenized problem (2). In the cover $\Omega_{+}$ of the periodic region of $\Omega^{\varepsilon}$, the $H^{1}(\Omega)$ ellipticity of the operator to the homogenized problem (2) may be violated. For

[TABLE]

and $h(x)$ tends to zero as $x$ approaches $M=\{x:h(x)=0\}$, which might be nonempty. The properties

[TABLE]

ensure that $W(\Omega,\Gamma)$ is a Hilbert space when equipped with the inner product $(u,v)=\int_{\Omega}A^{0}\nabla u\cdot\nabla v\,dx$, and that $C^{\infty}_{0}(\Omega)$ is embedded into $W(\Omega,\Gamma)$. We denote by $C^{\infty}(\Omega,\Gamma)$, the $C^{\infty}(\Omega)$ functions vanishing in a neighborhood of $\Gamma$.

The existence of a solution $u^{0}\in W(\Omega,\Gamma)$ to $\eqref{eq:limitproblem}$ is obtained by weak compactness in the proof of Theorem 6.1 below, with uniqueness by linearity. One might also obtain it in the direct way as follows.

Lemma 6.3**.**

There exists a unique solution $u^{0}\in W(\Omega,\Gamma)$ to $\eqref{eq:limitproblem}$.

Proof.

The variational form of (2) is: Find $u^{0}\in W(\Omega,\Gamma)$ such that

[TABLE]

for all $\psi\in W(\Omega,\Gamma)$.

The Hilbert space structure on $W(\Omega,\Gamma)$ has been chosen such that

[TABLE]

By the Hölder and Poincaré inequalities, the left hand side of (15) defines an inner product on $W(\Omega,\Gamma)$. The right hand side of (15) is a bounded linear functional on $W(\Omega,\Gamma)$. By the Riesz theorem, there exists a unique $u^{0}\in W(\Omega,\Gamma)$ satisfying (15). ∎

Proof of Theorem 6.1.

From Lemma 6.1 and Lemma 6.2 we have the following a priori estimates for the solutions $u^{\varepsilon}$ to problem (2) and their unfoldings:

[TABLE]

By weak compactness, there exist $u^{0}_{-}\in H^{1}(\Omega_{-},\Gamma)$, $u^{0}_{+}\in L^{2}(\Omega_{u})$, $p\in L^{2}(\Omega_{u})$, and a subsequence of $\varepsilon$ which we still denoted by $\varepsilon$, such that

[TABLE]

where the equality of second component of the weak limit of $T^{\varepsilon}\nabla u^{\varepsilon}$ and $\frac{\partial u^{0}_{+}}{\partial x_{2}}$, and that $u^{0}_{+}$ does not depend on $y$, follow from the boundedness of the sequences and

[TABLE]

By (H1) and (19), $u^{0}_{+}$ does not depend on $y$.

**Claim 1: The average of $p$ in $y$ is zero: $\int_{Y(x)}p\,dy=0$, a.e. $x\in\Omega_{+}$.

**

Information about $p$ may be obtained by using oscillating test functions in the equation for $u^{\varepsilon}$ (c.f. [40, 54]). The prototype $\varphi^{\varepsilon}=\varepsilon\phi(x)\{x_{1}/\varepsilon\}$, with $\phi\in C^{\infty}_{0}(\Omega_{+})$ and $\{\cdot\}$ denoting fractional part, would serve the purpose in view of (22) because $T^{\varepsilon}\varphi^{\varepsilon}\to 0$ and $T^{\varepsilon}\nabla\varphi^{\varepsilon}\to(\phi,0)$ strongly in $L^{2}(\Omega_{u})$. As $\varphi^{\varepsilon}$ are not necessarily continuous on $\Omega^{\varepsilon}_{+}$, appropriate shifts are introduced as follows. Note that

[TABLE]

Replace $x^{\varepsilon}_{k}=k\varepsilon$ with a grid where the graphs of $\eta(x_{1},x_{1}/\varepsilon)$ and $\eta_{-}(x_{1})$ are close:

[TABLE]

Let $\phi\in C^{\infty}_{0}(\Omega_{+})$. Then with $x_{k}^{\varepsilon}$ as in (20),

[TABLE]

belongs to $C^{\infty}(\Omega_{+}^{\varepsilon})$ for all small enough $\varepsilon$. Indeed,

[TABLE]

by the Lipschitz continuity of $\eta_{-}$, and $\phi$ is compactly supported.

With $\varphi^{\varepsilon}$ given by (20), (21) as test functions in the equation (22) for $u^{\varepsilon}$, and using that

[TABLE]

strongly in $L^{2}(\Omega_{u})$ as $\varepsilon$ tends to zero because $|x_{1}-x_{k}^{\varepsilon}|\leq\varepsilon$, one obtains in the limit

[TABLE]

which gives the first claim.

**Claim 2: $\chi_{\Omega_{-}}u^{0}_{-}+\chi_{\Omega_{+}}u^{0}_{+}\in W(\Omega,\Gamma)$ and is the solution to the homogenized problem (2).

**

To verify that $\chi_{\Omega_{-}}u^{0}_{-}+\chi_{\Omega_{+}}u^{0}_{+}\in W(\Omega,\Gamma)$ it suffices to check that $u^{0}_{-}=u^{0}_{+}$ in $L^{2}(\Gamma_{-})$. Note that $u^{0}_{+}\in L^{2}(\Omega_{u})$ is continuously traced into $L^{2}(\Gamma_{-})$ because $h>0$ in a neighborhood of $\Gamma_{-}$ in $\Omega$ (see Remark 6.1).

Because $u^{\varepsilon}\rightharpoonup u^{0}_{-}$ weakly in $H^{1}(\Omega_{-})$, the weak continuity of the trace gives $u^{\varepsilon}\rightharpoonup u^{0}_{-}$ weakly in $L^{2}(\Gamma_{-})$, and so $T^{\varepsilon}u^{\varepsilon}(x_{1},\eta_{-}(x_{1}))\rightharpoonup u^{0}_{-}(x_{1},\eta_{-}(x_{1}))$ weakly in $L^{2}((0,1)\times\mathbb{T})$, as $\varepsilon$ tends to zero. Denote $\Gamma_{\varepsilon}=\Omega_{\varepsilon}\cap\Gamma_{-}$. Let $\phi\in C^{\infty}_{0}(\Omega)$. On the one hand,

[TABLE]

as $\varepsilon$ tends to zero, because $T^{\varepsilon}\chi_{\Gamma_{\varepsilon}}(x_{1},\eta_{-}(x_{1}))$ converges to $\chi_{Y(x_{1},\eta_{-}(x_{1}))}(y)$ strongly in $L^{2}((0,1)\times\mathbb{T})$. On the other hand,

[TABLE]

as $\varepsilon$ tends to zero. Thus

[TABLE]

It follows that $u^{0}_{-}=u^{0}_{+}$ in $L^{2}(\Gamma_{-})$, for by the Lipschitz continuity of $\eta_{-}$,

[TABLE]

A split of the variational form (14) of problem (2) into the periodic $\Omega^{\varepsilon}_{+}$ and the fixed $\Omega_{-}$ reads

[TABLE]

After unfolding $\Omega^{\varepsilon}_{+}$ and using Lemma 5.1 one arrives with $\psi\in C^{\infty}(\overline{\Omega},\Gamma)$ at

[TABLE]

By passing to the limit in (22) with $\psi\in C^{\infty}(\overline{\Omega},\Gamma)$, using the weak convergence of $u^{\varepsilon}$, $T^{\varepsilon}u^{\varepsilon}$, $T^{\varepsilon}\nabla u^{\varepsilon}$ (16)–(18), and that $p$ is of average zero in $y$, one obtains that $\chi_{\Omega_{-}}u^{0}_{-}+\chi_{\Omega_{+}}u^{0}_{+}$ satisfies

[TABLE]

as $\varepsilon$ tends to zero.

Let

[TABLE]

and

[TABLE]

By the density of $C^{\infty}(\overline{\Omega},\Gamma)$ in $W$ under hypothesis (H1), (23) holds for any test function in $W(\Omega_{U},\Gamma\times\mathbb{T})$:

[TABLE]

for any $\psi\in W(\Omega_{U},\Gamma\times\mathbb{T})$. The equation (24) is well-posed in the Hilbert space $W(\Omega_{U},\Gamma\times\mathbb{T})$. Under (H1), the problem (24) is equivalent to (2), and $v$ belongs to the weighted space $W(\Omega,\Gamma)$ if and only if it belongs to $W(\Omega_{U},\Gamma\times\mathbb{T})$. One concludes that $u^{0}=\chi_{\Omega_{-}}u^{0}_{-}+\chi_{\Omega_{+}}u^{0}_{+}$.

**Claim 3: $\widetilde{u^{\varepsilon}}\rightharpoonup hu^{0}$ weakly in $L^{2}(\Omega)$.

**

The uniqueness of the solution $u^{0}$ to the homogenized problem (2) ensures that the full $\varepsilon$ sequences (16)–(18) converge.

The weak limit of a sequence is obtained from the weak unfolding limit (weak two-scale limit) by taking the average over the cell of periodicity. Because $T^{\varepsilon}u^{\varepsilon}$ converges weakly to $u^{0}$, and $u^{0}$ does not depend on $y$,

[TABLE]

weakly in $L^{2}(\Omega_{+})$, as $\varepsilon$ tends to zero. It follows that $\widetilde{u^{\varepsilon}}\rightharpoonup hu^{0}$ weakly in $L^{2}(\Omega)$, as $\varepsilon$ tends to zero.

**Claim 4: $\widetilde{\nabla u^{\varepsilon}}\rightharpoonup hA^{0}\nabla u^{0}$ weakly in $L^{2}(\Omega)$.

**

By the same property of weakly converging unfolding, as was used in the previous paragraph, and that $p$ is of average zero in $y$,

[TABLE]

weakly in $L^{2}(\Omega_{+})$, as $\varepsilon$ tends to zero. It follows that $\widetilde{\nabla u^{\varepsilon}}\rightharpoonup hA^{0}\nabla u^{0}$ weakly in $L^{2}(\Omega)$, as $\varepsilon$ tends to zero. ∎

Remark 6.1** (About the $L^{2}$ trace).**

Let $\gamma(v)(x)=v(x)$ for $x\in\Gamma_{-}$, $v\in C^{\infty}(\overline{\Omega_{+}})$. Then $\gamma:\{v\in L^{2}(\Omega_{+}):\frac{\partial v}{\partial x_{2}}\in L^{2}(\Omega_{+})\}\to L^{2}(\Gamma_{-})$ is linear and bounded, where the natural Hilbert space structures are employed. Indeed, for $v\in C^{\infty}(\overline{\Omega_{+}})$,

[TABLE]

By the triangle inequality for the integral, the Bunyakovsky-Cauchy-Schwarz inequality, and the inequality of arithmetic and geometric means,

[TABLE]

An integration in $x_{2}$ over the interval $(\eta_{-}(x_{1}),\eta_{+}(x_{1}))$ gives

[TABLE]

An integration in $x_{1}$ over the interval $(0,1)$ gives

[TABLE]

The set $C^{\infty}(\overline{\Omega_{+}})$ is dense in the domain of $\gamma$.

7. Justification

In this section it is shown that the error in approximating the solutions $u^{\varepsilon}$ to (2) and their flows with the solution $u^{0}$ to the homogenized problem (2) and its flow tends to zero, measured in $L^{2}(\Omega^{\varepsilon},h)$. The method we use is the convergence of energy.

Throughout this section (H1) is assumed to hold.

Theorem 7.1**.**

Suppose that (H1) holds. Let $u^{\varepsilon}\in H^{1}(\Omega^{\varepsilon},\Gamma)$ be the solutions to (2), and let $u^{0}\in W(\Omega,\Gamma)$ be the solution to (2). Then

[TABLE]

as $\varepsilon$ tends to zero.

Proof.

By the weak convergence of $u^{\varepsilon}$, $T^{\varepsilon}u^{\varepsilon}$, $T^{\varepsilon}\nabla u^{\varepsilon}$ (16)–(18), the property that sum of lim inf is less than or equal to lim inf of sum, using that $u^{\varepsilon}$, $u^{0}$ solve (22), (15),

[TABLE]

This shows that the energies converge. It follows that each weak convergence in (16)–(18) is strong, and $p=0$. ∎

As $T^{\varepsilon}\frac{\partial u^{\varepsilon}}{\partial x_{1}}$ tends to zero strongly in $L^{2}(\Omega_{u})$, the proof of Theorem 7.1 shows that

[TABLE]

as $\varepsilon$ tends to zero, which justifies the first term in the asymptotic expansion under hypothesis (H1).

8. Average preserving extension

In this section we prove the homogenization Theorem 6.1 for an average preserving extension (see [38, 26]). The method we use is the strong unfolding convergence obtained in the proof of Theorem 7.1, which is stated explicitly in the form of Lemma 8.1 below.

Throughout this section (H1) is assumed to hold. By the end of the section, a remark about (H2) is included.

Let the local average of a function $v$ in the $x_{1}$ direction over $\Omega^{\varepsilon}_{+}$ be denoted by

[TABLE]

where $m_{\varepsilon}(v)(x)$ is set to zero at points where $h(x)=|Y(x)|=0$. The average preserving extension is defined by

[TABLE]

For $v\in H^{1}(\Omega^{\varepsilon})$, $m_{\varepsilon}(v)$ and $m_{\varepsilon}(\nabla v)$ are well-defined, using the Sobolev space property that $v$ has a representative that is absolutely continuous on almost all line segments parallel to the coordinate axes and with square integrable partial derivatives. In particular, $v$ belongs to $H^{1}(\{(t,x_{2}):t\in\mathbb{R}\}\cap\Omega^{\varepsilon}_{+})$ for almost all $x_{2}$ of relevance.

Theorem 8.1**.**

Suppose that (H1) holds. Let $u^{\varepsilon}\in H^{1}(\Omega^{\varepsilon},\Gamma)$ be the solutions to (2), and let $u^{0}\in W(\Omega,\Gamma)$ be the solution to (2). Then

[TABLE]

as $\varepsilon$ tends to zero, where $\sim{\!m}$ denotes the extension (25), (26).

In the proof of Theorem 7.1, the following strong convergence of the unfolded sequences were obtained, which in terms other than unfolding means that the oscillation spectrum of the sequences belong to the integers if not empty, and that the sequences converge strongly two-scale.

Lemma 8.1**.**

Suppose that (H1) holds. Let $u^{\varepsilon}\in H^{1}(\Omega^{\varepsilon},\Gamma)$ be the solutions to (2), and let $u^{0}\in W(\Omega,\Gamma)$ be the solution to (2). Then

[TABLE]

as $\varepsilon$ tends to zero.

Proof of Theorem 8.1.

Because the extension does not alter the functions $u^{\varepsilon}$ in $\Omega_{-}$, the $\Omega_{-}$ parts of (i) and (ii) are included in Theorem 7.1. In $\Omega_{+}$, by definition,

[TABLE]

The first term on the right hand side tends to zero as $\varepsilon$ tends to zero according to Theorem 7.1. The second term may be estimated as follows. Because $u^{0}$ is independent of $y$, and $T^{\varepsilon}u^{\varepsilon}$ converges to $u^{0}$ strongly in $L^{2}(\Omega_{u})$, by Theorem 6.1 and Lemma 8.1(i), the Hölder inequality gives

[TABLE]

as $\varepsilon$ tends to zero, which gives part (i).

Part (ii) is obtained by repeating the above two steps for the components of $\nabla u^{\varepsilon}$ in place of $u^{\varepsilon}$, using the second parts of Theorem 7.1 and Lemma 8.1. ∎

Remark 8.1** (Hypothesis (H2)).**

The theorems 6.1, 7.1, 8.1, and Lemma 8.1, hold under the slightly weaker hypothesis (H2). The weak limit $u_{+}^{0}$ of $T^{\varepsilon}u^{\varepsilon}$ is the first component of the unique solution $(u_{+}^{0},u_{-}^{0})$ in the Sobolev space

[TABLE]

to the following problem: Find $(u_{+},u_{-})\in W$ such that

[TABLE]

for any $(\varphi,\phi)\in W$. This statement relies on the density of $C^{\infty}(\overline{\Omega},\Gamma)$ in $W$ (c.f. Section 10 below).

By (H2), $(u_{+}^{0}(x,k/N+y),u_{-}^{0})$, $k\in\mathbb{Z}$, and $(u_{+}^{0}(x,y_{0}-y),u_{-}^{0})$ solve the problem (27). By uniqueness of solution, $u_{+}^{0}(x,y)=u_{+}^{0}(x,k/N+y)$, $k\in\mathbb{Z}$, and $u_{+}^{0}(x,y)=u_{+}^{0}(x,y_{0}-y)$. Because all $Y_{0}(x)$ are connected, $\frac{\partial u^{0}_{+}}{\partial y}=0$ implies that $u^{0}_{+}$ is constant in $y$.

9. A numerical example

To illustrate the rate of convergence in Theorem 7.1 under hypothesis (H2), we consider the following example. Let $u^{\varepsilon}$ be the solutions to (2) with $f=1$ and

[TABLE]

where a bump is introduced at $y_{0}$ by

[TABLE]

with $\delta=1/10$ and $y_{0}=1/2$, $y\in[0,1)$. Then $\eta$ satisfies (H2) but not (H1). Moreover, $h=1$ on the graph $\Gamma_{-}$ of $\eta_{-}$, and $h=0$ on the graph of $\eta_{+}$. The domains $\Omega^{\varepsilon}$ and $\Omega$ are illustrated in Figure 3.

The solutions $u^{\varepsilon}$, $u^{0}$ to the problems (2), (2) are approximated by means of the finite element method using piecewise linear Lagrange elements. The numerical approximations are denoted by $u^{\varepsilon}_{s}$, $u^{0}_{s}$. The numerically computed rates of convergence for the approximation in Theorem 6.1, are illustrated in Figure 5, obtained using the numerical tool FreeFEM [34]. The data points were obtained using the values of $\varepsilon$ and the number of degrees of freedom given in Table 1.

One observes that the rate of convergence for $u^{\varepsilon}$ appears to be close to $\varepsilon^{3/4}$, and for $\nabla u^{\varepsilon}$ close to $\varepsilon^{1/4}$, for the approximation measured in $L^{2}(\Omega^{\varepsilon},h)$ for the selected values of $\varepsilon$.

10. A case of non-connected sections

In this section we consider a case where the sections $Y(x)$ are allowed to be disconnected, that is the hypotheses (H1) and (H2) are mildly relaxed. First we describe the weak unfolding limit of the zero extended solutions to problem (2), and then we describe the weak limit in terms of local domain densities using the ideas of Mel’nyk, splitting the oscillating part of the domain into branches (c.f. [27, 28, 43]).

Here, the full unfolded limit domain will be used, denoted by

[TABLE]

and the effective matrix is given by

[TABLE]

The limit problem reads

[TABLE]

where $\nu_{x}$ denotes the projection of the outward unit normal to $\Omega_{U}$ onto $\mathbb{R}^{2}$. The limit problem (10) is well-posed in the Sobolev space

[TABLE]

equipped with the natural Hilbert space structure.

Theorem 10.1**.**

Let $u^{\varepsilon}\in H^{1}(\Omega^{\varepsilon},\Gamma)$ be the solutions to problem (2). Let $u^{0}\in W(\Omega_{U},\Gamma\times\mathbb{T})$ be the solution to (10). Then

[TABLE]

as $\varepsilon$ tends to zero.

Proof.

The following a priori estimate holds

[TABLE]

There exist

[TABLE]

such that, along a subsequence still denoted by $\varepsilon$,

[TABLE]

At this point $u_{+}^{0}$ and $p$ depend in general on the fast variable $y$.

The unfolded variational form of problem (7) is by Lemma 5.1,

[TABLE]

as $\varepsilon$ tends to zero, for any sufficiently smooth $\psi\in H^{1}(\Omega^{\varepsilon},\,\Gamma)$.

**Step 1: $p=0$ a.e. in $\Omega_{u}$.

**

We will show that $p=0$ almost everywhere in $\Omega_{u}$. Let $\phi\in C^{\infty}_{0}(\Omega_{u})$ and consider the sequence of test functions \varphi^{\varepsilon}(x)=\varepsilon\phi\big{(}x,\frac{x_{1}}{\varepsilon}\big{)}. Then

[TABLE]

strongly in $L^{2}(\Omega_{u})$, as $\varepsilon$ tends to zero. By passing to the limit in (31) with test functions $\varphi^{\varepsilon}$ one obtains

[TABLE]

as $\varepsilon$ tends to zero. That is, $\nabla_{y}p=0$ in $\Omega_{u}$.

Let $\phi\in C^{\infty}_{0}(\Omega_{u})$ be such that $\nabla_{y}\phi=0$. Then consider the sequence of test functions $\varphi^{\varepsilon}$ defined by

[TABLE]

where $x^{\varepsilon}_{k}$ is a tagging such as in (20). Then

[TABLE]

strongly in $L^{2}(\Omega_{u})$, as $\varepsilon$ tends to zero. By passing to the limit in (31) with test functions $\varphi^{\varepsilon}$ one obtains

[TABLE]

as $\varepsilon$ tends to zero. By the density of $C^{\infty}_{0}(\Omega_{u},\nabla_{y})$ in the Sobolev space $\{v\in L^{2}(\Omega_{u}):\nabla_{y}v=0\text{ in }\Omega_{u}\}$, one concludes that (32) holds for $\phi=p$, that is

[TABLE]

**Step 2: Transmission condition.

**

To determine the transmission condition for $u^{0}_{-}$ and $u^{0}_{+}$ on the internal interface

[TABLE]

we will trace $u^{\varepsilon}$ from either side.

Let $\phi\in C^{\infty}_{0}(\Omega_{U})$, and set $\phi^{\varepsilon}(x,\frac{x_{1}}{\varepsilon})$.

On the one hand, from below,

[TABLE]

as $\varepsilon$ tends to zero, because $T^{\varepsilon}\chi_{\Gamma_{\varepsilon}}(x_{1},\eta_{-}(x_{1}))$ converges to $\chi_{Y(x_{1},\eta_{-}(x_{1}))}(y)$ strongly in $L^{2}((0,1)\times\mathbb{T})$.

On the other hand, from above,

[TABLE]

as $\varepsilon$ tends to zero.

Thus

[TABLE]

It follows that

[TABLE]

for by the Lipschitz continuity of $\eta_{-}$,

[TABLE]

The transmission condition (33) means here that the trace of $u^{0}_{+}$ on each connected component of the sections $Y(x)$ for $x\in\Gamma_{-}$ is equal to the trace of $u^{0}_{-}$. That is, zero jump condition into each branch from below.

**Step 3: Limit problem.

**

Let $\varphi\in C^{\infty}(\overline{\Omega_{U}})$ be such that $\nabla_{y}\varphi=0$, and $\varphi$ vanishes on $\Gamma\times\mathbb{T}$. Set $\psi^{\varepsilon}(x)=\varphi(x,\frac{x_{1}}{\varepsilon})$, which belongs to $H^{1}(\Omega^{\varepsilon},\Gamma)$. Using $\psi^{\varepsilon}$ as test functions in (31), one obtains

[TABLE]

in the limit as $\varepsilon$ tends to zero, where $u=\chi_{\Omega_{u}}u_{+}^{0}+\chi_{\Omega_{-}\times\mathbb{T}}u_{-}^{0}\in W(\Omega_{U},\Gamma\times\mathbb{T})$ by Step 1. In the upper part $\Omega_{u}$ one uses $p=0$ by Step 2. By the density of the considered set of test functions in $W(\Omega_{U},\Gamma\times\mathbb{T})$, as $\Omega_{U}$ is Lipschitz, one has that $u=u^{0}$ is the unique solution to problem (10). One concludes that the weak convergences (i) and (ii) have been established for the full sequences. ∎

Remark that Lemma 8.1 holds without hypotheses (H1), (H2), by the same argument, and it takes the following form, which is the result that corresponds to Theorem 7.1 in this case.

Lemma 10.1**.**

Let $u^{\varepsilon}\in H^{1}(\Omega^{\varepsilon},\Gamma)$ be the solutions to (2), and let $u^{0}\in W(\Omega_{U},\Gamma\times\mathbb{T})$ be the solution to (10). Then

[TABLE]

as $\varepsilon$ tends to zero.

We will conclude this section by translating the limit problem (10) into a system in the original coordinates $x\in\Omega$, ($\overline{\Omega}$ is the Hausdorff limit of $\overline{\Omega^{\varepsilon}}$), in a particular case of a finite number of ’bumps’ on the boundary.

Lemma 10.2**.**

Let $\omega$ be an open set in $\mathbb{R}^{N}$. Then there exists an open cover $K^{\prime}$ of $\omega$ such that for all $v\in L^{2}(\omega)$ with $\nabla_{x^{\prime}}v=0$ almost everywhere in $\omega$, $(x,x^{\prime})\in\omega$, one has that $v$ is constant in $x^{\prime}$ almost everywhere in each $V\in K^{\prime}$.

Proof.

A subset $V$ of $\omega$ is said to have the property $(P)$ if the intersection of $V$ and any $k$-plane parallel to the $x^{\prime}$-coordinate $k$-planes is connected, $x^{\prime}\in\mathbb{R}^{k}$. For $(x,x^{\prime})\in\omega$, let $V(x,x^{\prime})$ be a neighborhood of $(x,x^{\prime})$ in $\omega$ that is not contained in any distinct neighborhood of $(x,x^{\prime})$ in $\omega$ with the property (P). By the Zorn lemma, at least one such maximal element exists for the set inclusion partial order on the set of all neighborhoods of $(x,x^{\prime})$ with the property (P), because $\omega$ is open and if $F$ is a totally ordered subset, $\cup_{V\in F}V$ is an upper bound for $F$. Then $K^{\prime}=\{V(x,x^{\prime}):(x,x^{\prime})\in\omega\}$ is an open cover of $\omega$.

Let $v\in L^{2}(\omega)$ be such that $\nabla_{x^{\prime}}v=0$ almost everywhere in $\omega$. Then $v$ has a representative $\tilde{v}$ that is absolutely continuous on almost all line segments parallel to the $x^{\prime}$-coordinate axes and whose classical partial derivatives parallel to the $x^{\prime}$-coordinate axes belong to $L^{2}(\omega)$. It follows from Fubini’s theorem that $\tilde{v}$ assumes constant value on each connected component of each $k$-plane parallel to the $x^{\prime}$-coordinate $k$-planes when intersected with $\omega$. Therefore $v$ is constant in $x^{\prime}$ almost everywhere on each $V\in K^{\prime}$. ∎

Now we restrict to domains such that there is a finite open cover of $\Omega_{U}$ as in Lemma 10.2 with $\omega=\Omega_{U}$ and $x^{\prime}=y$, effectively discarding domains with a countably infinite number of ’bumps’ on the boundary. For instance, excluding $\eta(x,y)=x_{1}^{2}\sin(\frac{y}{x_{1}})$, and $\eta(x,y)={\mathop{\mathrm{dist}}}(y,K)$, where $K$ is a Cantor set.

We say that $K$ is a partition of an open set $\omega$ in $\mathbb{R}^{N}$ if $K$ consists of disjoint nonempty open subsets of $\omega$ such that $\overline{\omega}=\cup_{V\in K}\overline{V}$. Under the assumption that there is a finite open cover $K^{\prime}$ of $\Omega_{U}$ as in Lemma 10.2, we can construct a finite partition $K$ of $\Omega_{U}$ as follows. Given $K^{\prime}=\{V_{1},\ldots,V_{k}\}$, let $K=\{V_{1},V_{2}\setminus\overline{V_{1}},\ldots,V_{k}\setminus\cup_{j=1}^{k-1}\overline{V_{j}}\}$.

Let $K$ be a finite partition of $\Omega_{U}$ with the same property as the open cover in Lemma 10.2. Denote the measure of the connected components of $Y(x)$ by

[TABLE]

for $V\in K$, that is the density of $\Omega^{\varepsilon}$ in $\Omega$. Denote the common boundaries of the subdomains $U,V\in K$ of the partition by

[TABLE]

Let the effective matrix be given by

[TABLE]

Let

[TABLE]

where $\pi_{\mathbb{R}^{2}}$ denotes the projection onto $\mathbb{R}^{2}$, equipped with the natural Hilbert space structure.

The limit problem

[TABLE]

is then well-posed in $W_{K}$.

Corollary 10.1**.**

Let $u^{\varepsilon}\in H^{1}(\Omega^{\varepsilon},\Gamma)$ be the solutions to (2), and let $u^{0}\in W(\Omega_{U},\Gamma\times\mathbb{T})$ be the solution to (10). Let $K$ be a finite partition of $\Omega_{U}$ with the same property as the open cover in Lemma 10.2 with $\omega=\Omega_{U}$ and $x^{\prime}=y$. Then

[TABLE]

as $\varepsilon$ tends to zero. Moreover, $u^{0}=u_{V}^{0}$ in $V\in K$, where $\{u_{V}^{0}\}_{V\in K}$ in $W_{K}$ is the solution to (10).

Proof.

It follows from the weak unfolding convergence $T^{\varepsilon}u^{\varepsilon}\rightharpoonup u^{0}$ in $L^{2}(\Omega_{U})$, that

[TABLE]

as $\varepsilon$ tends to zero. The same computation gives the convergence of the flows. By writing $u^{0}=\sum_{V\in K}\chi_{V}u^{0}$, one verifies that the limit problem (10) is equivalent to the limit problem (10) with $u_{V}^{0}=\chi_{V}u^{0}$ on $V$, for any suitable partition $K$. ∎