arxiv: v1 [math-ph] 22 Apr 2015

Transcription

1 On the Beavers-Joseph-Saffman boundary condition for curved interfaces Sören Dobberschütz March 1, 18 arxiv: v1 [math-ph] Apr 15 Abstract This document is an extended version of the results presented in S. Dobberschütz: Effective Behaviour of a Free Fluid in Contact with a Flow in a Curved Porous Medium, SIAM Journal on Applied Mathematics, 15. The appropriate boundary condition between an unconfined incompressible viscous fluid and a porous medium is given by the law of Beavers and Joseph. The latter has been justified both experimentally and mathematically, using the method of periodic homogenisation. However, all results so far deal only with the case of a planar boundary. In this work, we consider the case of a curved, macroscopically periodic boundary. By using a coordinate transformation, we obtain a description of the flow in a domain with a planar boundary, for which we derive the effective behaviour: The effective velocity is continuous in normal direction. Tangential to the interface, a slip occurs. Additionally, a pressure jump occurs. The magnitude of the slip velocity as well as the jump in pressure can be determined with the help of a generalised boundary layer function. The results indicate the validity of a generalised law of Beavers and Joseph, where the geometry of the interface has an influence on the slip and jump constants. Contents 1 Introduction The Problem on the Microscale 4.1 Description of the Flow using Coordinate Transformations The ε-periodic Problem The effective behavior of a fluid at a curved porous medium 8 4 Derivation of the general law of Beavers and Joseph Correction of the velocity Elimination of the forces Continuity of the traces Second correction of the velocity Correction of the pressure Correction of the divergence Nano-Science Center, University of Copenhagen, Universitetsparken 5, 1 København, Denmark. sdobber@nano.ku.dk 1

2 4..1 Compressibility due to the boundary layer functions Compressibility due to the auxiliary functions Some Estimates for the Corrected Velocity and Pressure Proof of Theorem A Coordinate Transformations of Differential Operators 3 B Very weak solutions of the transformed Stokes system 6 C Various existence theorems 7 C.1 Transformed Stokes equation C. Boundary Layer Functions C..1 The Main Auxiliary Problem C.. Exponential Decay C..3 Application to the Stokes Boundary Layer Problems C..4 Dependence on the parameter x C.3 Functions for the Correction of the Divergence Introduction A now classical result in the theory of homogenization states that, starting with the Stokes or Navier-Stokes equation, the effective fluid flow in a porous medium is given by Darcy s law (see the works of Tartar in [SP8], Allaire in [Hor97] and Mikelić [Mik91]). When dealing with porous bodies inside another fluid, the boundary condition coupling the free fluid flow and the Darcy flow at the porous-liquid interface is of great interest. However, due to the different nature of the governing equations, the derivation of a natural boundary condition is difficult: While the equation for the Darcy velocity consists of a second order equation for the pressure and a first order equation for the velocity, the system of equations governing the free fluid velocity (e.g. the Stokes or Navier-Stokes equation) is of second order for the velocity and of first order for the pressure. For an incompressible fluid, the flow in the direction normal to the interface has to be continuous due to mass conservation. However, additional conditions at the interface are not clearly available. From a mechanical point of view, Beavers and Joseph [BJ67] concluded by practical experiments that a jump in the effective velocity appears in tangential direction. This jump is given by αk 1 (vf v D ) τ = ( v F ν) τ, (1.1) where v F denotes the velocity of the free fluid, v D denotes the effective Darcy velocity in the porous medium and K is the permeability of the porous medium. The factor α is the so-called slip coefficient which has to be determined experimentally. Moreover, ν and τ are the unit normal and tangential vector with respect to the interface separating the porous medium and the free fluid. The Darcy velocity in the absence of outer forces for given fluid viscosity µ is given by v D = 1 µ K p, where p denotes the pressure. Note that the condition mentioned above gives a relation between the velocity of the free fluid at the interface and the effective velocity inside the porous medium it does not impose a condition on the actual fluid velocity inside the porous medium at the interface. Figure 1 contains a schematic illustration.

3 Later, Saffman used a statistical model to derive the boundary condition of Beavers and Joseph. In [Saf71], he argued that v D τ is of lesser order than the other terms and arrived at a jump given by v F τ = 1 α K 1 ( vf ν) τ + O(K). (1.) Other boundary conditions were proposed as well: Ochoa-Tapia and Whitaker for example used the REV-method to obtain that the velocity and pressure as well as the normal stress are continuous over the porous-liquid interface, but a jump appears in the tangential stress in the form ( vd ν v F ν ) τ = βk 1 vd τ. Here v F denotes the averaged free fluid velocity, which is given by a Stokes equation, and v D is the averaged velocity in the porous medium, which in this case fulfills a Darcy law with Brinkman correction, v D = 1 µ K( p µ B v D ). µ B is a known constant, and the dimensionless factor β has to be determined experimentally. For details see [OTW95a] and [OTW95b]. However, a rigorous mathematical derivation of the effective fluid behavior at the boundary was not available until Jäger and Mikelić applied the theory of homogenization to the problem. In [JM96] they developed a mathematical boundary layer together with several corrector terms, which allowed them to justify a jump boundary condition. The main tool was the construction of several boundary layer functions : These functions have a given value at the interface and decay exponentially outside it. They can be used to correct the influence of spurious terms at the boundary, stemming from the contributions of other functions to the fluid velocity and pressure. In [JM], this theory was applied to give a mathematical proof of the Saffman modification of the boundary condition of Beavers and Joseph (see also Section 4 of the Chapter Homogenization Theory and Applications to Filtration through Porous Media in [EFM] for a more comprehensible, simplified version of the proofs), yielding the condition ε ( v F ν) τ = αv F τ + O(ε ) where α = 1 εc D can be calculated explicitely. The constant C D stems from a boundary layer problem for the Stokes equation, cf. [JM]. Numerical simulations of the boundary layer functions can be found in [JMN1]. These results suffer from several drawbacks: First, only a planar boundary in the form of a line or a plane is considered (this also applies to the results of Beavers, Joseph and Saffman). Therefore, the effect of a possible curvature of the interface is not known. Second, the external force on the fluid, appearing as a right hand side in the Navier-Stokes equation, had to be. This issue was adressed in the recent paper [MCM1], together with the derivation of the next corrector for the pressure. Generalizations of the boundary layers in [JM96] were developed by Neuss-Radu in [NR]. However, applications only treat reaction-diffusion systems without flow, and explicit results can only be obtained in the case of a layered medium, see [NR1]. The main problem which makes the treatment of general settings infeasible is the loss of exponential decay of the boundary layer functions (cf. Section C.): With the generalized definition, Neuss-Radu was able to show in [NR] that an exponential stabilization is not possible in a general setting. However, all available tools for the treatment of these problems depend on this type of decay 1. 1 Maria Neuss-Radu, private communications. 3

4 Figure 1: Schematic illustration of the velocity profile for a horizontal flow in a domain consisting of an impermeable upper boundary (with no-slip condition), a free fluid part and a porous region. v F denotes the velocity in the free fluid domain, whereas v D is the effective Darcy velocity. The quantity v = v F Σ v D corresponds to the jump across the interface as discussed in Equation (1.1). In this work, we use the approach developed in [Dob9] and [DB1] for providing a generalization of the law of Beavers and Joseph for curved interfaces. We closely follow [MCM1] for investigating the effective behaviour of a free fluid in contact with flow in a curved porous medium. The main idea is to transform a reference geometry with a straight interface to a domain with a curved interface. It is assumed that the porous part in the reference geometry consists of a periodic array of scaled reference cells and that the flow in the transformed geometry is governed by the stationary Stokes equation. Therefore, one obtains a set of transformed differential equations in the reference configuration. Boundary layer functions for these equations are constructed such that due to the straight transformed boundary their exponential decay can be assured. The difference to [NR] and [NR1] is that in theses works, a periodic geometry was intersected by a curved interface. In this work, a periodic geometry with planar interface is transformed to give the geometry in which the fluid flow takes place. These results have been announced in [Dob14] and [Dob15]. In comparison to the latter paper, we will present the transformed equations in a more general formulation, hopefully facilitating generalizations. The Problem on the Microscale.1 Description of the Flow using Coordinate Transformations In this section we describe the main geometrical setting which is used throughout this work. Let L, K, h >. Then Ω := (, L) ( K, h) is a rectangular domain in R (later corresponding to the transformed domain) with parts Ω 1 := (, L) (, h) (later the transformed free fluid 4

5 Figure : Illustration of the coordinate transformation and the periodic domain. domain), Ω := (, L) ( K, ) (the transformed porous medium) and Σ = (, L) {} (later the transformed interface). Let g C (R) be a given function such that g(x + L) = g(x) for all x R. We consider g to describe a periodic curved structure in our domain of interest. Define the coordinate transformation ψ : Ω Ω ) ( ) ( ) x1 z = 1. z x z + g(z 1 ) ( z1 such that Ω = ψ(ω), Ω 1 := ψ(ω 1 ), Ω := ψ(ω ) and Σ := ψ(σ) = {(x, g(x)) x (, L)}. We are interested in the behavior of a fluid flowing through the curved channel Ω, where Ω 1 represents a domain with a free fluid flow, and Ω is a porous medium. We are especially interested in the behavior of the fluid at the curved boundary Σ. See Figure for an illustration. Let there be given a solid inclusion Ω S Ω. We will later use a sequence of such inclusions to create a porous medium via homogenization theory. For a given volume force f L ( Ω), we assume that a mathematical description of the fluid is given by the steady state Stokes equation with no slip condition on the boundary of the solid inclusion and on the outer walls µ x ũ(x) + x p(x) = f(x) in Ω\ Ω S (.1a) div x (ũ(x)) = in Ω\ Ω S (.1b) ũ(x) = on Ω S Ω\({x 1 = } {x 1 = L}) (.1c) ũ, p are L-periodic in x 1 (.1d) Here µ > denotes the dynamic viscosity; we will set µ = 1 in the sequel. We are looking for a velocity field ũ H 1 ( Ω) and a pressure p L ( Ω)/R. The Stokes equation is an approximation of the full Navier-Stokes equation which is valid for low Reynolds number flows. Using the transformation rules for the differential operators (see Appendix A), we obtain the following equation for the transformed quantities u(z) = ũ(ψ(z)), p(z) = p(ψ(z)) and f(z) = f(ψ(z)) in 5

6 Figure 3: Example of a periodic domain to which our methods can be applied. The upper part corresponds to the free fluid domain, and the lower (gray) part to a porous medium. the rectangular domain Ω: div z (F 1 (z)f T (z) z u(z)) + F T (z) z p(z) = f(z) in Ω\Ω S div z (F 1 (z)u(z)) = in Ω\Ω S u(z) = on Ω S Ω\({z 1 = } {z 1 = L}) u, p are L-periodic in z 1 (.a) (.b) (.c) (.d) Here Ω S := ψ 1 ( Ω S ) is the transformed solid inclusion, and F is defined as the Jacobian matrix of ψ given by [ ] 1 F (z) = g. (.3) (z 1 ) 1 Since det F = 1, ψ is a volume preserving C -coordinate transformation. In this connection, please note that we define the gradient of a vector field column-wise, i.e. u is the transpose of the Jacobian matrix of u. Defining u row-wise leads to slightly different transformed equations. Now, the crucial assumption is that Ω S is given as an ε-periodic structure. This is described in the next subsection.. The ε-periodic Problem We assume an ε-periodic geometry in Ω : Define a reference cell as Y := [, 1], containing a connected open set Y S (corresponding to the solid part of the cell). Its boundary Y S is assumed to be of class C with Y S Y =. Let Y := Y \Y S be the fluid part of the reference cell. For given ε > such that L ε N, let χ be the characteristic function of Y, extended by periodicity to the whole R. Set χ ε (x) := χ( x ε ) and define the fluid part of the porous medium as Ω ε = {x Ω χ ε (x) = 1}. 6

7 Figure 4: The reference cell, consisting of the solid part Y S with boundary Y S, and the fluid part Y. The fluid domain is then given by Ω ε = Ω 1 Σ Ω ε. In order to be able to obtain the effective fluid behavior near Σ, we have to define a number of so-called boundary layer problems, see Section C.. To this end, we introduce the following setting: We consider the domain [, 1] R subdivided as follows: Z + = [, 1] (, ) corresponds to the free fluid region, whereas the union of translated reference cells ( ) Z = {Y }\S k k=1 is considered to be the void space in the porous part. Here S = [, 1] {} denotes the interface between Z + and Z. Finally, let and Z = Z + Z Z BL = Z + S Z be the fluid domain without and with interface. For illustrations, the reader is referred to Figure, Figure 4 and Figure 5. Using the constructions given above, we are interested in the limit behavior of the problem: For ε >, find a velocity u ε and a pressure p ε such that div(f 1 F T u ε ) + F T p ε = f in Ω ε div(f 1 u ε ) = in Ω ε u ε = on Ω ε \({x 1 = } {x 1 = L}) u ε, p ε are L-periodic in x 1 (.4a) (.4b) (.4c) (.4d) Note that here and in the sequel, we use the variable x to designate points in Ω, while y is used for the strip Z or the reference cell Y. There exists a solution (u ε, p ε ) H 1 (Ω ε ) L (Ω ε )/R of the problem above, see Appendix C.1. 7

8 Figure 5: The boundary layer strip Z BL. 3 The effective behavior of a fluid at a curved porous medium The effective velocity field in the free fluid domain Ω 1 is given by: Find a velocity u eff and a pressure p eff such that div(f 1 F T u eff ) + F T p eff = f in Ω 1 (3.1a) div(f 1 u eff ) = in Ω 1 (3.1b) p eff dx = Ω 1 (3.1c) u eff = on (, L) {h} (3.1d) u eff, p eff are L-periodic in x 1 (3.1e) u eff = εc bl on Σ (3.1f) C bl is the decay function of the boundary layer function β bl defined in Section It holds C bl F T e = (see (C.5)). In transformed tangential direction, a slip velocity εc bl F e 1 occurs. Its magnitude can be calculated using the result in Lemma C.3 as u eff (x) F (x)e 1 = εc bl (x) F (x)e 1 = ε The effective Darcy pressure in Ω is given by 1 β bl (x, y 1, +) F (x)e 1 dy 1. (3.) div(f 1 A(f F T p eff )) = in Ω (3.3a) A(f F T p eff ) F T e = on (, L) { K} (3.3b) p eff = p eff + C bl ω on Σ (3.3c) p eff is L-periodic in x 1 (3.3d) 8

9 A is given in (4.5), and Cω bl is the pressure stabilization function defined in Section 4.1., given by C bl ω (x 1 ) = 1 ω bl (x 1, y 1, +) dy 1. Define the effective mass flow rates in transformed tangential direction as M eff := u eff F e 1, dx Ω 1 M ε := u ε F e 1 dx. Ω 1 We obtain the following estimates: 3.1 Theorem. Let f C ( Ω) be L-periodic in the first variable. For (u ε, p ε ) as defined in (.4) and (u eff, p eff ) defined in (3.1) the estimates hold. u ε u eff + M L ε M eff Cε 3 (Ω 1 ) u ε u eff 1 + H (Ω1 p ε p eff ) L + (u ε u eff ) Cε 1 (Ω 1 ) L 1 (Ω 1 ) 4 x 1 (u ε u eff ) + x 1 (p ε p eff ) Cε L (Ω 1 ) 4 In the porous medium Ω, we arrive at the following results: L (Ω 1 ) 3. Theorem. For the effective pressure in the porous medium defined by (3.3), we have for all δ > 3.3 Theorem. At the interface Σ, it holds 1 ε uε A(f F T p eff ) in L ((, L) ( H, δ)) p ε p eff in L (Ω ) p ε p eff H 1 (Σ) Cε 1 (3.4) 1 ε (uε u eff ) in L (Σ) (3.5) u ε u eff H 1 Cε 4 3 (Σ). These results show that the following behavior of a free fluid in contact with a flow in a curved porous medium can be expected for low Reynolds number flows: In the free fluid domain, the velocity and pressure are given by the Stokes equation. In the porous medium, the flow is pressure driven and given by Darcy s law. At the interface, a slip-condition occurs. The velocity is given with the help of the decay function C bl (x) of an auxiliary boundary layer problem. This depends on parameters from the geometry of the interface. 9

10 In transformed normal direction C bl (x) F T (x)e = holds, which is an approximation of continuity of the velocity in that direction. In tangential direction, a jump between the velocities and pressures occurs. Looking at the form of the boundary condition, we see that the generalized boundary condition of Beavers and Joseph has to incorporate effects stemming from the geometry of the fluid-porous interface. The estimates (3.4) and (3.5) corresponds to a generalized rigorous version of (1.). In [JM] and [MCM1], the jump boundary condition is given by u eff 1 + εc1 bl =, whereas in our case it reads u eff + εc bl =. The missing velocity gradient appears in the equations for β bl, see Section Derivation of the general law of Beavers and Joseph In this section we carry out the steps that are necessary to derive the generalized boundary condition of Beavers and Joseph (3.) as well as the theorems given above. We will successively correct the velocity u ε and the pressure p ε given by (.4) with the help of auxiliary functions. This will give us insight into the effective behavior, while at the same time introducing problems with div(f 1 u ε ). Therefore we have to correct this term as well. 4.1 Correction of the velocity Elimination of the forces We start by eliminating the right hand side of (.4a) in Ω 1. Let u, π be a solution of u eff 1 x div(f 1 F T u ) + F T π = f in Ω 1 div(f 1 u ) = in Ω 1 u = on Ω 1 \({x 1 = } {x 1 = L}) u, π are L-periodic in x 1 (4.1a) (4.1b) (4.1c) (4.1d) There exists a unique solution u H 1 (Ω 1 ), π L (Ω 1 /R) by the results for the transfomed Stokes equation. By regularity results (see e.g. [Tem77]), this solution is smooth for smooth f. We extend the velocity u by in Ω and the pressure in a smooth manner to a pressure π defined in all of Ω. As it will turn out, π will be given by the Darcy pressure in Ω. Details of this extension procedure will be given below (see Section 4.4). To obtain estimates for u ε u, we need the following Lemma from [JM96], which is a variant of the Poincaré inequality: 4.1 Lemma. Let φ {g H 1 (Ω ε ) g = on Ω ε \ Ω }. Then φ L (Σ) C ε φ L (Ω ε ) as well as φ L (Ω ε ) Cε φ L (Ω ε ). Define the space W ε = {g H 1 (Ω ε ) g = on Ω ε \({x 1 = } {x 1 = L}), g is L-periodic in x 1 }; then by subtraction the weak formulations, one easily sees that u ε u satisfies the variational equation F T ( u ε u ) : F T φ dx (p ε π ) div(f 1 φ) dx Ω ε Ω ε (4.) = F 1 F T u e φ dσ x F 1 [ π ] Σ φ e dσ x + (f F T π )φ dx Σ Σ 1 Ω ε

11 for all φ W ε. Here [q] Σ := q Ω1 q Ω denotes the jump of the function q across the boundary Σ. 4. Proposition. We have the estimate ε u ε u L (Ω ε ) ε u ε L (Ω ε ) + uε L (Σ) Cε. Proof. Choose φ = u ε u in (4.). Since for u, v H 1 (Ω ε ) the form a(u, v) := F T u : Ω ε F T v dx is bounded and coercive, we obtain u ε u C F T ( u ε u ) dx L (Ω ε ) 4 Ω ε F 1 F T u e (u ε u ) dσ x + F 1 [ π ] Σ (u ε u ) e dσ x Σ Σ + (f F T π )(u ε u ) dx. Ω ε Since u and [ π ] Σ are smooth, the first two terms on the right hand side are bounded by C u ε u L (Σ) C ε u ε u L (Ω ε )4 C(δ)ε + Cδ u ε u L (Ω ε )4. The function (f F T π ) is smooth as well by elliptic regularity results for the Darcy pressure, so the last term is bounded by C u ε u L (Ω ε ) Cε u ε u L (Ω ε )4 C(δ)ε + Cδ u ε u L (Ω ε )4. Choosing δ small enough we arrive at u ε u Cε. Moreover, we have L (Ω ε ) 4 as well as 1 ε uε L (Ω ε ) Cε uε L (Ω ε )4 = Cε u ε u L (Ω ε )4 Cε u ε L (Σ) Cε uε u L (Ω ε )4 Cε. This finishes the proof. In order to derive estimates for the pressure and u ε u in Ω 1, we use the theory of very weak solutions for the transformed Stokes equations, see Appendix B. 4.3 Lemma. It holds ε p ε π L (Ω 1 ) + u ε u L (Ω 1 ) Cε. Proof. u ε u is a very weak solution with G 1 = G =, ξ = u ε ΣT, Σ T = Σ {x = h}. Using the estimate (B.4) from the theory of very weak solutions of the transformed Stokes equation, we obtain u ε u C ξ = C L (Ω 1 ) L uε (ΣT ) L (Σ) Cε by the previous Proposition. Using the Nečas inequality yields ε p ε π C ε( L (Ω 1 ) u ε u + L (Ω 1 ) u ε u ) Cε. L (Ω 1 ) 4 11

12 4.1. Continuity of the traces Looking at the right hand side of equation (4.) and the proof of Proposition 4., one can see that the expression Ω ε (f F T π )φ dx allows estimates for u ε u L (Ω ε ) 4 on the order of ε, whereas the other two integrals only allow for estimates on the order of ε 1. Later we will choose [ π ] Σ to allow for better estimates. This leaves us with the expression Σ F 1 F T u e φ dσ x, which we are going to correct next: Construct the following parameter dependent boundary layer functions (β bl, ω bl ) satisfying div y (F 1 (x)f T (x) y β bl (x, y)) + F T (x) y ω bl (x, y) = div y (F 1 (x)β bl (x, y)) = in Ω Z in Ω Z [β bl (x, y)] S = on Ω S [(F 1 (x)f T (x) y β bl (x, y) F 1 (x)ω bl (x, y))e ] S = F 1 (x)f T (x) u (x)e on Ω S β bl (x, y) = β bl (x, ), ω bl (x, ) are 1-periodic in y 1 on Ω k=1 { Y S ( k) } and define β bl,ε (x) = εβ bl (x, x ε ) as well as ωbl,ε (x) = ω bl (x, x ε ). By the theory of the boundary layer functions (see Appendix C.), there exists decay functions C bl : R R, Cω bl : R R such that 1 β bl,ε εc bl (x)h(x ) + ω L bl,ε C bl ε q (Ω) ω (x)h(x ) + L q (Ω) y β bl,ε L q (Ω) 4 + x (β bl (x, x ε ) H(x )C bl (x)) + div x (F 1 F T y β bl (x, x L q (Ω) 4 ε )) Cε 1 q L q (Ω) Here H denotes the Heaviside function. This correction introduces problems due to the stabilization towards C bl. Therefore we define the following counterflow: div(f 1 (x)f T (x) u σ (x)) + F T (x) π σ (x) = in Ω 1 div(f 1 (x)u σ (x)) = in Ω 1 u σ (x) = on (, L) {h} u σ (x) = C bl (x) u σ, π σ are L-periodic in x 1 Since Σ Cbl (x) F T (x)e dx = (see Lemma C.1), there exists a unique velocity u σ and a pressure π σ, unique up to constants. Consider u (β bl,ε εc bl H(x )) + εu σ H(x ) as a first macroscopic approximation. Define the error between this approximation and the microscopic velocity (similarly for the pressure) as on Σ U ε = u ε u + (β bl,ε εc bl H(x )) + εu σ H(x ), P ε = p ε π H(x ) π H( x ) + (ω bl,ε C bl ω H(x )) + επ σ H(x ). 1

13 Then we have F T U ε : F T φ dx P ε div(f 1 φ) dx Ω ε Ω ε = F 1 [ π ] Σ e φ dσ F 1 Cω bl e φ dσ + (f F T π )φ dx Σ Σ Ω ε + F T x (β bl,ε εc bl H(x )) : F T φ dx div x (F 1 F T y β bl (x, x ))φ dx Ω ε Ω ε ε + ε F 1 F T u σ e φ dσ + F T x (ω bl,ε H(x )Cω bl ) φ dx + ε F 1 π σ e φ dσ Σ Ω ε Σ + F 1 F T x (β bl,ε εc bl H(x ))e φ dσ {x =h} {x = H} F 1 F T x (β bl,ε εc bl H(x ))e dσ. We will later eliminate the first two terms on the right hand side by requiring that [ π ] Σ = Cω bl (denoted as (P), see below). Due to the exponential decay of the boundary layer functions, the last two terms can be chosen arbitrarily small (denoted (BL)). The terms on the right hand side can be estimated by the following orders of ε with respect to U ε (the signs of the terms are kept in order to facilitate the allocation of the terms, and \\ denotes a line break): F T U ε : F T φ dx P ε div(f 1 φ) dx = (P) (P) + O(ε) Ω ε Ω ε + O(ε ) O(ε ) \\ + O(ε ) + O(ε ) + O(ε ) \\ + (BL) \\ (BL). Now U ε is of order ε 1, and U ε is of order ε 3 in Ω ε. For the proof of the main result, we need U ε to be on the order of ε, and U ε should be on the order of ε. We therefore need to correct the velocity again Second correction of the velocity Define the following boundary layer function: div y (F 1 (x)f T (x) y γ bl (x, y)) + F T (x) y µ bl (x, y) div y (F 1 (x)γ bl (x, y)) = = div x (F 1 (x)f T (x) y β bl (x, y)) in Ω Z in Ω Z [γ bl (x, y)] S = on Ω S [(F 1 (x)(f T (x) y γ bl (x, y) µ bl (x, y)))e ] S = γ bl (x, y) = γ bl (x, ), µ bl (x, ) are 1-periodic in y 1 on Ω S on Ω k=1 { Y S ( k) } 13

14 The decay functions are denoted by C bl γ, C bl µ. Define the corresponding counterflow div(f 1 (x)f T (x) b(x)) + F T (x) q(x) = in Ω Correction of the pressure div(f 1 (x)b(x)) = in Ω 1 b(x) = on (, L) {h} b(x) = C bl γ (x) b, q are L-periodic in x 1 on Σ For the correction of the pressure, define the following boundary layer function div y (F 1 (x)f T (x) y λ bl (x, y)) + F T (x) y κ bl (x, y) = F T (x) x (ω bl (x, y) H(x )C bl ω (x)) div y (F 1 (x)λ bl (x, y)) = in Ω Z in Ω Z [λ bl (x, y)] S = on Ω S [(F 1 (x)(f T (x) y λ bl (x, y) κ bl (x, y)))e ] S = λ bl (x, y) = λ bl (x, ), κ bl (x, ) are 1-periodic in y 1 on Ω S on Ω k=1 { Y S ( k) } The decay functions are denoted by C bl λ, Cbl κ. Define the corresponding counterflow div(f 1 (x)f T (x) d(x)) + F T (x) l(x) = in Ω 1 div(f 1 (x)d(x)) = in Ω 1 d(x) = on (, L) {h} d(x) = C bl λ (x) d, l are L-periodic in x 1 on Σ We extend our approximation to u ε and p ε by adjusting U ε and P ε. Define Ũ ε = U ε + (ε γ bl (x, x ε ) ε H(x )C bl γ ) (ε λ bl (x, x ε ) ε H(x )C bl λ ) + ε H(x )b ε H(x )d P ε = P ε + (εµ bl (x, x ε ) εh(x )C bl µ ) (εκ bl (x, x ε ) εh(x )C bl κ ) + ε H(x )q ε H(x )l The weak formulation is now given by F T Ũ ε : F T φ dx Pε div(f 1 φ) dx Ω ε Ω ε = F 1 [ π ] Σ e φ dσ F 1 Cω bl e φ dσ + (f F T π )φ dx Σ Σ Ω ε + F T x (β bl,ε εc bl H(x )) : F T φ dx + ε F 1 π σ e φ dσ Ω ε Σ + ε F 1 F T u σ e φ dσ Σ 14

15 + ε F T x (γ bl (x, x Ω ε ) H(x )Cγ bl ) : F T φ dx ε div x (F 1 F T y γ bl (x, x ))φ dx ε Ω ε ε + ε F T x (λ bl (x, x Ω ε ) H(x )Cλ bl ) : F T φ dx ε div x (F 1 F T y λ bl (x, x ))φ dx ε Ω ε ε ε F 1 F T (b d)e φ dσ + εf T x (µ bl (x, x Σ Ω ε ) H(x )Cµ bl )φ dx ε + εf 1 Cµ bl e φ dσ εf T x (κ bl (x, x Σ Ω ε ) H(x )Cκ bl )φ dx ε εf 1 Cκ bl e φ dσ + ε F 1 (q l)e φ dx + (B). Σ Σ Here (B) designates terms stemming from the boundary layer functions at the outer boundary of Ω. Again, due to their exponential decay, they can be chosen arbitrarily small. Similar to the calculations above, we arrive at the following estimates of Ũ ε in terms of ε: F T Ũ ε : F T φ dx Pε div(f 1 φ) dx = (P) (P) + O(ε) Ω ε Ω ε + O(ε 3 ) + O(ε 3 ) \\ + O(ε 3 ) \\ + O(ε 5 ) O(ε 3 ) \\ + O(ε 5 ) O(ε 3 O(ε 5 ) + O(ε 3 ) \\ + O(ε 3 ) O(ε 3 ) \\ O(ε 3 ) + O(ε 5 ) + (BL). Thus, one arrives at an approximation Ũ ε O(ε), Ũ ε O(ε ) in Ω ε and Ũ ε Σ O(ε 3 ). In order to prove these estimates rigorously, we need to correct the divergence of Ũ ε next. 4. Correction of the divergence 4..1 Compressibility due to the boundary layer functions A calculation shows that div(f 1 Ũ ε )(x) = ε div x [F 1 (β bl (x, x ε ) H(x )C bl (x))] + ε div x [F 1 (γ bl (x, x ε ) H(x )C bl γ (x))] ε div x [F 1 (λ bl (x, x ε ) H(x )C bl λ (x))] We construct the following convergence of the divergence: Let Q β satisfy ) div y (F 1 (x)q β (x, y)) = div x (F 1 (x) [ β bl (x, y) H(y )C bl (x) ]) in Ω Z Q β (x, y) = [Q β ] S (x, y) = C Q β Q β (x, y) is 1-periodic in y 1 on Ω k=1 {Y S ( k) } (x) on Ω S where ( ) C Q β Z (x) = F (x) div x F 1 (x)(β bl (x, y) H(y )C bl (x)) e dy. BL 15

16 Set Q ε β (x) := ε Q β (x, x ε ). Similarly, we define Q γ and Q λ div y (F 1 (x)q γ (x, y)) = div x ( F 1 (x) [ γ bl (x, y) H(y )C bl γ (x) ]) as well as in Ω Z Q γ (x, y) = on Ω k=1 {Y S ( k) } [Q γ ] S (x, y) = Cγ Q (x) on Ω S Q γ (x, y) is 1-periodic in y 1 div y (F 1 (x)q λ (x, y)) = div x ( F 1 (x) [ λ bl (x, y) H(y )C bl λ (x) ]) Q λ (x, y) = [Q λ ] S (x, y) = C Q λ Q λ (x, y) is 1-periodic in y 1 in Ω Z on Ω k=1 {Y S ( k) } (x) on Ω S with decay functions ( ) Cγ Q (x) = F (x) div x F 1 (x)(γ bl (x, y) H(y )Cγ bl (x)) e dy Z BL ( ) C Q λ Z (x) = F (x) div x F 1 (x)(λ bl (x, y) H(y )Cλ bl (x)) e dy. BL Again, define Q ε γ(x) = ε 3 Q γ (x, x ε ) and Qε λ (x) = ε3 Q λ (x, x ε ). The tools and techniques for proving existence and exponential decay for these type of problems can be found in Appendix C.3. The jump across the boundary of the functions above is corrected with the help of the following counterflows: For i {β, γ, λ} define (u Qi, π Qi ) as the solution of div(f 1 (x)f T (x) u Qi (x)) + F T (x) π Qi (x) = in Ω 1 div(f 1 (x)u Qi(x) ) = in Ω 1 u Qi (x) = u Qi (x) = C Q i (x) u Qi, π Qi are L-periodic in x 1 on (, L) {h} on Σ These problems have a solution if Σ CQ i F T e dσ =. For C Q β we obtain ( ) C Q β F T e dσ = div x F 1 (x)(β bl (x, y) H(y )C bl (x)) dy dσ x Σ Σ Z BL ) = div x (F 1 (x)[ β bl (x, y) H(y )C bl (x) dy] dσ x Σ Z BL ( ) = F 1 (x)[ β bl (x, y) H(y )C bl (x) dy] e 1 Σ x 1 Z BL + ( ) F 1 (x)[ β bl (x, y) H(y )C bl (x) dy] e dσ x. x Z BL Now the first term on the right hand side vanishes due to the periodic boundary funcion g, which leads to periodic boundary conditions in x for F, β bl and C bl. The second term vanishes, 16

17 since the latter functions do not depend on x. The argument for the remaining counterflows is analogous. Now define which leads to U ε = Ũ ε Q ε β Q ε γ + Q ε λ + ε u Q β + ε 3 u Qγ ε 3 u Q λ P ε = P ε + ε π Q β + ε 3 π Qγ ε 3 π Q λ, div(f 1 (x)u ε (x)) = ε div x (F 1 (x)q β (x, x ε )) ε3 div x (F 1 (x)q γ (x, x ε ))+ε3 div x (F 1 (x)q λ (x, x ε )) 4.. Compressibility due to the auxiliary functions We correct the divergence even further by defining the functions φ 1,ε i and φ,ε i via and div(f 1 (x)φ 1,ε i (x)) = ε i div x (F 1 (x)q i (x, x ε )) in Ω 1 φ 1,ε i (x) = on (, L) {h} φ 1,ε i (x) = F (x) ε i div x (F 1 (x)q i (x, x Σ Ω 1 ε ))e dx on Σ φ 1,ε i is L-periodic in x 1 div(f 1 (x)φ,ε i (x)) = ε i div x (F 1 (x)q i (x, x ε )) in Ω φ,ε i (x) = on (, L) { K} φ,ε i (x) = F (x) ε i div x (F 1 (x)q i (x, x Σ Ω ε ))e dx on Σ φ,ε i is L-periodic in x 1 Here Q i is extended by into k=1 {Y S ( k) }, i {β, γ, λ}, and ε for i = β ε i := ε 3 for i = γ. ε 3 for i = λ It holds φ 1,ε i H 1 (Ω1 ) C ε i, φ,ε i H 1 (Ω ) C ε i and φ 1,ε i = φ,ε i on Σ. We define the final correction by making use of the following restriction operator (see [Dob9] for details): 4.4 Proposition. Let H = {u H 1 (Ω ) : u = on (, L) { H}, u is L periodic in the first variable} and set H ε = {u H : u = on Ω ε \ Ω }. There exists a linear restriction operator R ε : H Hε such that for w H: and w = on Ω ε \ Ω = R ε w = w Ω ε div(f 1 w) = in Ω = div(f 1 R ε w) = in Ω ε R ε w L (Ω ε ) + ε (R ε w) L (Ω ε ) 4 C( w L (Ω) + ε w L (Ω) 4). 17

18 Using an explicit characterization of R ε (see the Appendix of [SP8] by Luc Tartar), one arrives at the following identity div(f 1 R ε φ) = div(f 1 χ ε(y φ) + +k) ε Y div(f 1 φ) dx. (4.3) k {h Z ; ε(y S +k) Ω } ε(y S +k) Now define U ε = Ũ ε H(x ) P ε = P ε. i {β,γ,λ} φ 1,ε i H( x ) i {β,γ,λ} R ε φ,ε i Since div x (F 1 (x)q i (x, x ε )) = in ε(y S + k), k Z, equation (4.3) leads to div(f 1 U ε ) =, and we can use U ε as a test function. We look at the weak formulation F T U ε : F T φ dx P ε div(f 1 φ) dx. Ω ε Ω ε By calculating the right hand side of this equality and inserting φ = U ε, we obtain (similar to the proof of Proposition 4.) U ε L (Ω ε ) 4 Cε, U ε L (Ω ε ) Cε and U ε L (Σ) Cε 3. Similarly, one arrives at P ε H 1 (Ω ε ) Cε. 4.5 Corollary. Let π be a smooth function satisfying [ π ] Σ = C bl ω. Then the estimates hold. U ε H 1 (Ω1 ) + ε 1 P ε L (Ω1 ) Cε 3 Proof. The result follows as in [MCM1] by transforming the equations in Ω 1 back to a Stokes system in Ω Some Estimates for the Corrected Velocity and Pressure Note that U ε differs from U ε by O(ε ) in the L -norm, and by O(ε) in the H 1 -norm. This yields the estimates ε P ε H 1 (Ω ε ) + ε U ε L (Ω ε ) 4 + U ε L (Ω ε ) + ε 1 U ε L (Σ) Cε as well as ε P ε L (Ω1 ) + U ε H 1 (Ω1 ) Cε 3. We now collect some more estimates. They are analogous to what is known in the literature [MCM1], and the proofs follow along the same lines. 4.6 Lemma. u ε u eff L (Ω 1 ) Cε 3 18

19 Proof. Note that u ε u eff = U ε (β bl,ε εh(x )C bl ). By the theory of very weak solutions we have u ε u eff L (Ω 1 ) C U ε (β bl,ε εh(x )C bl ) L (Σ) Cε Lemma. ε (u ε u eff ) + (u ε u eff ) Cε, L (Ω 1 ) L 1 (Ω 1 ) Proof. First note that (u ε u eff ) L (Ω 1 ) U ε L (Ω1 ) + (β bl,ε εh(x )C bl ) L (Ω 1 ) Cε + Cε 1 Cε 1. The second assertion follows via (u ε u eff ) L 1 (Ω 1 ) U ε L (Ω1 ) + (β bl,ε εh(x )C bl ) L 1 (Ω 1 ) Cε. 4.8 Lemma. Proof. By interpolation we obtain u ε u eff 4.9 Lemma. 1 H (Ω 1 ) C u ε u eff 1 u ε u eff 1 Cε H (Ω1 ) L (Ω 1 ) (u ε u eff ) 1 L (Ω 1 ) Cε 3 1 x 1 (u ε u eff ) + x 1 (p ε p eff ) L (Ω 1 ) L (Ω 1 ) Cε 1 ε 1 Cε. Proof. We have x 1 (u ε u eff ) Cε + x 1 U ε + x 1 (β bl,ε εh(x )C bl ) L (Ω 1 ) L (Ω 1 ) Cε + C U ε L (Σ) + x 1 (β bl,ε εh(x )C bl ) L Cε. 1 (Ω 1 ) L (Ω 1 ) For the second part, use Proposition 16 in [MCM1] to obtain x 1 (p ε p eff ) C u ε u eff + Cε Cε L (Σ) L (Ω 1 ) L (Ω 1 ) 4.1 Lemma. Proof. p ε p eff L 1 (Ω 1 ) Cε p ε p eff L 1 (Ω 1 ) P ε (ω bl H(x )Cω bl ) εh(x )π σ L 1 (Ω 1 ) P ε + L (Ω1 ) ω bl H(x )C bl + Cε Cε ω L 1 (Ω 1 ) 4.11 Lemma. M ε M eff Cε 3 19

20 Proof. M ε M eff = (u ε u eff ) F e 1 dx Ω 1 C u ε u eff L (Ω 1 ) Cε Lemma. 1 ε (uε u eff ) in L (Σ) u ε u eff H 1 Cε 9 5 (Σ). Proof. It holds u ε u eff = U ε (β bl,ε εc bl ) on Σ. Since ε 1 U ε L (Σ) Cε and βbl (, ε ) 1 βbl (, y 1, +) dy 1 = C bl ( ) in L (Σ) (cf. [Amo97] or the remarks in [All9]), we obtain the first result. Now let v H 8 5 # (Σ). To begin with, we have β bl (, ε ) Cbl, v H 8 5 # (Σ), H 8 5# (Σ) 1 ε = 1 (β bl (x 1, +, x 1 ε, +) Cbl (x 1, +))v(x 1 ) dx 1 = ε (β bl (εy 1, +, y 1, +) C bl (εy 1, +))v(εy 1 ) dy 1 = ε (β bl (εy 1, +, y 1, +) C bl (εy 1, +))v(εy 1 ) dy 1, i I ε Y i where I ε N is some finite index set depending on ε, and Y i is an interval of the form [k, k + 1) for some k N. For each i I ε, choose a yi Y i. Since β bl (yi, +,, +) Cbl (, +) has mean value over Σ, we have that (β bl (εy Y i i, +, y 1, +) C bl (εyi, +)) v(εy i ) dy 1 =. Thus the expression above is equal to ε (β bl (εy 1, +, y 1, +) C bl (εy 1, +))v(εy 1 ) i I ε By Taylor series expansion, there exists a ξi Y i ε i I ε Y i = ε i I ε Y i (β bl (εy i, +, y 1, +) C bl (εy i, +))v(εy i ) dy 1. [ (β bl (εξi, +, y 1, +) C bl (εξi, +))v(εξi ) x 1 Y i such that the above expression is equal to ] ε(y 1 yi ) dy 1 [ ] (β bl (εξi, +, y 1, +) C bl (εξi, +))v(εξi ) εy 1 dy 1 x 1 by the vanishing mean value of β bl C bl. Now transformation of the integral and rearrangement yields ε i I ε εy i [ x 1 [β bl (ξ i, +, x 1 ε, +) Cbl (ξ i, +)]v(ξ i ) +v (ξ i )β bl (ξ i, +, x 1 ε, +) Cbl (ξ i, +) ] x 1 dx 1

21 By the embedding H 8 5 (R) C 1 (R) as well as by the boundedness of the terms containing β bl C bl and x 1, the integral can be estimated by εc v 8 H 5 (εyi, which upon summation gives ) a bound of εc v 8 H 5 (Σ) for the whole term. This implies that β bl,ε εc bl H 5 8 Cε. (Σ) Since L (Σ) H 8 5 (Σ), we have that U ε 8 H Cε. 5 (Σ) u ε u eff H 5 8 Cε, and by interpolation we obtain (Σ) u ε u eff H 1 (Σ) C u ε u eff Proof of Theorem 3. L (Σ) u ε u eff 4 5 H 3 (Σ) Thus similar to the first part Cε 1 5 ε 4 5 = Cε 9 5. We need the following auxiliary constructions: Fix x Ω and let (w j, π j ) H 1 # (Y ) L (Y )/R be a solution of the parameter dependent cell problem div y (F 1 (x)f T (x) y w j (x, y)) + F T (x) y π j (x, y) = e j div y (F 1 (x)w j (x, y)) = w j (x, y) = w j (x, y), π j (x, y) are Y -periodic in y. in Y in Y in Y S (4.4a) (4.4b) (4.4c) (4.4d) Define the matrix A by [A(x)] ji = Y w j i (x, y) dy i, j = 1,. (4.5) One can show that A is symmetric and positive definite. The Darcy pressure is given by the following problem: Find π H 1 (Ω ) such that div(f 1 A(f F T π )) = in Ω (4.6a) A(f F T π ) F T e = on (, L) { K} (4.6b) π = π + Cω bl on Σ (4.6c) π is L-periodic in x 1 (4.6d) Here π is given in (4.1). Using the weak formulation for Ũ ε, P ε, one arrives at the estimates Ũ ε L (Ω ε ) Cε, Ũ ε L (Ω ε ) 4 Cε, P ε L (Ω ε ) C. By the theory of two scale convergence (see [Ngu89], [All9]) there exist limits U lim L (Ω ; H# 1 (Y )) as well as P lim L (Ω ; L (Y )) with Ũ ε ε U lim, Ũ ε ε y U lim, Pε P lim. Choose a test function φ C (Ω ; H# 1 (Y F )) with div y (F 1 (x)φ) = and test the weak formulation for Ũ ε. In the -scale limit, one obtains F T (x) y U lim (x, y) : F T (x) y φ(x, y) dy dx Ω Y F P lim (x, y) div x (F 1 (x)φ(x, y)) dy dx = (f(x) F T (x) π (x))φ(x, y) dy dx. Ω Y F Ω Y F (4.7) 1

22 By choosing a test function φ C (Ω ; H# 1 (Y F )), testing with εφ(x, x ε ) and taking the -scale limit, one arrives at Ω Y F P lim (x, y) div y (F 1 (x)φ(x, y)) =, thus P lim does not depend on the variable y. Upon a separation of scales similar to the derivation of Darcy s law (see [Hor97] and [Dob9] for the case of the transformed equations), one obtains the representation [ ( )] U lim (x, y) = w j (x, y) f(x) F T (x) π (x) + P lim (x) j=1 [ ( j )] U lim (x, y) dy = A(x) f(x) F T (x) Y f π (x) + P lim (x). Choosing Ũ ε Ω1 as a test function yields Ũ ε L (Ω1 ) 4 Cε 3 and thus Ũ ε ε in L (Ω 1 ). Choose a φ C (Ω; H# 1 (Y F )) such that div y (F 1 (x)φ) =. Similar to the above calculations, testing with φ(x, x ε ) gives in the limit F T (x) y U lim (x, y) : F T (x) y φ(x, y) dy dx Ω Y F P lim (x) div x (F 1 (x)φ(x, y)) dy dx = (f(x) F T (x) π (x))φ(x, y) dy dx. Ω Y F Y F Using (4.7), one arrives at Ω 1 Y F P lim (x) div x (F 1 (x)φ(x, y)) dy dx =. Upon an integration by parts, this is equivalent to F T (x) P lim (x)φ(x, y) dy dx + P lim (x)φ(x, y) F T (x)e =. Y F Σ Y F Ω 1 Since φ is arbitrary, we conclude that F T P lim = in Ω 1 and P lim Σ =. Thus P lim = and Y f U lim (x, y) dy = A(x)(f(x) F T (x) π (x)), which is Darcy s law. It remains to compare π with p eff. Note that p eff π H k 1 (Ω ) C p eff π H k 1 (Σ) C p eff π H k (Ω 1 ). After transformation to Ω 1, the functions u eff, p eff, u and p satisfy a Stokes equation. By standard regularity results for this equation (see [Tem77], [MCM1]) we obtain p eff π H k (Ω 1 ) Cε for all k N. This gives the first statement of the theorem. Additionally, since Ω P ε Ω = p ε π + O(ε 1 ) = p ε p eff + O(ε 1 ), we have that p ε p eff in L (Ω ), which is the second assertion. Finally, we estimate p ε π H 1 (Σ). Let w H 1 (Ω ε ) such that w = on Ω ε \({x 1 = } {x 1 = L} Σ) and w is L-periodic in x 1. Then by trace estimates [Tem77] p ε π, w Σ H 1, H 1 ((pε π )e, w) L (Ω ε ) + ( x (p ε π ), w) L (Ω ε ) For the first term on the right hand side, notice that p ε π = (p ε p eff ) + ( p eff π ) = O(ε 1 ) + O(ε). For the second term, observe that ( x (p ε π ), w) L (Ω ε ) Cε w H 1 (Ω ε ) + Pε H 1 (Ω ε ) w H 1 (Ω ε ) + ( (ωbl,ε H(x )C bl ω ), w) L (Ω ε ) Cε w H 1 (Ω ε ) + Cε 1 ε w L (Ω ε ). This shows that p ε π H 1 (Σ) C ε, which implies the last statement of the theorem.

23 A Coordinate Transformations of Differential Operators We recall the definition of coordinate transformations and some differential operators and investigate their relations: Let Ω R n with n N be a Lipschitz domain; let c : Ω R be a scalar function, j : Ω R n a vector field and M : Ω R n n a matrix function. They are assumed to be sufficiently smooth. A.1 Definition. The gradient of a vector field is defined as ( j) ik = j k x i, for i, k = 1,..., n (i.e. j is the transpose of the Jacobian matrix of j); the divergence of a matrix-valued function is defined column-wise, thus (div( M)) k = n i=1 M ik x i, for k = 1,..., n; and the Laplacian of a vector field is given by For n = we define the two operators and Curl( c) = j = div( j). ( c ) x c = x 1 [ ] 1 c 1 curl( j) = j x 1 j 1 x. The curl -operators above are two-dimensional variants of the well-known curl operator describing the rotation of three-dimensional vector fields. We have the relations curl c = and div Curl c =, and curl is the formal adjoint of Curl (see [Ver7] and [DL9] for details concerning these operators). We call ψ a regular orientation- A. Definition. Let Ω, Ω R n be Lipschitz domains and let ψ : Ω Ω. preserving C k -coordinate transformation if 1. ψ is a C k -diffeomorphism, and. There exist c, C > such that where F denotes the Jacobian matrix of ψ. If det F 1, we call ψ volume preserving. c det F (z) C z Ω, 3

24 We will indicate coordinates in Ω by z = (z 1,..., z n ) and those in Ω by x = (x 1,..., x n ). Define c(z) := c(ψ(z)) j(z) := j(ψ(z)) M(z) := M(ψ(z)). A.3 Lemma. Let ψ : Ω Ω be a C 1 -coordinate transformation. Denote by F the Jacobian matrix of ψ, and let J(z) := det(f (z)). With the notations and definitions above it holds 1. F T z c = x c.. div z (JF 1 j) = (J ψ 1 ) div x ( j). 3. div z (JF 1 M) = (J ψ 1 ) div x ( M). Proof. The first assertion is a simple application of the chain rule, whereas the second one is known as the Piola-transformation (see [Zei88], Chapter 61. Note that Zeidler defines vectors and gradients row-wise, leading to slightly different formulas.) For the matrix divergence the second statement holds column-wise. Application of this lemma yields: A.4 Lemma. Let ψ be a volume-preserving C 1 -coordinate transformation. The operators from Definition A.1 transform according to 1. x ( c) = div z (F 1 F T z c).. x ( j) = div z (F 1 F T z j). 3. div x ( j) = div z (F 1 j). 4. div x ( M) = div z (F 1 M). 5. Curl x ( c) = Curl z (c), with 6. curl x ( j) = curl z (j), with Curl z (c) = [ ] 1 F T 1 z c. curl z (j) = curl z (F T j). Proof. For volume-preserving coordinate transformations it holds J 1, thus in that case by the preceding lemma we have div z (F 1 j) = div x ( j) and div z (F 1 M) = div x ( M), which gives the third and the fourth statement. The first and the second statement follow by the equalities x ( c) = div x ( x c) and x ( j) = div x ( x j) and application of the above results to the right hand sides. The fifth statement follows along the same lines, whereas the sixth can be obtained by a direct calculation of the effect of the transformation on the defining equation. 4

25 A simple computation shows that [ ] 1 F 1 T = F Curl z (c) = F Curl z (c). [ ] 1 ; thus it holds 1 A.5 Lemma (Transformed Differential Identities). Let ψ be a volume-preserving C 1 -coordinate transformation as above. Then the following identities hold: 1. div z (F 1 c) = F T z c.. div z (F 1 F T z (div z (F 1 j))) = div z (F 1 div z (F 1 F T z j)). 3. div z (F 1 Curl z (c)) =. 4. div z (F 1 (cj)) = c div z (F 1 j) + F T c j. 5. curl z (F T z c) =. 6. curl z (div z (F 1 F T z j)) = div z (F 1 F T z curlz (j)). 7. If div z (F 1 j) =, then F T z ( curl z (j)) = [ ] 1 div 1 z (F 1 F T z j). 8. F T z (div z (F 1 F T c)) = div z (F 1 F T z (F T z c)). Proof. To obtain the first statement transform the well-known equation div x ( ci) = x c. The second follows from x (div x j) = div x ( x j). Next transform div x (Curl x ( c)) = and div x ( c j) = c div x ( j) + x c j. Finally observe that curl x ( x c) = as well as curl x ( x j) = x (curl j). If div x ( j) =, a simple calculation together with the fact that in this case j 1 x 1 = j x shows that x (curl x j) = [ ] 1 1 x ( j), which upon transformation yields the result. For the last statement consider x ( x c) = x ( x c). A.6 Remark. Let ν(x) be the unit normal vector at x Ω. Then the corresponding transformed unit normal vector is given by ν(x) = F T (x)ν(x) 1 F T (x)ν(x). If n =, the unit tangential vector τ(x) has the direction [ ] 1 1 F T (x)ν(x) = F (x) [ ] 1 1 ν(x) = F (x)τ(x), thus it holds indicates the chosen norm in R n. τ(x) = F (x)τ(x) 1 F (x)τ(x). 5

26 B Very weak solutions of the transformed Stokes system In this appendix we develop the theory of very weak solutions for the transformed Stokes equations in Ω 1, which has been suggested for F = Id by Conca in [Con87]. We then derive estimates of the velocity and the pressure. We are looking for functions (B, β) L (Ω 1 ) H 1 (Ω 1 ) as solution of div(f 1 F T B) + F T β = G 1 + div(f 1 G ) in Ω 1 div(f 1 B) = in Ω 1 B = ξ on Σ T B, β are L-periodic in x 1. (B.1a) (B.1b) (B.1c) (B.1d) Here Σ T := Σ (, L) {h}, and ξ L (Σ T ), G 1 L (Ω 1 ) as well as G L (Ω 1 ) 4 are given functions. We require the compatibility condition F 1 ξ e dσ + F 1 ξ e dσ =. Σ {x =h} Define the space W 3 = {g H 1 (Ω 1 ), g is L-periodic in x 1, Ω 1 g dx = } and let g L (Ω 1 ), u W 3. The auxiliary problem div(f 1 F T φ) + F T π = g in Ω 1 div(f 1 φ) = u in Ω 1 φ = on Σ T φ, π are L-periodic in x 1. (B.a) (B.b) (B.c) (B.d) can be transformed to a Stokes system in Ω 1 and solved on that domain. This yields the existence of a unique solution (φ, π) H (Ω 1 ) H 1 (Ω 1 )/R with φ + π C( g + u H (Ω1 ) L (Ω 1 ) L (Ω 1 ) H 1 (Ω 1 ). (B.3) ) Using integration by parts, one calculates that B g dx β, u H 1, H 1 = Ω 1 (G 1 + div(f 1 G ))φ dx Ω 1 F 1 (F T φ πi)ν ξ dσ Σ T Define for (g, u) L (Ω 1 ) W 3 l(g, u) = (G 1 + div(f 1 G ), φ) Ω1 (F 1 (F T φ πi)ν, ξ) ΣT, where φ, π is a solution of (B.). B.1 Definition. (B, β) L (Ω 1 ) H 1 (Ω 1 ) is called a very weak solution of the Stokes problem (B.1), if for all (g, u) L (Ω 1 ) W 3 the identity holds. (B, g) Ω1 β, u H 1, H 1 = l(g, u) 6

27 B. Lemma. The functional l : L (Ω 1 ) W 3 R is linear and continuous. Proof. The solution operator to (B.) is linear, as is the gradient and the duality pairing. Thus, one obtains the linearity of l. For the continuity, we estimate l(g, u) (G 1 + div(f 1 G ), φ) Ω1 + (F 1 (F T φ πi)ν, ξ) ΣT G 1 L (Ω1 ) φ L (Ω 1 ) + G L (Ω1 ) 4 φ H 1 (Ω 1 ) + C φ L (Σ T ) 4 ξ L (Σ T ) + C π L (Σ T ) ξ L (Σ T ) ( G 1 L (Ω1 ) + G L (Ω1 ) 4) φ H (Ω 1 ) + C φ H 3 (Ω1 ) ξ L (Σ T ) + C π L (Ω 1 ) ξ L (Σ T ) C( G 1 L (Ω1 ) + G L (Ω1 ) 4 + ξ L (Σ T ) )( φ H (Ω 1 ) + π L (Ω 1 ) ) C( G 1 L (Ω1 ) + G L (Ω1 ) 4 + ξ L (Σ T ) )( g L (Ω 1 ) + u H 1 (Ω 1 ) ) by (B.3). This shows that l is bounded and thus continuous. B.3 Proposition. There exists a unique very weak solution of (B.1). Proof. The Lemma above shows that l (L (Ω ) W 3 ). Since L (Ω ) W 3 is a Hilbert space, the Riesz representation theorem yields the existence of a unique B L (Ω 1 ) and a unique β W 3 such that ( B, g) Ω1 + β, u H 1, H 1 = l(g, u) g, u L (Ω 1 ) W 3. By choosing B = B, β = β, we see that (B, β) is a very weak solution of (B.1). B.4 Lemma. Let (B, β) be a very weak solution of (B.1). There exists a constant C > such that the estimate holds. B C( G L (Ω1 ) 1 + G L (Ω1 ) + ξ L (Ω1 ) 4 L (Σ T ) ) (B.4) Proof. Choose g = B and u =. Using the estimates from the proof of the previous lemma together with the scaled Young s inequality yields the result. C Various existence theorems C.1 Transformed Stokes equation In this section we prove the existence and uniqueness of the solution of Problem (.4) for fixed ε. Basically, the approach is the same as in the functional-analytic treatment of the Stokes equation (see e.g. [SP8]). By multiplying (.4a) with φ H 1 div (Ωε ) where H 1 div(ω ε ) := {w H 1,#(Ω ε ) div(f 1 w) = }, integrating by parts and noting that F T (x) p(x) φ(x) dx = p(x) div(f 1 (x)φ(x)) dx =, Ω ε Ω ε 7

28 we obtain the weak formulation of Problem (.4) in the form F T (x) u ε (x) : F T (x) φ(x) dx = Ω ε f(x) φ(x) dx Ω ε φ Hdiv(Ω 1 ε ) (C.1) Note that Hdiv 1 (Ωε ) is a Banach space with respect to the norm L (Ω ε ). following lemma for the estimation of the left hand side: We need the C.1 Lemma. There exist constants < k F < K F such that for the eigenvalues λ(x) of F 1 (x)f T (x) holds k F < λ(x) < K F x Ω, i.e. F 1 (x)f T (x) is symmetric and positive definite. Proof. A calculation of the eigenvalues λ 1 (x), λ (x) of F 1 (x)f T (x) yields λ 1 (x) = 1 + g (x 1 ) + g (x 1 ) + g (x 1 ) 4, 4 λ (x) = 1 + g (x 1 ) g (x 1 ) + g (x 1 ) 4. 4 Because of the smoothness of g there exists an M > 1 such that g (x) < M for all x Ω. Obviously λ i (x) 1 + M =: K F, i = 1,. Choose a k F small enough such that M + 1 k F. Another calculation shows that λ (x) k F g (x 1 ) 1 k F + k F, which gives the desired result since λ 1 (x) λ (x). C. Proposition. Let ε > be fixed and let F be given by (.3). For given f L (Ω), the Problem (C.1) has a unique solution u ε H 1 div (Ωε ). Proof. Define for u, v Hdiv 1 (Ωε ) the (bi-)linear forms a(u, v) = F T (x) u(x) : F T (x) v(x) dx Ω ε and b(v) = f(x) v(x) dx. Ω ε The continuity of b for f L (Ω) is standard. In order to apply the lemma of Lax-Milgram, we have to show that a is continuous and coercive. First note that as a pointwise estimate we have F T (x) v(x) : F T (x) v(x) = tr( v(x) T F 1 (x)f T (x) v(x)) = e T i v(x) T F 1 (x)f T (x) v(x)e i i=1 v(x)e i F 1 (x)f T (x) v(x)e i i=1 K F e T i v(x) T v(x)e i i=1 = K F v(x) : v(x) = K F v(x), 8

29 with being the Euclidean vector- and matrixnorm. This gives the continuity of a due to F T (x) u(x) : F T (x) v(x) dx Ω ε ( = F T (x) u(x) ) ( ij F T (x) v(x) ) ij dx Ω ε i,j=1 Ω ε ( ( F T ) ) 1 ( (x) u(x)) ij i,j=1 i,j=1 ( F T (x) v(x)) ij ) ) 1 ( = F T (x) u(x) : F T (x) u(x) ) 1 ( F T (x) v(x) : F T (x) v(x) ) 1 dx Ω ε K F Ω ε u(x) v(x) dx K F u L (Ω ε ) v L (Ω ε ), where the Cauchy-Schwarz inequality in L has been used in the last step. For the coercivity consider k F v L (Ω ε ) λ (x) v(x) : v(x) dx Ω ε λ i (x)e T i v(x) T v(x)e i dx Ω ε i=1 F T (x) v(x) : F T (x) v(x) dx. Ω ε Now the Lax-Milgram lemma implies the proposed result. Due to (C.1), the solution u ε fullfills div(f 1 F T u ε ) f (H 1 div(ω ε ) ). We will now characterize the orthogonal complement (H 1 div (Ωε ) ) of H 1 div (Ωε ) in order to reintroduce the pressure. We remind the reader of the following results, the proofs of which can be found in [Wlo9]: C.3 Theorem (Generalized Trace Theorem). Let k R, k, and let Λ be a bounded domain in R n, n N with boundary Λ C k+1. There exists a continuous linear operator T : H k (Λ) H k 1 ( Λ) H k 1 1 ( Λ) with dx T (φ) = (φ Λ, φ ν Λ) for all φ C k ( Λ). C.4 Theorem (Generalized Inverse Trace Theorem). Let Λ be a bounded domain in R n, n N, with boundary Λ C k+1 with a given k R, k. Let T be defined as above. There exists a continuous linear extension operator E : H k 1 ( Λ) H k 1 1 ( Λ) H k (Λ) such that T E = Id. 9