Научная статья на тему 'Carleman Formula for the Maxwell Equation in a Cap Type Domain'

Carleman Formula for the Maxwell Equation in a Cap Type Domain Текст научной статьи по специальности «Математика»

CC BY
165
49
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
Carleman formula / scattering / elliptic complex / Green formulas / Stratton-Chu formulas / Cauchy problem. / формула Карлемана / рассеяние / эллиптический комплекс / формулы Грина / фор- мулы Стрэттона-Чу / задача Коши.

Аннотация научной статьи по математике, автор научной работы — Komil O. Makhmudov

We consider a Cauchy problem for the Maxwell equations in a cap type domain X in R3. We find reasonable solvability conditions and a Carleman formula for its solution.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Формула Карлемана для уравнения Максвелла в области типа шапки

Рассмотрим задачу Коши для уравнений Максвелла в области типа шапки X в R3. Мы указываем разумное условие разрешимости и формулу Карлемана для ее решения.

Текст научной работы на тему «Carleman Formula for the Maxwell Equation in a Cap Type Domain»

УДК 517.5

Carleman Formula for the Maxwell Equation in a Cap Type Domain

Komil O. Makhmudov*

Department of Mechanics and Mathematics University of Samarkand University Boulevard, 15, Samarkand, 140104

Uzbekistan

Received 28.01.2019, received in revised form 06.02.2019, accepted 06.03.2019

We consider a Cauchy problem for the Maxwell equations in a cap type domain X in R3. We find reasonable solvability conditions and a Carleman formula for its solution.

Keywords: Carleman formula, scattering, elliptic complex, Green formulas, Stratton-Chu formulas, Cauchy problem.

DOI: 10.17516/1997-1397-2019-12-3-317-322.

Introduction

The moral of the Cauchy problem for elliptic equations is that, after J.Hadamard, a Cauchy problem for a Laplace equation is ill-posed. Of course, the instability refers here to the standard setting, for in spaces with the two-norm convergence the Cauchy problem for elliptic equations proves to be stable, see [1].

The classical equations of electrodynamics are usually formulated in terms of the de Rham complex in R3. In [4], we formulated Maxwell's equations within the framework of arbitrary elliptic complexes on a compact manifold X with boundary. Using the methods of [9] we studied both the Hilbert problem and the Cauchy problem with data on a part of the boundary dX for solutions of the corresponding stationary equations. The character of instability, solvability criteria and regularisation methods of the Cauchy problem for elliptic equations are studied in [9, 5]. For the complete bibliography we refer the reader to these books. Much of the theory developed in [9] extends immediately to other ill-posed problems of complex analysis or partial differential equations. N. Tarkhanov and A. Shlapunov has published a continuation formula for a larger class of boundary value problems not only for elliptic systems but also for elliptic complex.

Examples of bases with double orthogonality can be found in [6]. When working with domains such as a spherical layer, the main technique is to decompose the elements of a suitable space into a series of homogeneous harmonic functions that form a basis on the sphere. The main difficulty in the implementation of the proposed idea in practice is the choice of a basis with the property of double orthogonality in a given domain.

In this paper we will describe a simpler Carleman formula for the Maxwell equations, i.e. we consider the problem of analytic continuation of a solution of Maxwell equations in a spatial bounded domain from data on part of the boundary of the domain.

An explicit Carleman function for this problem was first described in [10]. Based on the ideas and methods developed by Sh.Yarmukhamedov, we will construct an example of the Carleman function of the Cauchy problem for the Maxwell equations In the limit, it and its derivative disappear outside an arbitrary fixed cone. It decreases rapidly enough at infinity.

* komil.84@mail.ru © Siberian Federal University. All rights reserved

1. Maxwell's equations

The Maxwell equations on an n-dimensional compact manifold X with boundary have the form

eE[ = - aE + d*H, ¡iH[ = - dE,

with E and H being unknown functions of t = [0, T] with values in differential forms of degree i and i + 1, respectively, and e, n and a positive constants. Substituting

E (x,t) = (e + ia/u)-1/2e-lut u(x,t), H (x,t) = p-1/2e-lut f (x,t)

we get

(e^u/k) u't = iku + d*f, - (k/u) ft = du - ikf

in X x [0,1], where u is a nonzero real constant and k2 = (e + ia/u)^u2. The sign of k is chosen from the condition Qk ^ 0. Equations (2) still make sense within the framework of arbitrary complexes of differential operators on X, see [4]. The operator on the right-hand side of (2) proves to be elliptic in the sense of Douglis-Nirenberg (see Example below). The Cauchy data on dX of forms u and f of Sobolev class H 1(X) with respect to the stationary Maxwell operator consist of t(u), the tangential part of u, and n(f), the normal part of f. Both t(u) and n(f) are differential forms of degree i on the boundary. Further we discuss the first mixed problem for solutions of Maxwell's equations (2) in the cylinder X x [0, T]. It stems from scattering of incident electromagnetic waves by a perfectly conducting body. In this case the tangential part t(u) of the "electric" field u must vanish on the body surface dX. Hence, we pose the initial conditions on the lower basis of the cylinder and a Dirichlet condition on the lateral surface. Let X be a compact differentiable manifold of dimension 3 with or without boundary. More explicitly, we use the de Rham complex in R3

0 Q0(R3) Q 1(R3) A Q2(R3) A n3(R3) 0,

where Q0(R3) is the space of smooth functions, Q 1(R3) is the space of 1-forms, and so forth. Forms which are the image of other forms under the exterior derivative, plus the constant 0 function in Q0(R3) are called exact and forms whose exterior derivative is 0 are called closed (see closed and exact differential forms); the relationship dod = 0 then says that exact forms are closed. Then we can think of E as a differential form u of degree 1, H as a differential form f of degree 2, thus identifying curl E with du and curl H with d*f. Here, d* stands for the formal adjoint operator of d. Let us denote by A = d*d + dd* is the Laplacian of the de Rham complex. In this way Maxwell's equations in stationary form can be written as

iku + d*f = 0,

(3)

-ik f + du = 0,

which already make sense not only for differential forms u and f of degree 1 and 2 in R3, respectively, but also for differential forms u and f of degree i and i + 1 in R", where -1 ^ i ^ n.

Definition 1. Let -1 ^ i ^ 3. By the Maxwell operator for the de Rham complex at step i is meant

m=(d

A characteristic and important example of a system elliptic in the sense of Douglis-Nirenberg and not elliptic in the sense of Petrovskii is the Maxwell equations (1.2) in the domain X c R3. In the case of the de Rham complex in R3 at step i = 0, Maxwell's equations (1.2) have the form:

iku — divf = 0 —ikf + gradu = 0

or

—iku + divf = 0 —ikf + gradu = 0

(4)

The system (1.3) with constant coefficients and its operator is:

( —ik d

M0 =

d d

3x\ dx2 dx3

— ik 0 0

0 — ik 0

0 0 —ik

dxi d

dx2 d

\dx3

We try to find the vectors s and t with integer components that simultaneously satisfy the conditions of definitions.

For the system (1.3) we choose the weights , for example, the form

s = (2,1,1,1) and t = (0, -1,-1, -1).

The corresponding main part will have the form

M0AD)

0 d d

3x\ 3x2 8x3

d

—ik 0 0

dxi

d

0 —ik 0

dx2

d

\8x3 0 0 —ik

and its determinant

det M0tt(D) = k2(£ + $ + £3)=0, V£ e R3 and k = 0.

Using definition, we come to the conclusion that the system (1.3) is elliptic in the sense of Douglis-Nirenberg.

By definition, M1 is a first order differential operator from sections of F1 © F1+1 to sections of the same bundle over X. This operator fails to be elliptic of order 1 in the classical sense unless N = 2. On the other hand, applying d* to both sides of iku + d*f = 0 we conclude that di-1*u = 0 unless k = 0. Analogously, from -ikf + du = 0 it follows that di+1f = 0 unless k = 0. Complementing Maxwell's equations by their differential consequences di-1*u = 0 and di+1f = 0 yields a system of first order differential equations for u and f, whose classical symbol is injective. Another way of stating this is to say that there is a differential operator G1 from sections of F1 © F1+1 to sections of the same bundle, such that GrM1 is a second order differential operator on X elliptic in the classical case. As is usual in homological algebra, we will omit the index i of Ml i.e. Ml = M when it is clear from the context. An easy computation shows that

ik + (1/ik)dd* d

d*

—ik — (1/ik)d*d

Lemma 1. As defined above, C satisfies

CM = MC = ( A - k2 A - k2 ).

Proof.

CM = MC = ( -k2 + dd* + d*d ikd* + (1/ik)d*d* -

V ikd - ikd - (1/ik)d*dd dd* - k2 + d*d

( A - k2 0 \ I 0 A - k2 )■

Lemma 2. The pseudodifferential operator

= ( G (ik +(1/ik)dd*) Gd* \

* V Gd G (-ik - (1/ik)d*d) )

is a left fundamental solution of the Maxwell operator M on X'. Proof. From Lemma 1 it follows immediately that

*=(G G) o c

is a left fundamental solution of M. It remains to substitute the explicit expression (5) for C. □

2. The Carleman formula

We now turn to the classical Maxwell equation in a three-dimensional space, which can be a three-dimensional manifold X' as well. To demonstrate our constructions along more classical lines, we consider the case X' = R3. As mentioned in Section 1, the classical Maxwell equations have the form

ikE + d*H = 0, -ikH + dE = 0,

E and H being functions in a closed domain X c R3 with values in R3. If E is suitably specified within 1 -forms and H within 2 -forms, both the exterior derivative d and its formal adjoint d* can be identified with the operator curl on vector fields in R3.

Applying Corollary 3.3 [4] to the classical Maxwell equations we obtain the Stratton-Chu formula [7].

Theorem 1. Suppose (E, H) is an electromagnetic wave in X whose electric component E and magnetic component H are both continuous up to the boundary. Then

( (1/ik)d*d -d* ) [ -1 exp(ik|x - y\) ( in(H) \ = (E(x)\ I -d -(1/ik)dd* ) JdX 4n \x - y\ \-v A t(E)) S = \H(x)J

for all x £ X \ dX, and the left-hand side vanishes away from X.

Let a be a positive number. Consider the entire function K(w) = exp(aw2) of complex variable w £ C. The restriction of K to any vertical line w = u + iv just amounts to K(u + iv) = = K(u) exp(2iauv - av2), which is a rapidly decreasing function of v.

Assume that X is a bounded domain in the upper half-space {x3 > 0} of R3 whose boundary consists of a smooth surface S lying in the half-space {x3 > 0}, and a closed piece of the plane

{x3 = 0}. Such domains are usually referred to as cap type domains. Note that the unit outward normal vector on the piece dX \ S just amounts to (0,0, —1). We consider the problem of finding an electric field E and a magnetic field H in X with given tangential component E0 of E and normal component H0 of H on S. Given two different points x = (x',x3) and y = (y',y3) in R3, set r' = \y' — x'l and introduce the integral

,,,-11 ( K(w) N cos kê JQ

2n2 K(x3) Jo \w — x3> A/r'2 + $2 where w = y3 + i\/r'2 + $2. An easy calculation shows that

(x,y) = pa (x,y; $) cos k$d$, (6)

o

where pa (x, y; $) is given by

— 1 ea(y2-x2)e-a(r'2+^2) [, sin2ay3Vr'2 + $2 /—n-777\

--((y3 — x3) " — cos2aWr '2 + *2) ■

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Hence it follows that

1 e-a'xl e-°(r'2+^2)

Pa (x,y; $) =

2n2 ê2 +

r2

1 e-ax2 e-a(r'2+V2)

dV3Pa (x,y; $) = -($2 + r2)2-x3 i1 + a($2 + r2))

on the plane y3 = 0. If a = 0 and K(0) = 1 then kernels (x,y) is a classical fundamental solution of the Helmholtz equation. Substituting (2.5) into the formula for the fundamental solution &(x — y) of M we conclude that (x,y) is a Carleman function of the Cauchy problem in the domain X with data on S, parametrised by a.

T , , ( (1/ik)d*d —d* \ , , , (x>y)={( -d —(1/lk)dd* )** (x,y)■

Lemma 3. The xta(x,y) matrix-valued function on the set R3 x (R3 \{y = x}), satisfying

M(dx)^a (x,y) = 0, M' (dy )($a (x,y))T = 0.

Proof. The proof follows from the properties of the fundamental solution <P(x — y) of the Maxwell operator M. □

Theorem 2. Let (E, H) be an electromagnetic wave in X continuous up to S. Then the limit relation

f E(x)\ f _ , , ( m(

}TooJs K(x,y)(-nf(E))"S

\H (x),

holds uniformly on each compact subset of X .

Proof. As is known, the value of the regular solution of Maxwell equations at a point x inside X in terms of t(E) and n(H) on the boundary X. Using the Stratton-Chu formula of Theorem 2.1

(HSH <->(-ZHE)+!, *-<->( E>)

The convergence of the improper integral on the right-hand side of (x,y) is thus guaranteed by the factor . If a ^ œ the expression (x,y), pa(x,y; ê) and its partial derivative

dy3pa(x,y; ê) tends to zero exponentially on the plane {x3 > 0}. It follows that the part of the boundary integral over dX \ S tends to zero on {x3 = 0}. This establishes the desired formula. □

References

[1] L.A.Aizenberg, N.N.Tarkhanov, Conditionally stable problems and Carleman formulas, Siberian Math. J., 31(1990), no. 6, 9-15.

[2] L.A.Aizenberg, Carleman's Formulas in Complex Analysis. Theory and Applications, Kluwer Academic Publishers, Dordrecht, 1993.

[3] D.Colton, R.Kress, Inverse Acoustic and Electromagnetic Scattering Theory, SpringerVerlag, Heidelberg, 1998.

[4] K.Makhmudov, O.Makhmudov, N.Tarkhanov, Equations of Maxwell type, J. Math. Anal. Appl, 378(2011), no. 1, 64-75.

[5] A.Shlapunov, Green's Integrals and Their Applications to Elliptic systems, Tesi di per-fezionamento, Scuola Normale Superiore, Pisa, 1996.

[6] A.Shlapunov, N.Tarkhanov, Bases with double orthogonality in the Couchy problem for systems with injective symbols, Proc. Londom Math. Soc., 71(1995), no. 1, 1-52.

[7] J.A.Stratton, Electromagnetic Theory, McGraw-Hill, NY, 1941.

[8] N.Tarkhanov, Complexes of Differential Operators, Kluwer Academic Publishers, Dordrecht, NL, 1995.

[9] N.Tarkhanov, The Cauchy Problem for Solutions of Elliptic Equations, Akademie Verlag, Berlin, 1995.

[10] Sh.Yarmukhamedov, The Cauchy problem for the Laplace equation, Soviet Math. Dokl., 18(1977), no. 4, 939-942.

Формула Карлемана для уравнения Максвелла в области типа шапки

Комил О. Махмудов

Механико-математический факультет Самаркандский университет Университетский бульвар, 15, Самарканд, 140104

Узбекистан

Рассмотрим .задачу Коши для уравнений Максвелла в области типа шапки X в К3. Мы указываем разумное условие разрешимости и формулу Карлемана для ее решения.

Ключевые слова: формула Карлемана, рассеяние, эллиптический комплекс, формулы Грина, формулы Стрэттона-Чу, задача Коши.

i Надоели баннеры? Вы всегда можете отключить рекламу.