arxiv: v1 [math.ap] 13 Nov 2013

Transcription

1 CORRECTORS FOR THE NEUMANN PROBLEM IN THIN DOMAINS WITH LOCALLY PERIODIC OSCILLATORY STRUCTURE arxiv: v1 [math.ap] 13 Nov 13 MARCONE C. PEREIRA AND RICARDO P. SILVA Abstract. In this paper we are concerned with convergence of solutions of the Poisson equation with Neumann boundary conditions in a two-dimensional thin domain exhibiting highly oscillatory behavior in part of its boundary. We deal with the resonant case in which the height, amplitude and period of the oscillations are all of the same order which is given by a small parameter >. Applying an appropriate corrector approach we get strong convergence when we replace the original solutions by a kind of first-order expansion through the Multiple-Scale Method. 1. Introduction Given a small positive parameter, our goal in this work is to analyze the convergence of the family of solutions w H 1 R ) of the Neumann boundary problem w + w = f in R, w ν = on 1.1) R, where f L R ) is uniformly bounded which respect to, ν = ν1, ν) denotes the unit outward normal vector field to R and ν denotes the outside normal derivative. The domain R R is a thin domain with a highly oscillating boundary which is given by R = {x 1, x ) R : x 1, 1), bx 1 ) < x < G x 1 )}, 1.) where b and G are smooth functions, uniformly bounded by positive constants b 1, G, and G 1, independent of, satisfying bx) b 1 and < G G x) G 1 for all x [, 1]. Noticing that R, 1) b 1, G 1 ) for all >, we see that the domain R is actually thin on the x -direction. Therefore we will say that R is collapsing to the interval, 1) as. Here, the function b, independent on, establishes the lower boundary of the thin domain, and the function G, dependent on the parameter, defines its upper boundary. Also, allowing G to present oscillatory behavior we obtain an oscillating thin domain R, whose amplitude and period of the oscillations are of the same order, which also coincides with the order of its thickness. For instance, we may think the function G of the form G x) = mx) + kx)gx/), where m, k : [, 1] R are C 1 -functions and g : R R is a L-periodic and smooth function. In this case, it is clear that the functions m and b define the profile of the thin domain, as well as k the amplitude of the oscillations which is given by the periodic function g. This includes the case where the function G is a purely periodic function, as G x) = + sinx/), but also includes the case where the function G exhibit a profile and amplitude of oscillations dependent on x, 1). We can refer the reader to [16] for an example of this kind of locally periodic structure. In fact, we are able to treat even more general cases, see hypotheses H1) and H) in the Section. Figure 1 below illustrates an example of such a thin domain that we are able to treat in this work. 1 Mathematics Subject Classification. 35B5, 35B7, 74Kxx. Key words and phrases. thin domains, correctors, boundary oscillation, homogenization. Partially supported by CNPq 351/8-4 and FAPESP 8/ and 1/1879-, Brazil. Partially supported by FUNDUNESP and FAPESP 8/ and 1/6753-8, Brazil. 1

2 M. C. PEREIRA AND R. P. SILVA R Ε Figure 1. An oscillating thin domain R. Notice that the existence and uniqueness of solutions for problem 1.1) is guaranteed by Lax-Milgram Theorem. The aim of this work is to treat the convergence aspects of the solution w as, to an effective limit solution as the domain gets thinner and thinner although the oscillations becomes higher. In [4, 5, 6], the authors combine methods from homogenization theory and boundary perturbation problems for partial differential equations to obtain the homogenized effective) equation associated to 1.1). Since the thin domain R degenerates to the interval, 1) as goes to zero, an one-dimensional limiting equation is established. The respective solution, w H 1, 1), is the limit of the family w H 1 R ). We also have that the homogenized coefficient of the limit equation depends on a family of auxiliary solutions, capturing somehow, the oscillatory geometry of the thin domain. We refer the reader to [9, 18, 31, 3] for a general introduction to the theory of homogenization and to [3] for a general treatise on boundary perturbation problems for partial differential equations. It is important to notice that the convergence of the sequence w H 1 R ) to the solution w H 1, 1) of the homogenized equation established in [6, Theorem.3] is only weak with respect to H 1 -norms. In fact, weak convergence in H 1 -norms is the best one for these kind of problems exhibiting oscillatory behavior since in [, 9, 17], the authors show that in general strong convergence is not possible. For this reason, seeking for better convergence results, we apply an appropriate corrector approach developed by Bensoussan, Lions and Papanicolaou in [9], to obtain a sort of strong convergence. The ideia is adjust the homogenized solution w, adding suitable corrector functions κ H 1 R ), κ = o) in L R ), obtaining 1/ w w κ H1 R ), as. Here we proceed as in [, 9, 3] rescaling the Lebesgue measure of H 1 R ) by a factor 1 in order to preserve the relative capacity of mensurable subset O R for small values of. Our challenge here is to set a proper family of corrector functions in the way to get the desired strong convergence once the auxiliary solutions associated with the homogenized limiting) equation are posed in an one-parameter family of representative cells. See Theorem 3. for a complete statement on these results. In such problems, where the analytical solution demands a substantial effort or is unknown, it is of interest in discussing error estimates in norm) when we replace it by numerical approximations. In this sense, the analysis performed here is fundamental in order understanding the convergence properties of the solutions of the problem 1.1). There are several works in the literature dealing with boundary value problems posed in thin domains presenting oscillating behavior. Among others, we mention [5, 6, 8] where asymptotic approximations of solutions to elliptic problems in thin perforated domains defined by purely periodic functions with rapidly varying thickness were studied, as well as [1, 11], dealing with nonlinear monotone problems in a multidomain with a highly oscillating boundary. Let us also cite [1, 8] where the asymptotic description of nonlinearly elastic thin films with fast-oscillating profile was successfully obtained in a context of Γ-convergence [19]. We also can mention other works addressing the problem of the behavior of solutions of partial differential equations where just part of the boundary of the domain presents oscillation. This is the case in [15] where the authors deal with the Poisson equation with Robin type boundary conditions in the oscillating part of the boundary. In general, the differential equation is not affected by the perturbation but the boundary

3 CORRECTORS FOR THE NEUMANN PROBLEM IN THIN DOMAINS 3 condition. Some very interesting interplay among the geometry of the oscillations and Robin boundary conditions can be analyzed and depending on this balance, the limiting boundary condition may be of a different type. Similar discussion can be found in [14] for the Navier-Stokes equations with slip boundary conditions and in [3] for nonlinear elliptic equations with nonlinear boundary conditions as well the references given there. In the case of the present paper, the effect of the oscillations at the boundary is coupled with the effect of the domain being thin. As previously observed, the oscillatory behavior considered here affect completely the limit behavior of the equation. Let us also point out that thin structures with rough contours as thin rods, plates or shells, as well as fluids filling out thin domains in lubrication procedures, or even chemical diffusion process in the presence of grainy narrow strips are very common in engineering and applied sciences. The analysis of the properties of these processes taking place on them, and understanding how the micro geometry of the thin structure affects the macro properties of the material is a very relevant issue in engineering and material design. Therefore, being able to obtain the limiting equation of a prototype equation like the Poisson equation 1.1) in different structures where the micro geometry is not necessarily smooth and being able to analyze how the different micro scales affects the limiting problem, goes in this direction and will allow the research and deep understanding in more complicated situations. In this context we also mention [1, 13, 7] for some concrete applied problems. We also would like to observe that although we just treat the Poisson equation with Neumann boundary conditions, we also may consider other different conditions on the lateral boundaries of the thin domain while preserving the Neumann type boundary condition in the upper and lower boundary. Indeed, the limiting problem preserves the boundary condition. The structure of this paper is as follows. In Section, we establish the functional setting for the perturbed equation 1.1) and we introduce its limiting problem, and in Section 3 we obtain convergence in H 1 -norms using the corrector approach. In the appendix we show smoothness of the auxiliary solutions with respect to the perturbation on the representative cell.. Preliminaries In the following, will denote a small positive parameter which will converges to. In order to setup the thin domain 1.) we consider a function b : [, 1] [, ) and a family of positive functions G : [, 1], ), satisfying the following hypothesis H1) The function b belongs to C 1, 1), it is continuous in [, 1] and satisfies bx) b 1 for some positive constant b 1 and all x [, 1]; H) There exist two positive constants G, G 1 such that < G G x) G 1, x [, 1],, )..1) Furthermore, the functions G ) have the following structure G x) = Gx, x/), where the function G : [, 1] R, ) x, y) Gx, y).) is of class C 1 in, 1) R, continuous in [, 1] R, L-periodic in the second variable y, i.e., Gx, y + L) = Gx, y) for all x, y) [, 1] R. We stress for the fact that R varies in accordance with the positive parameter and, when goes to, the domains R collapse themselves to the unit interval of the real line, therefore, in order to preserve the relative capacity of a mensurable subset O R we rescale the Lebesgue measure by a factor 1/, i.e., we deal with the singular measure ρ O) = 1 O, widely considered in thin domains problems see for instance [, 9, 3]). With this measure we introduce the Lebesgue and Sobolev spaces L R ; ρ ) and H 1 R ; ρ ). The norms in these spaces will be denoted,

4 4 M. C. PEREIRA AND R. P. SILVA respectivelly, by L R ) and H1 R ), which are induced by the inner products u, v) = 1 R u v dx 1 dx, u, v L R ) and a u, v) = 1 R { u v + u v}dx 1 dx, u, v H 1 R ). Remark.1. The - norms and the usual ones in L R ) and H 1 R ) are equivalents and easily related u L R ) = 1/ u L R ), u L R ) u H 1 R ) = 1/ u H 1 R ), u H 1 R ). The variational formulation of 1.1) is: find w H 1 R ) such that { } w ϕ + w ϕ dx 1 dx = f ϕ dx 1 dx, ϕ H 1 R ),.3) R R which is equivalent to find w H 1 R ; ρ ) such that a ϕ, w ) = ϕ, f ), ϕ H 1 R ; ρ )..4) We also observe that the solutions w satisfy a priori estimates uniformly in. In fact, taking ϕ = w as a test function we derive from.3) and.4) that w L R ) + w L R ) f L R ) w L R ) w L R ) + w L R ) f L R ) w L R )..5) Now, in order to capture the limiting behavior of a w, w ), as, we consider the sesquilinear form a in H 1, 1) defined by a u, v) = 1 {r u x v x + p u v} dx, u, v H 1, 1),.6) where r :, 1) R and p :, 1) R, defined in.1), are positive and bounded functions. We also consider the following inner product, ) in L, 1) given by u, v) = 1 p u v dx, u, v L, 1)..7) Remark.. There exists a relationship among the function p, the period L of the function Gx, ) and the measure of the representative cell, defined in.1)), by the expression L px) =, x, 1)..8) With such considerations, one can show that the homogenized problem associated to 1.1) is the Neumann problem 1 { } 1 r w x ϕ x + p w ϕ dx = ˆf ϕ dx, ϕ H 1, 1),.9) where r, p :, 1) R + are positive functions defined by rx) = 1 { } 1 y Xx)y, z) dydz, L px) =, L.1)

5 CORRECTORS FOR THE NEUMANN PROBLEM IN THIN DOMAINS 5 Xx) is the unique solution of the auxiliary problem Xx) = in in the representative cell given by Xx) N y, gy)) = N 1x) on B 1 x) Xx) N y, z) = on B x) Xx), z) L - periodic Xx) dydz =.11) = {y, z) R : < y < L, bx) < z < Gx, y)},.1) where B 1 x) is the upper boundary, B x) is the lower boundary of for x, 1) and Nx) = N 1 x), N x)) is the unitary outside vector field to. The function ˆf L, 1) is such that the sequence defined by Gξ) ˆf ξ) = 1 f ξ, s) ds, ξ, 1), bξ) satisfies ˆf ˆf, w L, 1). Indeed, it follows from [6, Theorem.3] that if w H 1, 1) is the unique solution of.9), then the solutions w of 1.1) satisfies w w L R ), as..13) In the present paper we introduce a suitable corrector function κ H 1 R ), with κ = o) in L R ) insomuch w w κ H 1 R ), as. Remark.3. i) If the function rx) is continuous, the integral formulation.9) is the weak form of the homogenized equation 1 px) rx)w x) x + w = fx), x, 1).14) w x ) = w x 1) = with fx) = ˆfx)/px). ii) If initially the function f x, y) = f x), then it is not difficult to see that ˆf x) = G x) + bx))f x) and G x) + bx) L = px), w L, 1), as, and therefore, ˆfx) = px)f x)..1. Smoothness of X. The auxiliary functions Xx) introduced by the problem.11) are originally defined in the representative cell, that depends on x, 1). Such variable performs a perturbation on changing its height at the direction z which also has influence on the auxiliary solution Xx). Since such a domain perturbation is smooth so is Xx) H 1 Y ) with respect to x. In fact we are able to compute at least one derivative with respect to x, 1) by using Corollary A.3 proved in Appendix A. There we combine methods from [6, Proposition A.] and [3, Chapter ] where more general results can be found. This is the reason of our smoothness assumptions to the functions b and G in H1) and H). Thus, in order to consider Xx) as a function in the thin domain R, with the purpose to define the corrector term κ, we first use the periodicity at variable y to extend them to the band Y x) = {y, z) R : y R, bx) < z < Gx, y)}, for all x, 1). Then, we properly compose them with the diffeomorphisms defining, at the end, the function T : R R : x 1, x ) x 1 /, x /), for, ),.15) Xx) T : R R for all x, 1)..16)

6 6 M. C. PEREIRA AND R. P. SILVA Finally, with some notational abuse, we consider the function X :, 1) R R, Xy, y, z) = Xy) T y/, z/).17) correlating variables x and y of the original function X given by.11). In order to simplify our notation at the analysis below we still denote functions.17) by X = Xx 1 ) x 1 /, x /) whenever x 1, x ) R..18) Now, under these considerations, we can obtain some estimates for X on R. Indeed, X L R ) = R Xx 1 )x 1 /, x /) dx 1 dx = 1/ 1/L C Ĉ k=1 Gy1,y ) by 1) Y ) Xy 1 )y 1, y ) dy 1 dy X)y, z) dydz for all, ˆ), for ˆ > small enough with Ĉ > independent of ). This is provided by the continuity of the application [, 1] s Xs) T H 1 R ) at s =, for each >. Thus, using the notation introduced at.18), we have that X L R ) = O )..19) Remark.4. Analogously, we can introduce similar notations for the derivatives of X with respect to variables x, y and z, obtaining that the norms x X L R ), y X L R ), and z X L R ) satisfy.19) for all, ˆ). 3. First-order Corrector As we mentioned before the solutions w of 1.1) doesn t converge in general in H 1 -norms. Therefore, following Bensoussan, Lions and Papanicolaou in [9], we use the auxiliary solution X to introduce the first-order corrector term κ x1 x 1, x ) = Xx 1 ), x ) dw x 1 ), x 1, x ) R. 3.1) dx 1 Remark 3.1. Let us still observe that the function κ defined in 3.1) is not the standard corrector according to [9, 17, 8] since its dependence with respect to variable x 1 involves the changeable auxiliary solution X with respect to the representative cell whose smoothness has been discussed in Section.1. As a matter of fact, we believe that this is the main contribution of our work established here. Now, we have elements to proof the following result Theorem 3.. Let w be the solution of problem 1.1) with f L R ) satisfying f L R ) C for some C > independent of. Consider the family of functions ˆf L, 1) defined by Gx 1) ˆf x 1 ) = 1 f x 1, x ) dx 3.) bx 1) with G ) and b ) satisfying hypotheses H1) and H). If ˆf ˆf w-l, 1), then lim w w κ H1 R ) =, 3.3)

7 CORRECTORS FOR THE NEUMANN PROBLEM IN THIN DOMAINS 7 where κ is the first-order corrector of w defined in 3.1), and w H, 1) C 1, 1) is the unique solution of the homogenized equation.14) with fx) = ˆfx)/px). 3.4) Proof. We have a ϕ, w ) = ϕ, f ), for all ϕ H 1 R ). Taking w + κ H 1 R ) 1 as a test function, we obtain that w w κ H 1 R ) = a w w κ, w w κ ) = a w, w w κ ) a w + κ, w ) + a w + κ, w + κ ) = w w + κ ), f ) + a w + κ, w + κ ). Now, performing the change of variables x, y) x, y/) on problem 1.1), it follows from [6, Theorem.3], as observed in.13) that 1 w w L R ) as, i.e., w w L R ). Hence 3.5) w w, f ) w w L R ) f L R ), as. 3.6) From estimate.19), we also obtain uniformly for x, 1)) that κ, f ) 1 κ L R ) f L R ) Xx) L R ) f L R ) dw /dx L,1) 3.7) Ĉ, as. Moreover, from the convergence ˆf ˆf w-l, 1), we obtain that 1 w, f ) = 1 w x 1 ) = 1 1 w x) ˆf x) dx Gx1,x 1/) bx 1) Therefore, we get from 3.4), 3.6), 3.7), 3.8) and.7) that Next, we will show that f x 1, x ) dx dx 1 w x) ˆfx) dx, as. 3.8) w w + κ ), f ) w, f). 3.9) a w + κ, w + κ ) a w, w ), as. 3.1) In order to do that, we need to pass to the limit on the expression a w + κ, w + κ ) = R 1 { w + κ ) + w + κ } dx 1 dx { = 1 dw 1 y X) + z X ) } + w dx 1 dx 3.11) R dx 1 { dw 1 y X) x X + dw d ) } w R dx 1 dx 1 dx X dw + w X dx 1 dx 3.1) 1 dx 1 { )} dw + x X + dw d w dx 1 dx 1 dx X xx + X d w + dw dx 1 dx 3.13) 1 dx 1 dx 1 R 1 Note that here w is considered to be a function of both variables x 1 and x, simply by w x 1, x ) = w x 1 ).

8 8 M. C. PEREIRA AND R. P. SILVA when goes to zero. By estimates in.19), we shall have to consider the parcel { 1 dw x1 1 y Xx 1 ) R dx 1, x )) x1 + z Xx 1 ), x ) ) } + w x 1 ) dx 1 dx 1 Gx,x/) { dw x )) x ) ) } = 1 y Xx) bx) dx, z + z Xx), z + w x) dzdx 1 dw Gx,x/) { x )) x ) } = 1 y Xx) dx, z + z Xx), z dzdx 3.14) + bx) 1 w x) Gx, x/) + bx)) dx. 3.15) Since the function Gx, y) is L-periodic at variable y, we also have that Gx,y) { Φy) = 1 y Xx) y, z)) + z Xx) y, z) } dz bx) is L-periodic for each x, 1) fixed. Consequently, we obtain from [17, Theorem.6] that as. Therefore 1 dw R dx 1 1 Φx/) 1/L L { x1 1 y Xx 1 ) dw dx Φs) ds, w L, 1),, x as. See also [7, Lemma 4.]. Now, since Xx) satisfies.11), we have that y,z Xx) dydz = )) x1 + z Xx 1 ), x ) } dx 1 dx 1 { 1 y Xx) y, z)) + z Xx) y, z) } dydzdx 3.16) L B 1 N 1 Xx) ds where N = N 1, N ) is the unit outward normal to B 1, the upper boundary of. Hence, we obtain from ) Xx) N 1 Xx) ds = div y,z dydz = B 1 y Xx) dydz that, for all x, 1), y,z Xx) dydz = y Xx) dydz. 3.17) Thus, due to 3.17), we have = = { 1 y Xx) y, z)) + z Xx) y, z) } dydz { 1 y Xx) y, z) + y,z Xx) } dydz {1 y Xx) y, z)} dydz = rx) L, 3.18) where the positive function rx), called homogenization coefficient, is introduced at.1). This theorem is often called Average Convergence Theorem for Periodic Functions. See also [18].

9 CORRECTORS FOR THE NEUMANN PROBLEM IN THIN DOMAINS 9 Then, putting 3.16) and 3.18) together, we get the following limit 1 dw { x1 1 y Xx 1 ) R dx 1, x )) x1 + z Xx 1 ), x ) } dx 1 dx 1 rx) dw dx, as. 3.19) dx Back to the parcel 3.15), we notice that as in [17, Theorem.6] and [7, Lemma 4.] we have that 1 w x) Gx, x/) + bx)) dx 1 px) w x) dx, as, 3.) where the function p is defined at.1), and it is related to the measure of the representative cell. Therefore, the limit 3.1) follows from 3.19) and 3.). Finally, according to 3.5), we complete the proof getting w w κ H 1 R ) from 3.1), 3.9), and using that w satisfies.14). Appendix A. a w, w ) w, f) =, A perturbation result in the basic cell. In the proof of the main results in Section 3 we have used the smoothness of the auxiliary solution X given by.11) with respect to variable x, 1) evaluating the function x X. In order to get such a derivative we analyze here, without loss of generality, how the function X H 1 ), solution of Xx) = in in the representative cell defined by Xx) N y, gy)) = N 1x) on B 1 x) Xx) N y, z) = on B x) Xx), z) L - periodic Xx) dydz = A.1) = {y, z) R : < y < L, < z < Gx, y)}, A.) depends on the function G where B 1 x) and B x) are respectively the upper and lower boundary of for x, 1), and Nx) = N 1 x), N x)) is the unitary outside vector field to. For this we consider the following class of admissible functions AM) = {G C 1 R), L periodic, < G G ) G 1, G s) M} endowed with the norm G = max i 1 sup ξ R ξ igξ). We call by Y G) the basic cell A.) and by XG) the unique solution of problem A.1) for G varying in AM). Following the approach applied in [,, 3] we begin by making the transformation on the cell Y G) into the domain Y Ĝ) for each Ĝ AM) LĜ : Y G) Y Ĝ) z 1, z ) z 1, ˆF z 1 ) z ) = y 1, y ) where ˆF z) = Ĝz) Gz). The Jacobian matrix for this transformation is JLĜz 1, z ) = 1 ˆF z 1 )z ˆF z1 ) ) A.3)

10 1 M. C. PEREIRA AND R. P. SILVA whose determinant is given by JLĜz 1, z ) = ˆF z 1 ) = Ĝz 1) Gz 1 ) that is a rational function with respect to G and Ĝ satisfying A.4) < G G 1 ˆF z) G 1 G, for all z R. A.5) Consequently, using LĜ we have that problem A.1) in the cell Y Ĝ) is equivalent to 1ˆF divbĝw ) = in Y G) BĜW N = on B G) BĜW N = Ĝ on B 1+ Ĝ) 1 G) where BĜW )z 1, z ) = W, z ) = W, z ) Y G) W JL Ĝ dz 1dz = L-periodic ˆF z 1 ) W z 1 z 1, z ) ˆF W z 1 ) z z z 1, z ) ˆF W z 1 ) z z 1 z 1, z ) + ˆF z ˆF z z 1 )) ) W 1) z z 1, z ) ) A.6). A.7) More precisely, we have that XĜ) is the solution of A.1) in Y Ĝ) if and only if W Ĝ) = XĜ) L Ĝ satisfies equation A.6) in Y G). Further, if we define UĜ) = W Ĝ) 1 Y W G) Ĝ) Y G) then, UĜ) is the unique solution of the problem 1ˆF divbĝu) = in Y G) BĜU N = on B G) BĜU N = Ĝ on B 1+ Ĝ) 1 G) U, z ) = U, z ) Y G) U dz 1dz = L-periodic. A.8) We also have that, for each Ĝ AM), the variational formulation of A.8) is given by the bilinear form ρĝ : HG) HG) U, V ) R Y G) ) BĜU dz 1 dz VˆF where HG) is the Hilbert space defined by { HG) = U H 1 Y G)) : U, z ) = U, z ) is L-periodic, Y G) U = and endowed with the H 1 Y G)) norm. In fact, ρĝ is a coercive bilinear form since we have imposed the condition U = on the space HG). Y G) Therefore, XĜ) is the solution of A.1) in Y Ĝ) if and only if satisfies UĜ) = XĜ) L Ĝ 1 Y G) ρĝu, V ) = B 1G) Ĝ V 1 + Ĝ ) ˆF Y G) XĜ) L Ĝ } A.9) A.1) ds, for all V HG). A.11)

11 CORRECTORS FOR THE NEUMANN PROBLEM IN THIN DOMAINS 11 In particular the coercive bilinear form associated to A.1) in Y G) is given by in which the weak formulation is ρ G U, V ) = ρ G : HG) HG) B 1G) U, V ) R Y G) G 1 + G ) V ds, U V dz 1 dz for all V HG). Thus we have performed the boundary perturbation problem given by A.1), where the function spaces are varying with respect to Ĝ AM), into a problem of perturbation of linear operators defined by the bilinear form A.9) and condition A.11) in the fixed space HG). Remark A.1. In order to obtain the smoothness of auxiliary solutions X with respect to Ĝ, we first observe that the maps ρ and γ defined by and ρ : AM) HG) HG) R Ĝ, U, V ) ρ Ĝ U, V ) A.1) γ : AM) HG) Ĝ, V ) R B 1G) Ĝ 1+Ĝ ) VˆF ds A.13) are Fréchet differentiable in the Banach spaces AM) HG) HG) and AM) HG). Indeed, setting M > and G AM), we have for each Ĝ HG) that ρĝ,, ) is a coercive bilinear form in HG) HG), and γĝ, ) is a linear functional on HG). Moreover, due to expression A.4), A.5) and A.7), we have that ρ and γ are rational functions with respect to Ĝ AM) and its derivatives. Therefore, it follows from [1] that ρ and γ defined by A.1) and A.13) are analytic maps in AM) HG) HG) and AM) HG) respectively see [3, Cchapter ]). Proposition A.. Let us consider the family of solutions UĜ) of problem A.8) for each Ĝ AM) where AM) is the functional space defined in A.3) for some constant M >. Then, the map Ĝ UĜ) HG) is analytic in AM). Proof. Using the applications ρĝ,, ) is a coercive bilinear form, and γĝ, ) is a linear functional, we can define the linear operator ZĜ) : HG) HG), and the linear functional SĜ) : HG) R given by < ZĜ)φ, ϕ >= ρĝ, φ, ϕ) and SĜ)ϕ = γĝ, ϕ) such that ZĜ)UĜ) + SĜ) =, for all Ĝ AM) where HG) is the dual space of HG) and UĜ) is the unique solution of problem A.8) for each Ĝ AM) that satisfies expression A.11). Hence, due to Remark A.1, we have that the desired result is a direct consequence of the Implicit Functional Theorem applied to the map Φ : AM) HG) HG) Ĝ, U) ZĜ) U + SĜ). Consequently we have the following result: Corollary A.3. Let T : R R be the diffeomorphism defined in.15) for each >, and consider the family of solutions XĜ) given by the problem A.1) for Ĝ AM) where AM) is the functional space defined in A.3) for some constant M >. Then, the map Ĝ XĜ) T H 1 R ) is analytic in AM).

12 1 M. C. PEREIRA AND R. P. SILVA References [1] N. Ansini, A. Braides; Homogenization of oscillating boundaries and applications to thin films; J. Anal. Math. 83 1) [] J. M. Arrieta; Spectral properties of Schrödinger operators under perturbations of the domain; Ph.D. Thesis, Georgia Inst. of Tech. 1991). [3] J.M. Arrieta, S.M. Bruschi; Very rapidly varying boundaries in equations with nonlinear boundary conditions. The case of a non uniformly Lipschitz deformation; Discrete and Continuoud Dinamycal Systems B 14 ) 1) [4] J. M. Arrieta, A. N. Carvalho, M. C. Pereira and R. P. Silva; Semilinear parabolic problems in thin domains with a highly oscillatory boundary; Nonlinear Analysis: Theory Methods and Appl ) [5] J. M. Arrieta and M. C. Pereira; Elliptic problems in thin domains with highly oscillating boundaries; Bol. Soc. Esp. Mat. Apl. 51 1) [6] J. M. Arrieta and M. C. Pereira; Homogenization in a thin domain with an oscillatory boundary; J. Math. Pures et Appl ) [7] J. M. Arrieta and M. C. Pereira; The Neumann problem in thin domains with very highly oscillatory boundaries; J. Math. Anal. Appl ) [8] M. Baía and E. Zappale; A note on the 3D-D dimensional reduction of a micromagnetic thin film with nonhomogeneous profile; Appl. Anal. 86 5) 7) [9] A. Bensoussan, J. L. Lions and G. Papanicolaou; Asymptotic Analysis for Periodic Structures; North-Holland 1978). [1] D. Blanchard, A. Gaudiello, G. Griso; Junction of a periodic family of elastic rods with a thin plate. Part II ; J. Math. Pures et Appl. 88 ) 7) [11] D. Blanchard, A. Gaudiello, J. Mossino; Highly oscillating boundaries and reduction of dimension: the critical case; Anal. Appl. 5 ) 7) [1] M. Boukrouche and I. Ciuperca; Asymptotic behaviour of solutions of lubrication problem in a thin domain with a rough boundary and Tresca fluid-solid interface law; Quart. Appl. Math. 64 6) [13] D. Caillerie; Thin elastic and periodic plates; Math. Meth. Appl. Sci ) [14] J. Casado-Díaz, M. Luna-Laynez, F.J. Suárez-Grau; Asymptotic behavior of the Navier-Stokes system in a thin domain with Navier condition on a slightly rough boundary; SIAM J. Math. Anal. 45 3) 13) [15] G. A. Chechkin, A. Friedman, A. L. Piatnitski; The Boundary-value Problem in Domains with Very Rapidly Oscillating Boundary; J. Math. Analysis and Appl ) [16] G. A. Chechkin, A. L. Piatnitski; Homogenization of boundary-value problem in a locally periodic perforated domain, Applicable Analysis, v. 71, 1-4, 1998) [17] D. Cioranescu and P. Donato; An Introduction to Homogenization; Oxford lecture series in mathematics and its applications 1999). [18] D. Cioranescu and J. Saint J. Paulin; Homogenization of Reticulated Structures; Springer Verlag 198). [19] G. Dal Maso, An introduction to Γ-convergence; Birkhǎuser, Boston 1993). [] A. Damlamian, K. Pettersson, Homogenization of oscillating boundaries; Discrete and Continuous Dynamical Systems 3 9) [1] J. Dieudonné; Foundations of Modern Analysis; Academic Press, New York and London 1969). [] J. K. Hale and G. Raugel; Reaction-diffusion equation on thin domains; J. Math. Pures et Appl. 9) 71 1) 199) [3] D. B. Henry; Perturbation of the Boundary in Boundary Value Problems of PDEs; Cambridge University Press 5). [4] J. -L. Lions; Asymptotic expansions in perforated media with a periodic structure; Rocky Mountain J. Math. 1 1) 1998) [5] A. L. Madureira and F. Valentin; Asymptotics of the Poisson Problem in domains with curved rough boundaries; SIAM Journal on Mathematical Analysis 38 7) [6] T. A. Mel nyk and A. V. Popov; Asymptotic analysis of boundary value and spectral problems in thin perforated domains with rapidly changing thickness and different limiting dimensions; Mat. Sb. 3 8) 1) [7] I. Pazanin, F. J. Suárez-Grau; Effects of rough boundary on the heat transfer in a thin-film flow; Comptes Rendus Mécanique 341 8) 1) [8] M. C. Pereira and R. P. Silva, Error estimates for a Neumann problem in highly oscillating thin domains; Discrete and Continuous Dyn. Systems 33 ) 13) [9] M. Prizzi and K. P. Rybakowski; The effect of domain squeezing upon the dynamics of reaction-diffusion equations; Journal of Diff. Equations, 173 ) 1) [3] G. Raugel; Dynamics of partial differential equations on thin domains; Lecture Notes in Math. 169, Springer Verlag 1995). [31] E. Sánchez-Palencia; Non-Homogeneous Media and Vibration Theory; Lecture Notes in Phys. 17, Springer Verlag 198). [3] L. Tartar; The General Theory of Homogenization. A personalized introduction; Lecture Notes of the Un. Mat. Ital, 7, Springer-Verlag, Berlin 9).

13 CORRECTORS FOR THE NEUMANN PROBLEM IN THIN DOMAINS 13 Marcone C. Pereira) Escola de Artes, Ciências e Humanidades, Universidade de São Paulo, São Paulo SP, Brazil address: marcone@usp.br Ricardo P. Silva) Instituto de Geociências e Ciências Exatas, Universidade Estadual Paulista, Rio Claro SP, Brazil address: rpsilva@rc.unesp.br