LECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI

Similar documents
The Factorization Method for Inverse Scattering Problems Part I

The Helmholtz Equation

A coupled BEM and FEM for the interior transmission problem

On the Spectrum of Volume Integral Operators in Acoustic Scattering

Volume and surface integral equations for electromagnetic scattering by a dielectric body

ON THE EXISTENCE OF TRANSMISSION EIGENVALUES. Andreas Kirsch1

[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,

Asymptotic Behavior of Waves in a Nonuniform Medium

MATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Spring 2018 Professor: Jared Speck

ELECTROMAGNETIC SCATTERING FROM PERTURBED SURFACES. Katharine Ott Advisor: Irina Mitrea Department of Mathematics University of Virginia

MATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Fall 2011 Professor: Jared Speck

Laplace s Equation. Chapter Mean Value Formulas

The Factorization Method for a Class of Inverse Elliptic Problems

Partial Differential Equations

u xx + u yy = 0. (5.1)

Lecture 8: Boundary Integral Equations

Wave equation on manifolds and finite speed of propagation

SOLUTION OF POISSON S EQUATION. Contents

Traces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains

THE INTERIOR TRANSMISSION PROBLEM FOR REGIONS WITH CAVITIES

An eigenvalue method using multiple frequency data for inverse scattering problems

arxiv: v1 [math.ap] 21 Dec 2018

SOLUTION OF THE DIRICHLET PROBLEM WITH L p BOUNDARY CONDITION. Dagmar Medková

NEW RESULTS ON TRANSMISSION EIGENVALUES. Fioralba Cakoni. Drossos Gintides

MATH 425, FINAL EXAM SOLUTIONS

Inverse scattering problem with underdetermined data.

ESTIMATES FOR ELLIPTIC HOMOGENIZATION PROBLEMS IN NONSMOOTH DOMAINS. Zhongwei Shen

ON THE MATHEMATICAL BASIS OF THE LINEAR SAMPLING METHOD

Elliptic Problems for Pseudo Differential Equations in a Polyhedral Cone

LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,

A method for creating materials with a desired refraction coefficient

Boundary problems for fractional Laplacians

The oblique derivative problem for general elliptic systems in Lipschitz domains

Fractional order operators on bounded domains

Discreteness of Transmission Eigenvalues via Upper Triangular Compact Operators

Inverse Scattering Theory: Transmission Eigenvalues and Non-destructive Testing

Extended Foldy Lax Approximation on Multiple Scattering

6.3 Fundamental solutions in d

EXISTENCE OF GUIDED MODES ON PERIODIC SLABS

A new method for the solution of scattering problems

Green s Functions and Distributions

Liouville-type theorem for the Lamé system with singular coefficients

Scattering of electromagnetic waves by thin high contrast dielectrics II: asymptotics of the electric field and a method for inversion.

Integral Representation Formula, Boundary Integral Operators and Calderón projection

The Interior Transmission Eigenvalue Problem for Maxwell s Equations

Introduction to the Boundary Element Method

LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)

and finally, any second order divergence form elliptic operator

Fourier Transform & Sobolev Spaces

On spherical-wave scattering by a spherical scatterer and related near-field inverse problems

Weak Convergence Methods for Energy Minimization

RIEMANN-HILBERT PROBLEM FOR THE MOISIL-TEODORESCU SYSTEM IN MULTIPLE CONNECTED DOMAINS

4.2 Green s representation theorem

Chapter 1 Mathematical Foundations

BIHARMONIC WAVE MAPS INTO SPHERES

A Direct Method for reconstructing inclusions from Electrostatic Data

Chapter 7: Bounded Operators in Hilbert Spaces

Homogenization of the Transmission Eigenvalue Problem for a Periodic Media

THE L 2 -HODGE THEORY AND REPRESENTATION ON R n

Proceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005

When is the error in the h BEM for solving the Helmholtz equation bounded independently of k?

2. Function spaces and approximation

The Factorization Method for Maxwell s Equations

Recovery of anisotropic metrics from travel times

Mathematical Foundations for the Boundary- Field Equation Methods in Acoustic and Electromagnetic Scattering

arxiv: v3 [math.ap] 4 Jan 2017

GREEN S IDENTITIES AND GREEN S FUNCTIONS

Lecture 2: A Strange Term Coming From Nowhere

Factorization method in inverse

Elliptic Regularity. Throughout we assume all vector bundles are smooth bundles with metrics over a Riemannian manifold X n.

INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES. L q -solution of the Neumann, Robin and transmission problem for the scalar Oseen equation

Survey of Inverse Problems For Hyperbolic PDEs

Numerical Solutions to Partial Differential Equations

Math The Laplacian. 1 Green s Identities, Fundamental Solution

Preprint 1 1.pdf INTERACTION PROBLEMS OF METALLIC AND PIEZOELECTRIC MATERIALS WITH REGARD TO THERMAL STRESSES

The mathematics of scattering by unbounded, rough, inhomogeneous layers

We denote the space of distributions on Ω by D ( Ω) 2.

NOTES ON SCHAUDER ESTIMATES. r 2 x y 2

New Identities for Weak KAM Theory

Applied Numerical Mathematics. High-order numerical schemes based on difference potentials for 2D elliptic problems with material interfaces

The double layer potential

i=1 α i. Given an m-times continuously

Integral Equations in Electromagnetics

Hamburger Beiträge zur Angewandten Mathematik

The Imaging of Anisotropic Media in Inverse Electromagnetic Scattering

Lecture Notes on PDEs

A global solution curve for a class of free boundary value problems arising in plasma physics

Quantitative Homogenization of Elliptic Operators with Periodic Coefficients

FINITE ENERGY SOLUTIONS OF MIXED 3D DIV-CURL SYSTEMS

Creating materials with a desired refraction coefficient: numerical experiments

Velocity averaging a general framework


A Perron-type theorem on the principal eigenvalue of nonsymmetric elliptic operators

THE DEPENDENCE OF SOLUTION UNIQUENESS CLASSES OF BOUNDARY VALUE PROBLEMS FOR GENERAL PARABOLIC SYSTEMS ON THE GEOMETRY OF AN UNBOUNDED DOMAIN

The second-order 1D wave equation

Determination of thin elastic inclusions from boundary measurements.

Elliptic Kirchhoff equations

Elliptic PDEs of 2nd Order, Gilbarg and Trudinger

On uniqueness in the inverse conductivity problem with local data

THE GREEN FUNCTION. Contents

Transcription:

LECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI Georgian Technical University Tbilisi, GEORGIA 0-0

1. Formulation of the corresponding transmission problem 2. Basic integral relations 3. Equivalent reduction to LBDIE system 4. Uniqueness theorem 5. Investigation of the LBDIO and Existence result 0-1

1. Formulation of the corresponding transmission problem. 2. Basic integral relations. 3. Equivalent reduction to LBDIE system. 4. Uniqueness theorem. 5. Investigation of the LBDIO and Existence result O. Chkadua, S. Mikhailov, and D. Natroshvili, Localized boundary-domain integral equations for acoustic scattering by inhomogeneous anisotropic obstacles. (Prepared for publication) 0-2

O. Chkadua, S. Mikhailov, and D. Natroshvili, Localized boundary-domain singular integral equations based on harmonic parametrix for divergence-form elliptic PDEs with variable matrix coefficients. Integral Equations and Operator Theory, 76, No. 4 (2013), 509 547. O. Chkadua, S. Mikhailov, and D. Natroshvili, Analysis of direct boundary-domain integral equations for a mixed BVP with variable coefficient, Part I: Equivalence and invertibility. J. Integral Equations Appl. 21(2009), No. 4, 499 542. D. Colton and R. Kress, Inverse Acoustic and Electromagnetic Scattering Theory. Third Edition, Springer-Verlag, Berlin, 2013. P. Werner, Zur mathematischen Theorie akustischer Wellenfelder. Arch. Ration. Mech. Anal., 6 (1961), 231 260. P.A. Martin, Acoustic scattering by inhomogeneous obstacles. SIAM Journal of Applied Mathematics, 64, No. 1 (2003), 297 308. 0-3

INTRODUCTION We consider the time-harmonic acoustic wave scattering by a bounded anisotropic inhomogeneity embedded in an unbounded anisotropic homogeneous medium. We assume that the material parameters are functions of position within the inhomogeneous bounded obstacle. The problem is formulated as a transmission problem (TP) for a second order elliptic partial differential equation with variable coefficients A 2 (x, x )u(x) xk [ a (2) kj (x) x j u(x) ] + ω 2 κ(x)u(x) = f 2 (1) in the inhomogeneous anisotropic region Ω 2 and for the anisotropic Helmoltz type equation with constant coefficients A 1 ( x )u(x) xk [ a (1) kj x j u(x) ] + ω 2 u(x) = f 1 (2) in the unbounded homogeneous region Ω 1. 0-4

Since the material parameters a (q) kj and the refractive index κ are assumed to be discontinuous across the interface S = Ω 1 = Ω 2 between the inhomogeneous interior and homogeneous exterior regions, there are given standard transmission conditions relating the interior and exterior traces of the sought for wave functions and their co-normal derivatives on S. 0-5

Since the material parameters a (q) kj and the refractive index κ are assumed to be discontinuous across the interface S = Ω 1 = Ω 2 between the inhomogeneous interior and homogeneous exterior regions, there are given standard transmission conditions relating the interior and exterior traces of the sought for wave functions and their co-normal derivatives on S. In a particular case of isotropic inhomogeneity, when a (2) kj (x) = δ kj, a (1) kj = δ kj, i.e., when A 2 (x, x ) = + ω 2 κ(x), A 1 ( x ) = + ω 2, the similar transmission problems, is well investigated in the literature (e.g. [Colton, Kress](2013): Lippmann Schwinger Equation - Fredholm-Riesz type integral equation). 0-6

Another particular case of isotropic inhomogeneity, when a (2) kj (x) = a(x) δ kj, a (1) kj = δ kj, i.e., when A 2 (x, x )u(x) = xk [ a(x) xk u(x) ] + ω 2 κ(x)u(x), (3) A 1 ( x )u(x) = u(x) + ω 2 u(x), (4) by the indirect boundary-domain integral equation method is investigated by P.Werner (1960). P.Werner reduces the problem to the Fredholm-Riesz type integral equations system and proves its unique solvability. 0-7

Another particular case of isotropic inhomogeneity, when a (2) kj (x) = a(x) δ kj, a (1) kj = δ kj, i.e., when A 2 (x, x )u(x) = xk [ a(x) xk u(x) ] + ω 2 κ(x)u(x), (3) A 1 ( x )u(x) = u(x) + ω 2 u(x), (4) by the indirect boundary-domain integral equation method is investigated by P.Werner (1960). P.Werner reduces the problem to the Fredholm-Riesz type integral equations system and proves its unique solvability. The same problem for isotropic inhomogeneity by the direct method is considered by P.Martin (2003) (using essentially the existence results obtained by P.Werner). 0-8

Our goal is to consider the above described wave scattering problems for general inhomogeneous anisotropic case, applying the method of Localized Boundary-Domain Integral Equations. 0-9

FORMULATION OF THE TRANSMISSION PROBLEM Let Ω + = Ω 2 be a bounded domain in R 3 with a simply connected boundary Ω + = S, Ω + = Ω + S, and Ω 1 = Ω := R 3 \ Ω +. For simplicity, we assume that S C if not otherwise stated. n = (n 1, n 2, n 3 ) denotes the unit normal vector to S directed outward with respect to the domain Ω +. We assume that the propagation region of time harmonic acoustic waves is all of R 3 which consists of an anisotropic inhomogeneous part Ω 2 := Ω + and a anisotropic homogenous one Ω 1 := Ω. 0-10

Acoustic wave propagation in Ω 1 and Ω 2 is governed by the uniformly elliptic second order scalar partial differential equations ( ) A 1 ( x )u 1 (x) xk a (1) kj x j u 1 (x) +ω 2 u 1 (x) = 0, x Ω 1, (5) ( ) A 2 (x, x )u 2 (x) xk a (2) kj (x) x j u 2 (x) +ω 2 κ(x)u 2 (x) = = f(x), x Ω 2 ; (6) u 1, u 2 wave amplitudes, ω R wave number, κ(x) 0 refractive index, f L 2 (Ω 2 ) a given function. a (2) kj, κ C 2 (Ω 2 ), a (q) kj = a(q) jk, j, k = 1, 2, 3, q = 1, 2, (7) The matrices a (q) = [a (q) kj ] 3 3 are uniformly positive definite, i.e., there are positive constants c 1 and c 2 such that c 1 ξ 2 a (q) kj (x)ξ k ξ j c 2 ξ 2, x Ω q, ξ R 3, q = 1, 2. (8) 0-11

In the unbounded region Ω 1 = Ω we have a total wave field u = u tot = u ins + u sc, where u ins is a wave motion initiating known incident field and u 1 = u sc is a radiating unknown scattered field. 0-12

In the unbounded region Ω 1 = Ω we have a total wave field u = u tot = u ins + u sc, where u ins is a wave motion initiating known incident field and u 1 = u sc is a radiating unknown scattered field. The co-normal derivative operators read as {T 1 (x, x )v(x)} := {a (1) kj n k(x) xj v(x)}, x S, (9) {T 2 (x, x )v(x)} + := {a (2) kj (x) n k(x) xj v(x)} +, x S; (10) u ± {u} ±, T ± j u {T j u} ± one sided traces on S from Ω ±. 0-13

SOMMERFELD RADIATION CONDITIONS Denote by S ω the characteristic surface (ellipsoid) associated with the operator A 1 ( ) = a (1) kj k j, Φ 1 (ξ, ω) := a (1) kj ξ k ξ j ω 2 = 0, ξ R 3. For an arbitrary vector η R 3 with η = 1 there exists only one point ξ(η) S ω such that the outward unit normal vector n(ξ(η)) to S ω at the point ξ(η) has the same direction as η, i.e., n(ξ(η)) = η. Note that ξ( η) = ξ(η) S ω and n( ξ(η)) = η. It can be easily verified that ξ(η) = ω (bη η) 1 2 bη, (11) where b := [a (1) ] 1 is the matrix inverse to a (1) := [a (1) kj ] 3 3. 0-14

SOMMERFELD CLASS Z(Ω ) associated with the operator A 1 ( ) : A complex valued function w belongs to Z(Ω ) if for sufficiently large x 1 the following radiation conditions hold w(x) = O( x 1 ), (12) k w(x) iξ k (η)w(x) = O( x 2 ), k = 1, 2, 3, (13) where ξ(η) S ω corresponds to the vector η = x/ x. LEMMA 1. [Analogue of the Rellich-Vekua lemma, 1943] Let w be a solution of the homogeneous equation A 1 ( x )w = 0 in Ω and let lim w(x) 2 dσ R = 0, (14) R + Σ R where Σ R is the sphere centered at the origin and radius R. Then w = 0 in Ω. 0-15

Find complex valued functions TRANSMISSION PROBLEM: u 2 H 1, 0 (Ω +, A 2 ), u 1 H 1, 0 loc (Ω, A 1 ) Z(Ω ) satisfying the differential equations A 1 ( x )u 1 (x) = 0 in Ω 1, (15) A 2 (x, x )u 2 (x) = f(x) in Ω 2, (16) and the transmission conditions on the interface S where u + 2 u 1 = φ 0 on S, (17) T + 2 u 2 T 1 u 1 = ψ 0 on S, (18) φ 0 := {u inc } H 1 2 (S), ψ0 := {T 1 u inc } H 1 2 (S), (19) f H 0 (Ω + ). (20) 0-16

Basic integral relations in the bounded domain Ω 2 = Ω + Localized harmonic parametrix P (x) P χ (x) := χ(x) 4 π x, χ Xk +, k 3. (21) Green s third formula for u 2 H 1, 0 (Ω +, A 2 ) β(y) u 2 (y) + N 2 u 2 (y) V 2 (T + 2 u 2)(y)+ + W 2 (u + 2 )(y) = P 2( A2 u 2 ) (y), y Ω +, (22) β(x) = 1 3 [ a(2) 11 (x) + a(2) 22 (x) + a(2) 33 (x) ], (23) N 2 v(y) := v.p. [A 2 (x, x )P (x y)]v(x) dx, (24) Ω + 0-17

V 2, W 2 and P 2 are the localized single layer, double layer and Newtonian volume potentials V 2 (g)(y) := P (x y) g(x) ds x, (25) S [ W 2 (g)(y) := T2 (x, x ) P (x y) ] g(x) ds x, (26) S P 2 (h)(y) := P (x y) h(x) dx. (27) Ω + The principal homogeneous symbol S 0 (N 2 ; y, ξ) is rational in ξ S 0 (N 2 ; y, ξ) = a(2) kl (y) ξ k ξ l ξ 2 β(y) = A 2(y, ξ) ξ 2 β(y), (28) y Ω +, ξ R 3, If a (2) kl (y) = a(2) (y)δ kl, then S 0 (N 2 ; y, ξ) = 0 and the operator N 2 becomes weakly singular integral operator. 0-18

The interior trace of Green s third formula on S reads as N + 2 u 2 V 2 (T + 2 u 2)+(β µ+w 2 )u + 2 = P+ 2 ( A2 u 2 ) on S, (29) where µ(y) := 1 2 a(2) kj (y) n k(y) n j (y) > 0, y S, (30) V 2 g(y) := P (x y) g(x) ds x, y S, (31) S [ W 2 g(y) := T2 (x, x ) P (x y) ] g(x) ds x, y S, (32) N + 2 u 2 := {N 2 u 2 } + S, P+ 2 S ( A2 u 2 ) := {P2 ( A2 u 2 ) } +. (33) 0-19

LEMMA 2. Let Φ 2 H 1, 0 (Ω +, ), ψ 2 H 1 2 (S), φ2 H 1 2 (S). (34) Moreover, let u 2 H 1 (Ω + ) and the following equation hold in Ω + β(y) u 2 (y) + N 2 u 2 (y) V 2 (ψ 2 )(y) + W 2 (φ 2 )(y) = Φ(y). (35) Then u H 1, 0 (Ω +, A 2 ). 0-20

Basic relations in unbounded domain Ω 1 = Ω For any radiating solution u 1 H 1, 0 loc (Ω, A 1 ) Z(Ω ) of homogeneous equation A 1 ( )u 1 = 0 there holds the following Green s third formula: u 1 (y) + V 1 (T 1 u 1)(y) W 1 (u 1 )(y) = 0, y Ω, (36) where V 1 (g)(y) := γ(x y) g(x) ds x, y R 3 \ S, (37) S W 1 (g)(y) := [T 1 (x, x )γ(x y)] g(x) ds x, y R 3 \ S. (38) S γ(x, ω) is a radiating fundamental function of the operator A 1 ( x ): γ(x, ω) = exp{iω(bx x)1/2 } 4π(det a) 1/2 (bx x) 1/2, b = a 1, a = [a (1) kj ] 3 3.(39) 0-21

(i) Mapping properties Properties of radiating potentials V 1 : H 1 2 (S) H 1 loc (Ω, A 1 ) Z(Ω ), W 1 : H 1 2 (S) H 1 loc (Ω, A 1 ) Z(Ω ). (40) (ii) Jump relations for h H 1 2 (S) and g H 1 2 (S): {V 1 (h)} + = {V 1 (h)} = V 1 (h) on S, (41) {T 1 V 1 (h)} ± = (± 2 1 I + W 1 )h on S, (42) {W 1 (g)(y)} ± = ( 2 1 I + W 1 )g on S, (43) T + 1 W 1(g) = T 1 W 1(g) =: L 1 (g) on S, (44) 0-22

where I stands for the identity operator, and V 1 (h)(y) := γ(x y) h(x) ds x, y S, (45) S W 1 (g)(y) := [T 1 (x, x )γ(x y)] g(x) ds x, y S, (46) S W 1 (h)(y) := [T 1 (y, y )γ(x y)] h(x) ds x, y S. (47) S (iii) The operators V 1 : H 1 2 (S) H 1 2 (S), (48) W 1 : H 1 2 (S) H 1 2 (S), (49) W 1 : H 1 2 (S) H 1 2 (S), (50) L 1 : H 1 2 (S) H 1 2 (S), (51) are continuous. The operators (49)-(51) are compact, V 1 and L 1 are pseudodifferential operators of order 1 and 1 respectively. 0-23

Let α be a positive constant and define the operators K 1 and M 1 : K 1 := 2 1 I + W 1 i α V 1 {[T 1 iα]v 1 (g)} +, (52) M 1 := L 1 i α ( 2 1 I + W 1 ) {[T 1 iα]w 1 (g)} +. (53) For a solution u 1 H 1, 0 loc (Ω, A 1 ) Z(Ω ) of equation A 1 u 1 = 0 the following relation can be derived from Green s third formula: K 1 (T 1 u 1) M 1 (u 1 ) = 0 on S. (54) LEMMA 3. The operators K 1 : H 1 2 (S) H 1 2 (S), M1 : H 1 2 (S) H 1 2 (S), (55) are invertible. From (54) and Lemma3 the following Steklov-Poincaré relation for arbitrary u 1 H 1, 0 loc (Ω, A 1 ) Z(Ω ) follows T 1 u 1 = K 1 1 M 1 u 1 on S, (56) where K 1 1 is the inverse to the operator K 1. 0-24

Equivalent reduction to the system of integral equations Let us set φ 1 := u 1, φ 2 := u + 2, ψ 1 := T 1 u 1, ψ 2 := T + 2 u 2. (57) If a pair (u 1, u 2 ) solves the transmission problem, then in view of the notation (57), the relations obtained above and the transmission conditions can be rewritten as ( ) β u 2 + N 2 u 2 V 2 (ψ 2 ) + W 2 (φ 2 ) = P 2 f in Ω +, (58) N + 2 u 2 V 2 ψ 2 + (β µ + W 2 )φ 2 = P + ( ) 2 f on S, (59) u 1 + V 1 (ψ 1 ) W 1 (φ 1 ) = 0 in Ω, (60) K 1 ψ 1 M 1 φ 1 = 0 on S, (61) φ 2 φ 1 = φ 0 on S, (62) ψ 2 ψ 1 = ψ 0 on S. (63) 0-25

Rewrite this LBDIE system in the following equivalent form ( ) β u 2 + N 2 u 2 V 2 (ψ 2 ) + W 2 (φ 2 ) = P 2 f in Ω +, (64) N + 2 u 2 V 2 ψ 2 + (β µ + W 2 )φ 2 = P + ( ) 2 f on S, (65) K 1 ψ 2 M 1 φ 2 = K 1 ψ 0 M 1 φ 0 on S, (66) φ 1 = φ 2 φ 0 on S, (67) ψ 1 = ψ 2 ψ 0 on S, (68) u 1 + V 1 (ψ 1 ) W 1 (φ 1 ) = 0 in Ω. (69) Let us consider these relations as a system of equations with respect to unknowns where (u 2, ψ 2, φ 2, ψ 1, φ 1, u 1 ) H, (70) H := H 1,0 (Ω +, A 2 ) H 1 2 (S) H 1 2 (S) H 1 2 (S) H 1 2 (S) ( H 1,0 loc (Ω, A 1 ) Z(Ω ) ). (71) 0-26

EQUIVALENCE THEOREM The transmission problem (15)-(20) and the system of integral equations (64)-(69) are equivalent in the following sense: (i) If a pair (u 2, u 1 ) H 1, 0 (Ω +, A 2 ) ( H 1, 0 loc (Ω, A 1 ) Z(Ω ) ) solves the transmission problem (15)-(20), then the vector (u 2, ψ 2, φ 2, ψ 1, φ 1, u 1 ) H with ψ 2, φ 2, ψ 1, φ 1, defined by the equalities φ 1 := u 1, φ 2 := u + 2, ψ 1 := T 1 u 1, ψ 2 := T + 2 u 2, (72) solves the system (64)-(69), and vice versa, (ii) If a vector (u 2, ψ 2, φ 2, ψ 1, φ 1, u 1 ) H solves the system (64)- (69), then the pair (u 2, u 1 ) H 1, 0 (Ω +, A 2 ) ( H 1, 0 loc (Ω, A 1 ) Z(Ω ) ) solves the transmission problem (15)-(20) and the relations (72) hold true. 0-27

Proof. (i) The first part of the theorem is trivial. (ii) Let a vector (u 2, ψ 2, φ 2, ψ 1, φ 1, u 1 ) H solve the system (64)- (69). Taking the trace of (64) on S and comparing with (65) lead to the equation u + 2 = φ 2 on S. (73) Further, since u 2 H 1, 0 (Ω 2, A 2 ) we can write Green s third identity (22) which in view of (73) can be rewritten as β(y) u 2 (y) + N 2 u 2 (y) V 2 (T + 2 u 2)(y) + W 2 (φ 2 )(y) = From (64) and (74) it follows that implying = P 2 ( A2 u 2 ) (y), y Ω2. (74) V 2 (T + 2 u 2 ψ 2 ) + P 2 ( A2 u 2 f 2 ) = 0 in Ω2, (75) A 2 u 2 = f 2 in Ω 2, T + 2 u 2 = ψ 2 on S. (76) 0-28

Further, from equation (69) it follows that A 1 u 1 = 0 in Ω 1, K 1 ψ 1 M 1 φ 1 = 0 on S. (77) Now, let us consider the function It can be shown that v := V 1 (ψ 1 ) W 1 (φ 1 ) in Ω 2 = Ω +. (78) A 1 v = 0 in Ω 1, (79) T + 1 v i α v+ = 0 on S, α R \ {0}, (80) which implies that v vanishes identically in Ω 2 = Ω +, since the Robin problem (79)-(80) possess only the trivial solution. Therefore u 1 = u 1 + v+ = φ 1, T 1 u 1 = T 1 u 1 + T + 1 v = ψ 1. (81) which completes the proof. 0-29

Uniqueness Theorem. The homogeneous transmission problem possesses only the trivial solution. Proof. Let B(R) be a ball centered at the origin and radius R, such that Ω + B(R). Let a pair (u 1, u 2 ) be a solution to the homogeneous transmission problem. Write Green s formulas for the domains Ω + and Ω R = Ω B(R) [a (2) kj (x) ju 2 (x) k u 2 (x) ω 2 κ(x) u 2 (x) 2 ]dx = T + 2 u 2, u + 2 S, Ω + Ω R [a (1) kj ju 1 (x) k u 1 (x) ω 2 u 1 (x) 2 ] dx = T 1 u 1, u 1 S+ + T 1 u 1, u 1 Σ(R). Since the matrix a = [a kj ] 3 3 is symmetric and positive definite, in view of the homogeneous transmission conditions we get { } Im T 1 (x, )u 1 (x) u 1 (x) dσ R = 0. (82) Σ R 0-30

Since u 1 Z(Ω ) we get u 1 (x)t 1 (x, )u 1 (x) = i ω (bη η) 1 2 u1 (x) 2 + O( x 3 ). (83) Due to positive definiteness of the matrix b := [a (1) ] 1, there are positive constants δ 0 and δ 1 such that for all η Σ 1 0 < δ 0 (bη η) 1 2 δ1 <. (84) Therefore (82) lim R + u 1 (x) 2 dσ R = 0 Σ R (85) and by Rellich-Vekua Lemma u 1 = 0 in Ω. 0-31

Consequently, the function u 2 solves the homogeneous Cauchy problem in Ω + for the elliptic partial differential equation A 2 (x, )u 2 = 0 in Ω +, {u 2 } + = 0, {T 2 u 2 } + = 0 on S, with variable coefficients a (2) kj and κ 2 with a (2) kj, κ 2 C 2 (Ω 2 ). By the interior and boundary regularity properties of solutions to elliptic equations we have that u 2 C 2 (Ω + ) and therefore u 2 = 0 in Ω + due to the well known uniqueness theorem for the Cauchy problem (E.M. Landis - 1956, A.P.Calderon - 1958). 0-32

Consequently, the function u 2 solves the homogeneous Cauchy problem in Ω + for the elliptic partial differential equation A 2 (x, )u 2 = 0 in Ω +, {u 2 } + = 0, {T 2 u 2 } + = 0 on S, with variable coefficients a (2) kj and κ 2 with a (2) kj, κ 2 C 2 (Ω 2 ). By the interior and boundary regularity properties of solutions to elliptic equations we have that u 2 C 2 (Ω + ) and therefore u 2 = 0 in Ω + due to the well known uniqueness theorem for the Cauchy problem (E.M. Landis - 1956, A.P.Calderon - 1958). COROLLARY 4. The LBDIE system (64)-(69) possesses at most one solution in the space H. 0-33

Consequently, the function u 2 solves the homogeneous Cauchy problem in Ω + for the elliptic partial differential equation A 2 (x, )u 2 = 0 in Ω +, {u 2 } + = 0, {T 2 u 2 } + = 0 on S, with variable coefficients a (2) kj and κ 2 with a (2) kj, κ 2 C 2 (Ω 2 ). By the interior and boundary regularity properties of solutions to elliptic equations we have that u 2 C 2 (Ω + ) and therefore u 2 = 0 in Ω + due to the well known uniqueness theorem for the Cauchy problem (E.M. Landis - 1956, A.P.Calderon - 1958). COROLLARY 4. The LBDIE system (64)-(69) possesses at most one solution in the space H. REMARK 5: Due to the recent results, the uniqueness theorem for the Cauchy problem for a scalar elliptic operator holds if variable coefficients are Lipschitz continuous and Ω + is a Dini domain (X.X.Tao and S.Y.Zhang - 2007). 0-34

Investigation of the LBDIE system ( ) (β I + N 2 ) E 0 u 2 V 2 (ψ 2 ) + W 2 (φ 2 ) = P 2 f2 in Ω +, (86) N + 2 E 0 u 2 V 2 ψ 2 + (β µ + W 2 )φ 2 = P + ( ) 2 f2 on S, (87) K 1 ψ 2 M 1 φ 2 = K 1 ψ 0 M 1 φ 0 on S, (88) φ 1 = φ 2 φ 0 on S, (89) ψ 1 = ψ 2 ψ 0 on S, (90) u 1 + V 1 (ψ 1 ) W 1 (φ 1 ) = 0 in Ω. (91) where E 0 denotes the extension operator by zero from Ω + onto Ω. We need to investigate the matrix operator generated by the left hand side expressions in the first three equations. If the unknowns u 2, ψ 2, and φ 2 are found from the first three equations, then the unknowns ψ 1, φ 1,u 1 can be defined explicitly from the last three equations. 0-35

Let us rewrite the first three equations of the LBDIE system (86)-(91) in matrix form M := M U (2) = F, (92) r (β I + N Ω + 2)E 0 r V Ω + 2 r W Ω + 2 A + 2 E 0 V 2 (β µ)i + W 2 (93) 0 K 1 M 1 U (2) := (u 2, ψ 2, φ 2 ), F := ( F 1, F 2, F 3 ), (94) Applying the mapping properties of the layer potentials and pseudodifferential operators with rational type symbols we deduce that the following operator is continuous M : X 1 Y 1, (95) X 1 := H 1 (Ω + ) H 1 2 (S) H 1 2 (S), (96) Y 1 := H 1 (Ω + ) H 1 2 (S) H 1 2 (S). (97) 0-36

THEOREM 4. The operator is invertible. M : X 1 Y 1 0-37

THEOREM 5. The operator is invertible. M : X 1 Y 1 Proof. Invertibility of the operator M can be shown with the help of the Vishik-Eskin theory. The proof is performed into several steps. Step 1. The operator M 11 = r Ω + (β I + N 2)E 0 : H 1 (Ω + ) H 1 (Ω + ) is Fredholm with zero index. The proof is based on the fact that the principal homogeneous symbol S 0 (D 11 ; y, ξ) = a(2) kl (y)ξ kξ l ξ 2 > 0, y Ω +, ξ R 3 \ {0}, is positive, rational, homogeneous of order 0 in ξ. 0-38

Step 2. The operator M is Fredholm with zero index. (i) Fredholm property follows from the fact that for the operator M the generalized Šapiro-Lopatinskĭi condition holds which is verified by direct calculations of specific Cauchy type integrals. (ii) The zero index is established by showing that M is homotopic to an invertible operator. 0-39

Step 2. The operator M is Fredholm with zero index. (i) Fredholm property follows from the fact that for the operator M the generalized Šapiro-Lopatinskĭi condition holds which is verified by direct calculations of specific Cauchy type integrals. (ii) The zero index is established by showing that M is homotopic to an invertible operator. Step 3. The null space of the operator M is trivial. This is shown with the help of the Equivalence and Uniqueness Theorems. Consequently, the operator M is invertible. 0-40

COROLLARY 6. The operator M : H 1, 0 (Ω +, A 2 ) H 1 2 (S) H 1 2 (S) is invertible. H 1, 0 (Ω +, ) H 1 2 (S) H 1 2 (S) (98) COROLLARY 7. Let a cut-off function χ X+ Then the operator and r 0. M : H r+1 (Ω + ) H r 1 2 (S) H r+ 1 2 (S) is invertible. H r+1 (Ω + ) H r+ 1 2 (S) H r 1 2 (S) (99) 0-41

From the above results it follows that the LBDIE system (β I + N 2 ) E 0 u 2 V 2 (ψ 2 ) + W 2 (φ 2 ) = h 1 in Ω +, (100) A + 2 E 0 u 2 V 2 ψ 2 + (β µ + W 2 )φ 2 = h 2 on S, (101) K 1 ψ 2 M 1 φ 2 = h 3 on S, (102) ψ 2 ψ 1 = h 4 on S, (103) φ 2 φ 1 = h 5 on S, (104) u 1 + V 1 (ψ 1 ) W 1 (φ 1 ) = h 6 in Ω. (105) is uniquely solvable in the space X := H 1 (Ω + ) H 1 1 2 (S) H 2 (S) H 1 2 (S) ( H 1 2 (S) H 1 loc (Ω ) Z(Ω ) ) for arbitrary right hand side data (h 1,, h 6 ) Y with Y := H 1 (Ω + ) H 1 2 (S) H 1 2 (S) H 1 2 (S) H 1 2 (S) H 1 comp (Ω ). 0-42

THE PROBLEMS OF ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES are uniquely solvable. 0-43

THANK YOU! 0-44