Waveguides with fast oscillating boundary Текст научной статьи по специальности «Математика»

Waveguides with fast oscillating boundary

G. Cardone

University of Sannio, Department of Engineering, Corso Garibaldi, 107, 82100 Benevento, Italy gcardone@unisannio.it

DOI 10.17586/2220-8054-2017-8-2-160-165

We consider an elliptic operator in a planar waveguide with a fast oscillating boundary where we impose Dirichlet, Neumann or Robin boundary

conditions assuming that both the period and the amplitude of the oscillations are small. We describe the homogenized operator, establish the

norm resolvent convergence of the perturbed resolvent to the homogenized one, and prove the estimates for the rate of convergence. It is shown

that under the homogenization, the type of the boundary condition can change.

Keywords: elliptic operator, unbounded domain, norm resolvent convergence.

Received: 4 January 2017

Revised: 29 January 2017

1. Introduction

We study the problem of homogenization of boundary value problems in domains with a fast oscillating boundary when such boundary is given by the graph of the function x2 = n(e)b(x1e-1), where e is a small positive parameter, n(e) is a positive function tending to zero as e ^ +0, and b is a smooth periodic function. The parameter e describes the period of the boundary oscillations while n(e) is their amplitude.

In previous results, the weak or strong resolvent convergence of the solutions was proved and the resolvents were also treated in various possible norms. In some cases, the estimates for the convergence rate were proven. It was also shown that when constructing the next terms of the asymptotics for the perturbed solutions, one can get estimates for the convergence rate or improve it [1-8]. In some cases, complete asymptotic expansions were constructed [9-12].

One more type of established results is the uniform resolvent convergence for the problems. Such convergence was established just for few models, see [13, Ch. II, Sec. 4], [8]. The estimates for the rates of convergence were also established. In both papers, the amplitude and the period of oscillations were of the same order. At the same time, the uniform resolvent convergence for the models considered in the homogenization theory provided quite strong results. Moreover, recently a series of papers by M. Sh. Birman, T. A. Suslina, V. V. Zhikov and S.E. Pastukhova have stimulated interest in this aspect (see [14-27] and references therein and further papers by these authors). The uniform resolvent convergence was shown to hold true for the elliptic operators with fast oscillating coefficients and the estimates for the rates of convergence were obtained. There are also similar results for some problems in bounded domains, see [26]. Similar results but for the boundary homogenization were established in [28-32]. Here, the Laplacian in a planar straight infinite strip with frequently alternating boundary conditions was considered. Such boundary conditions were imposed by partitioning the boundary into small segments where Dirichlet and Robin conditions were imposed in turn. The homogenized problem involves one of the classical boundary conditions instead of the alternating ones. For all possible homogenized problems, the uniform resolvent and the estimates for the rates of convergence were proven and the asymptotics for the spectra were constructed.

In the present paper, we also consider the boundary homogenization for the elliptic operators in unbounded domains but the perturbation is a fast oscillating boundary. As the domain, we choose a planar straight infinite strip with a periodic fast oscillating boundary; the operator is a general self-adjoint second order elliptic operator. The operator is regarded as an unbounded one in an appropriate L2 space. On the oscillating boundary, we impose Dirichlet, Neumann, or Robin conditions. Apart from a mathematical interest in this problem, as a physical motivation, we can mention a model of a planar quantum or acoustic waveguide with a fast oscillating boundary.

Our main result is the form of the homogenized operator and the uniform resolvent convergence of the perturbed operator to the homogenized one. This convergence is established in the sense of the norm of the operator acting from L2 into W2. The estimates for the rate of convergence are provided. Most of the estimates are sharp. In the case of the Dirichlet or Neumann conditions on the oscillating boundary, the homogenized problem involves the same condition on the mollified boundary no matter how the period and amplitude of the oscillations

behave. Provided the amplitude is not greater than the period (in order), the Robin conditions on the oscillating boundary leads us to a similar condition but with an additional term in the coefficient. If the amplitude is greater than the period, the homogenization transforms the Robin conditions into those of Dirichlet. The last result is in a good accordance with a similar case, treated in [33]. The difference is that in [33], the strong resolvent convergence was proven provided the coefficient in the Robin conditions was positive, while we succeeded to prove the uniform resolvent convergence provided the coefficient is either positive or non-negative and vanishing on the set of zero measure. All the results stated in this paper are proved in [34].

2. Problem and main results

Let x = (xi,x2) be the Cartesian coordinates in R2, e be a small positive parameter, n = n(e) be a nonnegative function uniformly bounded for sufficiently small e, b = b(t) be a non-negative 1-periodic function belonging to C2(R). We define two domains, cf. Fig. 1:

fio := {x : 0 < X2 < d}, := {x : n(e)b(xie-1) < X2 < d}, where d > 0 is a constant, and its boundaries are indicated as:

r := {x : X2 = d}, To := {x : X2 = 0},

{x : X2 = n(e)b(xie 1)}.

By Aij = Aij (x), Aj = Aj(x), A0 = A0(x), i, j = 1, 2, we denote the functions defined on Q0 and satisfying the belongings Aij g W^ (Q0), Aj g W^ (Q0), A0 g Lto(Q0). Functions Aj j, Aj are assumed to be complex-valued, while A0 is real-valued. In addition, functions Aj satisfy the ellipticity condition:

A • • = A ••

ij Ji'

Y, AijZiZj > co(|zi|2 + |Z212), x G fio, Zj G C. (2.1)

i,j=1

By a = a(x) we denote a real function defined on {x : 0 < x2 < J} for some small fixed J, and it is supposed that a G Wl({x : 0 < x2 < d}).

Fig. 1. Domain with oscillating boundary

The main object of our study is the operator:

2 3d A . d d —

-E ^ijdx + E aj- + Ao in (2.2)

dxj ij dxi j dxj dxj j

i,J=1 j=1

subject to Dirichlet conditions on r. On the other boundary, we choose either Dirichlet conditions:

u = 0 on re,

or Robin conditions:

d \ d 2 d 2 _

dV^ + aJu = 0 on r" dV^ = Jdx AJvJ

i,J=1 J=1

where vf = (vf, vf) is the outward normal to re. In the case of Dirichlet conditions on re we denote this operator as , while for Robin conditions it is . The former also includes the case of Neumann conditions since the function a can be identically zero.

Rigorously, we introduce as the lower-semibounded self-adjoint operator in L2(ne) associated with the closed symmetric lower-semibounded sesquilinear form:

2 ( du dv \ 2 / du \ 2 / dv \

Ah^n (u,v):= ^ (Aij tA-XM + ^ (Aj tAh + ^ (u,Aj -XM + (Aou, v)L2(fie), i,j=1 ^ dXj dXi/L2(^E) j=1 V dXj /¿2 (ne) j=1 V dXj/L2(fiE)

in L2(Q£) with the domain D(hDn) := W2,0(Q£, dQe). Hereinafter D(^) is the domain of a form or an operator, and Wj 0(Q, S) denotes the Sobolev space consisting of the functions in (Q) with zero trace on a curve S lying in a domain Q c R2. The operator HRn is introduced in the same way via the sesquilinear form:

C (uv):= ¿(A, £ )l + ± (a, £ „

2 / dv \ + (M,Aj d- ) + (A0u,v)L2(Oe) + ^^(r^

, = 1 V dXj/L2(^)

with the domain S(hRn) := W210(Q£, r).

The main aim of the paper is to study the asymptotic behavior of the resolvents of and HRn as e ^ +0. To formulate the main results we first introduce some additional operators.

By HD we denote operator (2.2) in L2(Q0) subject to Dirichlet conditions. We introduce it by analogy with as associated with the form:

^ (A.,|U-)t (0) + ^ (A,,v (0>

v J ¿/£2(^0) , = 1 V J /¿2(^0)

^ (2.3)

0 — 1 V i/ L5(fio)

_ ^i

j = 1 ^ ^ ¿2(^0)

in L2(O0) with the domain D(hD) := W21,0(Q0, The domain of operator HD is W22,0(Q0, d^0) that can be

shown by analogy with [35, Ch. II, Sec. 7,8], [36, Lm. 2.2].

Our first main result (proved in section 2 in [34]) describes the uniform resolvent convergence for .

Theorem 2.1. Let f G L2(Q0). For sufficiently small e, the estimate:

ll(HeDn - i)-1f - (HD - i)-1f ywl(nE) < CV/2||f yL2(no),

holds true, where C is a constant independent of e and f.

The next four theorems describe the resolvent convergence for operator HRn. Given a0 g W^(r0), let HR

be the self-adjoint operator in L2(00) associated with the lower-semibounded sesquilinear symmetric form:

22

du dv \ / , du

hR(u,-):=£ Uj —+ £ Ai —,

¿,:=1 ^ : :=1 V

2 / dv \

+ X) (U'Ai dT" ) + (A0u,v)L2(^o) + (a0u7 v)L2(ro)^

:=1 V aX^L2(fio)

with the domain D(hR) := W210(Q0, r). It can be shown by analogy with [35, Ch. H, Sec. 7,8], [36, Lm. 2.2] that the domain of HR consists of the functions u G W220(^0, r) satisfying Robin conditions:

f Q \ d 2 d _

' -+ «0 u = 0 on r0, ^-0 := -V) Ai2 ^ - A2. (2.4)

v ' ¿=1

First, we consider the particular case of Neumann conditions on T£, i.e., a = 0. Operator HRn and associated quadratic form are re-denoted in this case by and hNn. By HN, we denote the self-adjoint lower-semibounded operator in L2(Q0) associated with the sesquilinear form hN which is taken for a0 = 0. Its domain is the set of the functions in W220(Q0, r) satisfying boundary conditions (2.4) with a0 = 0. The resolvent convergence in this case is given in following theorem (for the proof see section 3 in [34]).

Theorem 2.2. Let f G L2(Q£). Then for sufficiently small e the estimate

ll(HeNn - i)-1f - (HN - i)-1 f Uw^) < Cn1/2||f yL2(no holds true, where C is a constant independent of e and f.

Assume now a ^ 0. Here we consider separately two cases:

e-1n(e) ^ a = const > 0, e ^ +0, (2.5)

e-1n(e) ^ e ^ +0. (2.6)

The first assumption means that the amplitude of the oscillation of curve r£ is of the same order (or smaller) as the period. The other assumption corresponds to the case when the amplitude is much greater than the period. In what follows, the first case is referred to as a relatively slow oscillating boundary r£ while the other describes relatively high oscillating boundary r£.

We begin with the slowly oscillating boundary. We denote:

1 _

a0(x1) := a(x1, 0) J \J 1 + a2(b'(t))2 dt. (2.7)

o

The proof of the following theorem is given in section 3 in [34].

Theorem 2.3. Suppose (2.5) and let f G L2(fi£). Then, for sufficiently small e, the estimate

ll(H££, - i)-1f - (HR - i)-1f !w21(oe) < C(n1/2(e) + |e-2n2(e) - a2|)||f H^o)

holds true, where function a0 in (2.4) is defined in (2.7), and C is a constant independent of e and f.

We proceed to the case of the highly oscillating boundary r£. Here, the homogenized operator happens to be quite sensitive to the sign of a and zero level set of this function. In the paper, we describe the resolvent convergence as a is non-negative. We first suppose that a is bounded from below by a positive constant. Surprisingly, but here the homogenized operator has the Dirichlet condition on r0 as in Theorem 2.1. The proof of the following Theorem is given in section 4 in [34].

Theorem 2.4. Suppose (2.6),

a(x) ^ c1 > 0, c1 = const, (2.8)

and that the function b is not identically constant. Let f G L2(fi0). Then, for sufficiently small e, the estimate:

ll(HRn - i)-1f - (HD - i)-1f ||wi(nE) < C(n1/2 + e1/2n-1/2) |f ||L2(no) (2.9)

holds true, where C is a constant independent of e and f.

zeroes. An essential assumption is that zero level set of a is of zero measure. We let b* := max b.

[0,1]

In the next theorem, that is proved in section 4 in [34], we still suppose that a is non-negative but can have es. An essential assumpti

Theorem 2.5. Suppose (2.6),

a > 0, (2.10)

and that the function b is not identically constant. Assume also that for all sufficiently small S, the set {x : a(x) < S, 0 < x2 < (b* + 1)n} is contained in an at most countable union of the rectangles {x : |x1 - Xn| < «(S), 0 < x2 < (b* + 1)n}, where «(S) is a some nonnegative function such that «(S) ^ +0 as S ^ +0, and numbers Xn, n G Z, are independent of S, are taken in the ascending order, and satisfy uniform in n and m estimate:

|Xn - Xm| > c> 0, n = m. (2.11)

Let f G L2(fi0). Then, for sufficiently small e, the estimate:

||(H£Rn - i)-1f-(HD - i)-1f Hw^)

< C(n1/2 + e1/2n-1/2S-1/2 + «1/2(S)| In«(S)|1/2)Hf HL2(Oo) ( . )

holds true, where C is a constant independent of e and f, and S = S(e) is any function tending to zero as e ^ +0.

Let us discuss the main results. We first observe that under the hypotheses of all theorems we have the corresponding spectral convergence, namely, the convergence of the spectrum and the associated spectral projectors -see, for instance, [37, Thms. VDL23, VHI.24]. We also stress that in all Theorems 2.1-2.5 the resolvent convergence is established in the sense of the uniform norm of bounded operator acting from L2(fi0) into W21(fi£).

In the case of the Dirichlet conditions on r£, the homogenized operator has the same condition on r0 no matter how the boundary r£ oscillates, slowly or highly. The estimate for the rate of convergence is also universal being O(n1/2). Despite here we consider a periodically oscillating boundary, in the proof of Theorem 2.1 this fact

is not used. This is why its statement is also valid for a periodically oscillating boundary described by the equation x2 = e), where b is an arbitrary function bounded uniformly in e and such that e) G C(R). The estimate

is Theorem 2.1 is sharp, see the discussion in the end of Sec. 2 in [34].

A similar situation occurs if we have Neumann conditions on Ге. Here, Theorem 2.2 says that the homogenized operator is subject to Neumann conditions on Г0 and the rate of the uniform resolvent convergence is the same as in Theorem 2.1, namely, O(e1/2). This estimate is again sharp, as the example in the end of Sec. 3 in [34] shows.

Once we have Robin conditions on Ге, the situation is completely different. If the boundary oscillates slowly, the homogenized operator still has Robin conditions on Г0, but the coefficient depends on the geometry of the original oscillations, cf. (2.7). The estimate for the rate of the resolvent convergence in this case involves an additional term in comparison with the Dirichlet or Neumann cases, cf. Theorem 2.3. The estimate in this theorem is again sharp, see the example in the end of Sec. 3 in [34].

As boundary Ге oscillates relatively highly, the resolvent convergence changes dramatically. If coefficient a is strictly positive, the homogenized operator has Dirichlet conditions on Г0. A new term, e1/2n-1/2, appears in the estimate for the rate of the uniform resolvent convergence, cf. Theorem 2.4. We are able to prove that this term is sharp, see the discussion in the end of Sec. 4 in [34].

Provided function a is non-negative and vanishes only on a set of zero measure, the homogenized operator still has Dirichlet conditions on Г0, but the estimate for the rate of the uniform resolvent convergence becomes worse. Namely, the behavior of a in a vicinity of its zeroes becomes important. This is reflected by functions ^(5) and 5 in (2.12). The latter should be chosen so that 5 ^ +0, e1/2n-1/25-1/2 ^ +0, e ^ +0, that is always possible. The optimal choice of 5 is so that:

M1/2(5)| lnM(5)|1/2 - e1/2n-1/25-1/2,

5^(5) | ln^(5)|~ en-1. (2.13)

As we see, the choice of 5 depends on a particular structure of ^(5). The most typical case is ^(5) — 51/2, i.e., the function a vanishes by the quadratic law in a vicinity of its zeroes. In this case, condition (2.13) becomes:

53/2| ln 5| - en-1,

which implies:

5 - e2/3n-2/3| lnen-1|-2/3. Then, the estimate for the resolvent convergence in Theorem 2.5 is of order O ((n1/2 + e1/6n-1/6| ln en-111/3).

We are not able to prove the sharpness of estimate (2.12), but in the end of Sec. 4 in [34] we provide some arguments showing that estimate (2.12) is rather close to being optimal.

In conclusion, we discuss the case of Robin conditions on highly oscillating Ге when the coefficient a does not satisfy the hypotheses of Theorems 2.4, 2.5. If it is still non-negative but vanishes for a set of non-zero values, and at the end-points of this set the vanishing happens with certain rate like in Theorem 2.5, we conjecture that the homogenized operator involves mixed Dirichlet and Neumann conditions on Г0. Namely, if a(x1,0) = 0 on rN and a(x1,0) > 0 on rD, Г0 = rN U rD, it is natural to expect that the homogenized operator has Neumann conditions on rN and Dirichlet one on Г^. This conjecture can be regarded as the mixture of the statements of Theorems 2.2 and 2.5. The main difficulty of proving this conjecture is that the domain of such homogenized operator is no longer a subset of W22 (П0) because of the mixed boundary conditions. At the same time, this fact was essentially used in all our proofs. An even more complicated situation occurs once a is negative or sign-indefinite. If a is negative on a set of non-zero measure, it can be shown that the bottom of the spectrum of the perturbed operator goes to -то as e ^ +0. In such cases, one should study the resolvent convergence near this bottom, i.e., for the spectral parameter tending to -то. This makes the issue quite troublesome. We stress that under the hypotheses of all Theorems 2.1-2.5, the bottom of the spectrum is lower-semibounded uniformly in e.

References

[1] Friedman A., Hu B., Liu Y. A boundary value problem for the poisson equation with multi-scale oscillating Boundary. J. Diff. Eq., 1997, 137, P. 54-93.

[2] Friedman A., Hu B. A Non-stationary Multi-scale Oscillating Free Boundary for the Laplace and Heat Equations. J.Diff. Eq., 1997, 137, P. 119-165.

[3] Belyaev A.G., Pyatnitskii A.L., Chechkin G.A. Asymptotic behavior of a solution to a boundary value problem in a perforated domain with oscillating boundary. Siberian Math. J., 1998, 39, P. 621-644.

[4] Kozlov V.A., Nazarov S.A. Asymptotics of the spectrum of the Dirichlet problem for the biharmonic operator in a domain with a deeply indented boundary. St. Petersburg Math. J., 2011, 22, P. 941-983.

[5] Chechkin G.A., Chechkina T.P. On homogenization of problems in domains of the "Infusorium" type. J. Math. Sci., New York, 2004, 120, P. 1470-1482.

[6] Mikelic A. Rough boundaries and wall laws. Qualitative properties of solutions to partial differential equations, Lecture notes of Necas Center for mathematical modeling, ed. by E. Feireisl, P. Kaplicky and J. Malek, Vol. 5, Matfyzpress, Publishing House of the Faculty of Mathematics and Physics Charles University in Prague, Prague, 2009, P. 103-134.

[7] Nazarov S.A. The two terms asymptotics of the solutions of spectral problems with singular perturbations. Math. USSR-Sb., 1991, 69, P. 307-340.

[8] Nazarov S.A. Dirichlet problem in an angular domain with rapidly oscillating boundary: Modeling of the problem and asymptotics of the solution. St. Petersburg Math. J., 2008, 19, P. 297-326.

[9] Amirat Y., Chechkin G.A., Gadyl'shin R.R. Asymptotics for eigenelements of Laplacian in domain with oscillating boundary: multiple eigenvalues. Appl. Anal., 2007, 86, P. 873-897.

[10] Amirat Y., Chechkin G.A., Gadylshin R.R. Asymptotics of simple eigenvalues and eigenfunctions for the Laplace operator in a domain with oscillating boundary. Computational Mathematics and Mathematical Physics, 2006, 46, P. 97-110.

[11] Nazarov S.A. Asymptotics of solutions and modeling the problems of elasticity theory in domains with rapidly oscillating boundaries. Izv. Math., 2008, 72.

[12] Gobbert M.K., Ringhofer C.A. An Asymptotic Analysis for a Model of Chemical Vapor Deposition on a Microstructured Surface. SIAM Jour. Appl. Math., 1998, 58, P. 737-752.

[13] Olejnik O.A., Shamaev A.S., Yosifyan G.A. Mathematical problems in elasticity and homogenization. Studies in Mathematics and its Applications. 26, North-Holland, Amsterdam etc., 1992.

[14] Bunoiu R., Cardone G., Suslina T. Spectral approach to homogenization of an elliptic operator periodic in some directions. Math. Meth. Appl. Sci., 2011, 34, P. 1075-1096.

[15] Cardone G., Pastukhova S.E., Zhikov V.V. Some estimates for nonlinear homogenization. Rend. Accad. Naz. Sci. XL Mem. Mat. Appl., 2005, 29, P. 101-110.

[16] Cardone G., Pastukhova S.E., Perugia C. Estimates in homogenization of degenerate elliptic equations by spectral method. Asympt. Anal., 2013, 81, P. 189-209.

[17] Pastukhova S.E., Tikhomirov R.N. Operator Estimates in Reiterated and Locally Periodic Homogenization. Dokl. Math., 2006, 76, P. 548553.

[18] Pastukhova S.E. Some Estimates from Homogenized Elasticity Problems. Dokl. Math., 2006, 73, P. 102-106.

[19] BirmanM.S. On homogenization procedure for periodic operators near the edge of an internal gap. St. Petersburg Math. J., 2004, 15, P. 507-513.

[20] Birman M.S., Suslina T.A. Homogenization of a multidimensional periodic elliptic operator in a neighbourhood of the edge of the internal gap. J. Math. Sciences, 2006, 136, P. 3682-3690.

[21] Birman M.S., Suslina T.A. Homogenization with corrector term for periodic elliptic differential operators. St. Petersburg Math. J., 2006, 17, P. 897-973.

[22] Birman M.S., Suslina T.A. Homogenization with corrector for periodic differential operators. Approximation of solutions in the Sobolev class H 1(Rd). St. Petersburg Math. J., 2007, 18, P. 857-955.

[23] Vasilevskaya E.S., Suslina T.A. Threshold approximations for a factorized selfadjoint operator family with the first and second correctors taken into account. St. Petersburg Math. J., 2012, 23, P. 275-308.

[24] Zhikov V.V. On operator estimates in homogenization theory. Dokl. Math., 2005, 72, P. 534-538.

[25] Zhikov V.V. Spectral method in homogenization theory. Proc. Steklov Inst. Math., 2005, 250, P. 85-94.

[26] Zhikov V.V. Some estimates from homogenization theory. Dokl. Math., 2006, 73, P. 96-99.

[27] Zhikov V.V., Pastukhova S.E., Tikhomirova S.V. On the homogenization of degenerate elliptic equations. Dokl. Math., 2006, 74, P. 716720.

[28] Borisov D., Bunoiu R., Cardone G. On a waveguide with frequently alternating boundary conditions: homogenized Neumann condition. Ann. H. Poincare, 2010, 11, P. 1591-1627.

[29] Borisov D., Bunoiu R., Cardone G. On a waveguide with an infinite number of small windows. C.R. Mathematique, 2011, 349, P. 53-56.

[30] Borisov D., Bunoiu R., Cardone G. Homogenization and asymptotics for a waveguide with an infinite number of closely located small windows. J. of Math. Sci., 2009, 176, P. 774-785.

[31] Borisov D., Bunoiu R., Cardone G. Waveguide with non-periodically alternating Dirichlet and Robin conditions: homogenization and asymptotics. Z. Ang. Math. Phys., 2013, 64 (3), P. 439-472.

[32] Borisov D., Cardone G. Homogenization of the planar waveguide with frequently alternating boundary conditions. J. of Phys. A: Mathematics and General, 2009, 42, P. 365205 (21 pp.).

[33] Chechkin G.A., Friedman A., Piatnitski A.L. The Boundary-value Problem in Domains with Very Rapidly Oscillating Boundary. J. Math. Anal. Appl., 1999, 231, P. 213-234.

[34] Borisov D., Cardone G., Faella L., Perugia C. Uniform resolvent convergence for a strip with fast oscillating boundary. J. Diff Eq., 2013, 255, P. 4378-4402.

[35] Ladyzhenskaya O.A., Uraltseva N.N. Linear and quasilinear elliptic equations. Academic Press, New York, 1968.

[36] Borisov D. Asymptotics for the solutions of elliptic systems with fast oscillating coefficients. St. Petersburg Math. J., 2009, 20, P. 175-191.

[37] Reed M., Simon B. Methods of mathematical physics. Functional analysis, Academic Press, 1980.