Zum Hauptinhalt springen

Partial data inverse problems for quasilinear conductivity equations.

Kian, Yavar ; Krupchyk, Katya ; et al.
In: Mathematische Annalen, Jg. 385 (2023-04-01), Heft 3/4, S. 1611-1638
Online academicJournal

Partial data inverse problems for quasilinear conductivity equations  Introduction and statement of results

We show that the knowledge of the Dirichlet-to-Neumann maps given on an arbitrary open non-empty portion of the boundary of a smooth domain in R n , n ≥ 2 , for classes of semilinear and quasilinear conductivity equations, determines the nonlinear conductivities uniquely. The main ingredient in the proof is a certain L 1 -density result involving sums of products of gradients of harmonic functions which vanish on a closed proper subset of the boundary.

Let ΩRn , n2 , be a connected bounded open set with C boundary. Let us consider the Dirichlet problem for the following isotropic semilinear conductivity equation,

1.1 div(γ(x,u)u)=0inΩ,u=fonΩ.

Graph

Here we assume that the function γ:Ω¯×CC satisfies the following conditions,

  • the map Cτγ(·,τ) is holomorphic with values in the Hölder space C1,α(Ω¯) with some 0<α<1 ,
  • γ(x,0)=1 , for all xΩ .

The semilinear conductivity equation (1.1) can be viewed as a steady state semilinear heat equation where the conductivity depends on the temperature, and in physics, such models occur, for instance, in nonlinear heat conduction in composite materials, see [[23]].

It is shown in Theorem B.1 that under the assumptions (a) and (b), there exist δ>0 and C>0 such that when fBδ(Ω):={fC2,α(Ω):fC2,α(Ω)<δ} , the problem (1.1) has a unique solution u=ufC2,α(Ω¯) satisfying uC2,α(Ω¯)<Cδ . Let ΓΩ be an arbitrary non-empty open subset of the boundary Ω . Associated to the problem (1.1), we define the partial Dirichlet-to-Neumann map

ΛγΓ(f)=(γ(x,u)νu)|Γ,

Graph

where fBδ(Ω) with supp(f)Γ . Here ν is the unit outer normal to the boundary.

We are interested in the following inverse boundary problem for the semilinear conductivity equation (1.1): given the knowledge of the partial Dirichlet-to-Neumann map ΛγΓ , determine the semilinear conductivity γ in Ω¯×C . Our first main result gives a complete solution to this problem.

Theorem 1.1

Let ΩRn , n2 , be a connected bounded open set with C boundary, and let ΓΩ be an arbitrary open non-empty subset of the boundary Ω . Let γ1,γ2:Ω¯×CC satisfy the assumptions (a) and (b). If Λγ1Γ=Λγ2Γ then γ1=γ2 in Ω¯×C .

It is also of great interest and importance to consider inverse boundary problems for nonlinear conductivity equations with conductivities depending not only on the solution u but also on its gradient, u . Such equations occur, in particular, in the study of transport properties of non-linear composite materials, see [[36]], as well as in glaciology, when modeling the stationary motion of a glacier, see [[12]]. Furthermore, such equations can be considered as a simple scalar model of the nonlinear elasticity system, see [[44], Section 2]. To this end, we are able to solve partial data inverse boundary problems for a class of quasilinear conductivities of the form γ(x,u,ω·u) , depending on the space variable, the solution, as well as the derivative of the solution in a fixed direction ωSn-1={ωRn:|ω|=1} . To state the result, let ωSn-1 be fixed and let us consider the Dirichlet problem for the following isotropic quasilinear conductivity equation,

1.2 div(γ(x,u,ω·u)u)=0inΩ,u=λ+fonΩ.

Graph

Here we assume that the function γ:Ω¯×C×CC satisfies the following conditions,

  • the map C×C(τ,z)γ(·,τ,z) is holomorphic with values in C1,α(Ω¯) with some 0<α<1 ,
  • γ(x,τ,0)=1 , for all xΩ and all τC .

It is established in Theorem B.1 that under the assumptions (i) and (ii) for each λC , there exist δλ>0 and Cλ>0 such that when fBδλ(Ω):={fC2,α(Ω):fC2,α(Ω)<δλ} , the problem (1.2) has a unique solution u=uλ,fC2,α(Ω¯) satisfying u-λC2,α(Ω¯)<Cλδλ . Associated to the problem (1.2), we define the partial Dirichlet-to-Neumann map

ΛγΓ(λ+f)=(γ(x,u,ω·u)νu)|Γ,

Graph

where fBδλ(Ω) with supp(f)Γ and λC .

Our second main result is as follows.

Theorem 1.2

Let ΩRn , n2 , be a connected bounded open set with C boundary, and let ΓΩ be an arbitrary open non-empty subset of the boundary Ω . Let ωSn-1 be fixed. Assume that γ1,γ2:Ω¯×C×CC satisfy the assumptions (i) and (ii). Let ΣC be a set which has a limit point in C . Then if for all λΣ , we have

Λγ1Γ(λ+f)=Λγ2Γ(λ+f),fBδλ(Ω),supp(f)Γ,

Graph

then γ1=γ2 in Ω¯×C×C .

Note that in Theorem 1.2 the Dirichlet-to-Neumann maps ΛγjΓ map the Dirichlet data λ+f , which is not supported on Γ , unless λ=0 , to the Neumann data which is measured on Γ .

Remark 1.3

To the best of our knowledge, the partial data results of Theorem 1.1 and Theorem 1.2 are the first partial data results for nonlinear conductivity equations.

Remark 1.4

It might be interesting to note that an analog of the partial data results of Theorem 1.1 and Theorem 1.2 is still not known in the case of the linear conductivity equation in dimensions n3 . We refer to [[17]] for the corresponding partial data result for the linear conductivity equation in dimension n=2 .

Remark 1.5

An analog of Theorem 1.1 in the full data case, i.e. when Γ=Ω , was proved in [[42]] where instead of working with small Dirichlet data one considers small perturbations of constant Dirichlet data as in (1.2). Furthermore, it was assumed in [[42]] that the semilinear conductivity is strictly positive while no analyticity was required. The proof of [[42]] relies on a first order linearization of the Dirichlet-to-Neumann map at constant Dirichlet boundary values which leads to the inverse boundary problem for the linear conductivity equation and therefore, an application of results of [[47]] and [[35]] for the linear conductivity problem in dimensions n3 and in dimension n=2 , respectively, gives the recovery of the semilinear conductivity.

Remark 1.6

To the best of our knowledge Theorem 1.2 is new even in the full data case. Indeed, in the full data case, so far authors have only considered the recovery of conductivities of the form γ(x,u) , see e.g. [[42], [46]], or of the form γ(u,u) , see e.g. [[34], [41]], or conductivities which depend x and u in some specific way, see e.g. [[5]]. We obtain in Theorem 1.2, for what seems to be the first time, the recovery of some general class of quasilinear conductivities of the form γ(x,u,ω·u) , depending on the space variable, the solution, as well as the derivative of the solution in a fixed direction.

Remark 1.7

The assumption that the conductivity is holomorphic as a function Cτγ(·,τ,·) in Theorem 1.2 is motivated by the proof of the solvability of the forward problem and the differentiability with respect to the boundary data. This assumption could perhaps be weakened and one could show that the full knowledge of the partial Dirichlet-to-Neumann map ΛγΓ determines the conductivity γ . As the main focus of this paper is on establishing the partial data inverse results, we decided not to elaborate upon this issue further.

We remark that starting with [[27]], it has been known that nonlinearity may be helpful when solving inverse problems for hyperbolic PDE. Analogous phenomena for nonlinear elliptic equations have been revealed and exploited in [[10], [29]], see also [[24]–[26], [28], [30]]. A noteworthy aspect of all of these works is that the presence of a nonlinearity enables one to solve inverse problems for nonlinear PDE in situations where the corresponding inverse problems for linear equations are still open. The present paper is also concerned with illustrating this general phenomenon.

Let us proceed to discuss the main ideas of the proofs of Theorem 1.1 and Theorem 1.2. Using the technique of higher order linearizations of the partial Dirichlet-to-Neumann map, introduced in [[10], [29]], see also [[42], [46]] for the use of the second linearization, we reduce the proof of Theorem 1.2 to the following density result.

Theorem 1.8

Let ΩRn , n2 , be a connected bounded open set with C boundary, let ΓΩ be an open non-empty subset of Ω , let ωSn-1 be fixed, and let m=2,3,, be fixed. Let fL(Ω) be such that

1.3 Ωf(k=1mr=1,rkm(ω·ur)uk)·um+1dx=0,

Graph

for all functions ulC(Ω¯) harmonic in Ω with supp(ul|Ω)Γ , l=1,,m+1 . Then f=0 in Ω .

Similarly, using higher order linearizations of the partial Dirichlet-to-Neumann map, we show that Theorem 1.1 will follow from the following density result.

Theorem 1.9

Let ΩRn , n2 , be a connected bounded open set with C boundary, let ΓΩ be an open non-empty subset of Ω , and let m=2,3,, be fixed. Let fL(Ω) be such that

1.4 Ωf(k=1mr=1,rkmuruk)·um+1dx=0,

Graph

for all functions ulC(Ω¯) harmonic in Ω with supp(ul|Ω)Γ , l=1,,m+1 . Then f=0 in Ω .

Theorems 1.8 and 1.9 can be viewed as extensions of the results of [[8]] and [[24]]. Indeed, it was proved in [[8]] that the linear span of the set of products of harmonic functions in Ω which vanish on a closed proper subset of the boundary is dense in L1(Ω) , and this density result was extended in [[24]] by showing that the linear span of the set of scalar products of gradients of harmonic functions in Ω which vanish on a closed proper subset of the boundary is also dense in L1(Ω) .

To prove Theorem 1.8, we shall follow the general strategy of the work [[8]], see also [[24]]. We first establish a corresponding local result in a neighborhood of a boundary point in Γ assuming, as we may, that Γ is a small open neighborhood of this point, see Proposition 2.1 below. We then show how to pass from this local result to the global one of Theorem 1.8. The essential difference here compared with the works [[8], [24]] is that working with products of m+1 gradients in the orthogonality identity (1.3), we need to prove a certain Runge type approximation theorem in the W1,m+1 -topology for any m=2,3, fixed, as opposed to L2 and H1 approximation results obtained in [[8]] and [[24]], respectively.

We shall only prove Theorem 1.8 as the proof of Theorem 1.9 is obtained by inspection of that proof as the only difference between the orthogonality relations (1.3) and (1.4) is that (1.3) contains ω·ur with harmonic functions ur while (1.4) contains ur instead, and no new difficulties occur.

Remark 1.10

While the present paper was under review, the inverse boundary problem with full data, i.e. when the measurement are performed along the entire boundary Ω , was solved in [[6]] for quasilinear isotropic conductivity γ of the form γ(x,u,u) , showing that the quasilinear conductivity γ can indeed be uniquely determined from these measurements, provided that the map C×Cn(ρ,μ)γ(·,ρ,μ) is is holomorphic with values in C1,α(Ω¯) with some 0<α<1 , and 0<γ(·,0,0)C(Ω¯) . It would be interesting to solve the partial data inverse problem for such conductivities to be on par with the full data result of [[6]]. The difficulty here compared with the recovery of the conductivities of the form γ(x,u,ω·u) in Theorem 1.2 is that higher order linearizations of the partial Dirichlet-to-Neumann map lead to a density statement in the spirit of Theorem 1.8 where instead of working with a scalar function f one has to work with a function with values in the space of symmetric tensors of rank mN . Furthermore, a challenge in the proof of partial data result compared with the full data result of [[6]] is that one has to work with harmonic functions which vanish on an arbitrary portion of the boundary in the density statement. It is not quite clear how to extend the analytic microlocal analysis framework of [[8]] to prove the needed density result in this more general situation.

Let us finally remark that inverse boundary problems for nonlinear elliptic PDE have been studied extensively in the literature. We refer to [[4], [7], [10], [15], [18]–[22], [26], [29], [34], [41]–[43], [45]], and the references given there. In particular, inverse boundary problems with partial data were studied for a certain class of semilinear equations of the form -Δu+V(x,u)=0 in [[25], [30]] relying on the density result of [[8]], for semilinear equations of the form -Δu+q(x)(u)2=0 in [[24]], and for nonlinear magnetic Schrödinger equations in [[28]].

The paper is organized as follows. In Sect. 2 we establish Theorem 1.8. Theorem 1.2 in proven in Sect. 3. The proof of Theorem 1.1 occupies Sect. 4. In Appendix A we present an alternative simple proof of Theorem 1.2 in the full data case. In Appendix B we show the well-posedness of the Dirichlet problem for our quasilinear conductivity equation, in the case of boundary data close to a constant one.

Proof of Theorem 1.8

We shall proceed by following the general strategy of [[8]]. It suffices to assume that ΓΩ is a proper open nonempty subset of Ω , and even a small open neighborhood of some boundary point.

Local result

Theorem 1.8 will be obtained as a corollary of the following local result.

Proposition 2.1

Let ΩRn , n2 , be a bounded open set with C boundary, and let m=2,3,, be fixed. Let x0Ω , and let Γ~Ω be the complement of an open boundary neighborhood of x0 . Then there exists δ>0 such that if we have (1.3) for any harmonic functions ulC(Ω¯) satisfying ul|Γ~=0 , l=1,,m+1 , then f=0 in B(x0,δ)Ω .

Proof

It suffices to choose u1==um in (1.3). Hence, (1.3) implies that

2.1 Ωf(ω·v1)m-1v1·v2dx=0,

Graph

for all harmonic functions v1,v2C(Ω¯) satisfying vl|Γ~=0 , l=1,2 . Our goal is to show that (2.1) gives that f=0 in B(x0,δ)Ω with δ>0 . Using conformal transformations (in particular Kelvin transforms) of harmonic functions as in [[8], Section 3], and arguing as in that work, we are reduced to the following setting: x0=0 , the tangent plane to Ω at x0 is given by x1=0 ,

Ω{xRn:|x+e1|<1},Γ~={xΩ:x1-2c},e1=(1,0,,0),

Graph

for some c>0 .

Let p(ζ)=ζ2 , ζCn , be the principal symbol of -Δ , holomorphically extended to Cn . Let ζp-1(0) and let χC0(Rn) be such that supp(χ){xRn:x1-c} and χ=1 on {xΩ:x1-2c} . We shall work with harmonic functions of the form

2.2 v(x,ζ)=e-ihx·ζ+r(x,ζ),

Graph

where r is the solution to the Dirichlet problem,

-Δr=0inΩ,r|Ω=-(e-ihx·ζχ)|Ω.

Graph

By the boundary elliptic regularity, we have vC(Ω¯) , and furthermore v|Γ~=0 . Since in view of (2.1) we shall work with products of m+1 gradients of harmonic functions, we need to have good estimates for the remainder r in C1(Ω¯) . To that end, in view of Sobolev's embedding, we would like to bound rHk(Ω) with kN , k>n/2+1 . Boundary elliptic regularity gives that for k2 ,

2.3 rHk(Ω)Ce-ihx·ζχHk-1/2(Ω),

Graph

see [[9], Section 24.2]. Now by interpolation, we get

2.4 e-ihx·ζχHk-1/2(Ω)e-ihx·ζχHk(Ω)1/2e-ihx·ζχHk-1(Ω)1/2,

Graph

see [[14], Theorem 7.22, p. 189]. We have

e-ihx·ζχL2(Ω)Ce1hsupxKx·Imζ,

Graph

where K=suppχΩ , and therefore,

2.5 e-ihx·ζχHk(Ω)C(1+|ζ|h++|ζ|khk)e1hsupxKx·Imζ.

Graph

It follows from (2.4) and (2.5) that

2.6 e-ihx·ζχHk-1/2(Ω)C(1+|ζ|khk)e1hsupxKx·Imζ.

Graph

Using (2.3) and (2.6), we see that

rHk(Ω)C(1+|ζ|khk)e1hsupxKx·Imζ.

Graph

Taking k>n/2+1 and using the Sobolev embedding Hk(Ω)C1(Ω¯) , we get

2.7 rC1(Ω¯)C(1+|ζ|khk)e1hsupxKx·Imζ.

Graph

Using that supp(χ){xRn:x1-c} and χ=1 on {xΩ:x1-2c} , we obtain from (2.7) that

2.8 rC1(Ω¯)C(1+|ζ|khk)e-chImζ1e1h|Imζ|,

Graph

when Imζ10 .

Now the identity (2.1) implies that

2.9 Ωf(x)(ω·hDv(x,ζ))m-1hDv(x,ζ)·hDv(x,mη)dx=0,

Graph

for all ζ,ηp-1(0) . Here v(x,ζ) and v(x,mη) are harmonic functions of the form (2.2) and D=i-1 . Using that

(ω·hDv(x,ζ))m-1=(-ω·ζe-ihx·ζ+ω·hDr(x,ζ))m-1=(-ω·ζ)m-1e-(m-1)ihx·ζ+l=1m-1m-1l(ω·hDr(x,ζ))l(-ω·ζe-ihx·ζ)m-1-l,

Graph

we obtain from (2.9) that

2.10 Ωf(x)(-ω·ζ)m-1m(ζ·η)e-mihx·(ζ+η)dx=I1+I2,

Graph

where

I1=-Ωf(x)(-ω·ζ)m-1e-(m-1)ihx·ζ(-ζe-ihx·ζ·hDr(x,mη)-mηe-mihx·η·hDr(x,ζ)+hDr(x,ζ)·hDr(x,mη))dx,I2=-Ωf(x)l=1m-1m-1l(ω·hDr(x,ζ))l(-ω·ζe-ihx·ζ)m-1-l(mζ·ηe-ihx·(ζ+mη)-ζe-ihx·ζ·hDr(x,mη)-mηe-mihx·η·hDr(x,ζ)+hDr(x,ζ)·hDr(x,mη))dx.

Graph

We shall next proceed to bound the absolute values of I1 and I2 . To that end, first note that when Imζ10 , using the fact that Ω{xRn:|x+e1|<1} , we have

2.11

Graph

Using (2.8) and (2.11), we obtain that for all ζ,ηp-1(0) , Imζ10 , Imη10 ,

2.12 |I1|CfLem(|Imζ|+|Imη|)he-chmin(Imζ1,Imη1)|ζ|m-1(|ζ|h(1+|mη|khk)+m|η|h(1+|ζ|khk)+h2(1+|mη|khk)(1+|ζ|khk)),

Graph

and

2.13 |I2|CfLem(|Imζ|+|Imη|)he-chmin(Imζ1,Imη1)h(1+|ζ|khk)m-1(1+|ζ|m-2)(m|ζ||η|+|ζ|h(1+|mη|khk)+m|η|h(1+|ζ|khk)+h2(1+|mη|khk)(1+|ζ|khk)).

Graph

As noticed in [[8]], the differential of the map

s:p-1(0)×p-1(0)Cn,(ζ,η)ζ+η.

Graph

at a point (ζ0,η0) is surjective, provided that ζ0 and η0 are linearly independent. The latter holds if ζ0=γ and η0=-γ¯ with γCn given as follows. Recall that ω=(ω1,,ωn)Sn-1 is fixed. Then there exists ωk0 , and if 2kn we set γ=(i,0,,0,1,0,,0) where 1 is on the kth position. If ω10 then we set γ=(i,1,0,,0)Cn .

Note that γ·ω0 and ζ0+η0=2ie1 . An application of the inverse function theorem gives that there exists ε>0 small such that any zCn , |z-2ie1|<2ε , may be decomposed as z=ζ+η where ζ,ηp-1(0) , |ζ-γ|<C1ε and |η+γ¯|<C1ε with some C1>0 . We obtain that any zCn such that |z-2iae1|<2εa for some a>0 , may be decomposed as

2.14 z=ζ+η,ζ,ηp-1(0),|ζ-aγ|<C1aε,|η+aγ¯|<C1aε.

Graph

It follows from (2.14) that

2.15 |Imζ|<C1aε,|Imη|<C1aε,|ζ|Ca,|η|Ca.

Graph

We also conclude from (2.14) that for ε>0 small enough,

2.16 Imζ1>a/2,Imη1>a/2,|ζ·η|a2,|ω·ζ|>a2ω12+ωk2.

Graph

Hence, assuming that a>1 , we obtain from (2.10) with the help of (2.12), (2.13), (2.14), (2.15), (2.16) that

2.17 |Ωf(x)e-mihx·zdx|CfLe-ca2he2mC1aεh(ah)NCfLe-ca4he2mC1aεh,

Graph

for all zCn such that |z-2iae1|<2εa and ε>0 sufficiently small. Here N is a fixed integer which depends on k and m. The estimate (2.17) is completely analogous to the bound (3.8) in [[8]], and hence, the proof of Proposition 2.1 is completed by repeating the arguments of [[8]] exactly as they stand. The idea is to extrapolate the exponential decay to more values of the frequency variable z which is achieved in [[8]] by using a variant of the proof of the Watermelon theorem.

Next in order to pass from this local result to the global one of Theorem 1.8, we need a Runge type approximation theorem in the W1,m+1 -topology, m=2,3, , which will extend [[8], Lemma 2.2] and [[24], Lemma 2.2], where approximations in the L2 and H1 topologies were established, respectively. To prove such an approximation theorem, we need to recall some facts about Lp based Sobolev spaces which we shall now proceed to do.

Some facts about Lp based Sobolev spaces

Let ΩRn , n2 , be a bounded open set with C boundary, and let 1<p< . Then we have for the dual space of the Sobolev space W1,p(Ω) ,

(W1,p(Ω))=W~-1,p(Ω),

Graph

where

W~-1,p(Ω)={uW-1,p(Rn):supp(u)Ω¯},

Graph

and 1p+1p=1 , see [[3], page 163], [[38], Section 4.3.2]. The duality pairing is defined as follows: if vW~-1,p(Ω) and uW1,p(Ω) , we set

2.18 (v,u)W~-1,p(Ω),W1,p(Ω):=(v,Ext(u))W-1,p(Rn),W1,p(Rn),

Graph

where Ext(u)W1,p(Rn) is an arbitrary extension of u, see [[2], Theorem 9.7] for the existence of such an extension, and (·,·)W-1,p(Rn),W1,p(Rn) is the extension of L2 scalar product (φ,ψ)L2(Rn)=Rnφ(x)ψ(x)¯dx . One can show that the definition (2.18) is independent of the choice of an extension.

We shall also need the following fact, see [[38], Section 4.3.2, p. 318].

Proposition 2.2

C0(Ω) is dense in W~-1,p(Ω) with respect to W-1,p(Rn) topology.

We have the following result concerning the solvability of the Dirichlet problem for the Laplacian, see [[32], Theorem 7.10.2, p. 494].

Theorem 2.3

Let vW-1,p(Ω) and gW1-1/p,p(Ω) with 1<p< . Then the Dirichlet problem

-Δu=vinΩ,u|Ω=g,

Graph

has a unique solution uW1,p(Ω) . Moreover,

uW1,p(Ω)C(vW-1,p(Ω)+gW1-1/p,p(Ω)).

Graph

We shall also need the following result about the structure of distributions in W-1,p(Rn) supported by a smooth hypersurface in Rn . We refer to [[1], Theorem 5.1.13], [[31], Lemma 3.39] for this result in the case of distributions in H-1(Rn) . Since we did not find a reference for the case of distributions in W-1,p(Rn) with 1<p< , we shall present the proof of this result here.

Proposition 2.4

Let F be a smooth compact hypersurface in Rn . Let uW-1,p(Rn) , with some 1<p< , be such that supp(u)F . Then

u=vδF,v(W1-1/p,p(F))=Bp,p-(1-1/p)(F).

Graph

Here 1p+1p=1 and Bp,p-(1-1/p)(F) is the Besov space on the manifold F, see [[38], Section 2.3.1, p. 169], [[39]] for the definition, and for any φC0(Rn) , u(φ)=(vδF)(φ)=v(φ|F) .

Proof

Introducing a partition of unity and making a smooth change of variables, we see that it suffices to establish the following local result: let uW-1,p(Rn) , 1<p< , such that supp(u){xn=0} , then u=vδxn=0 with v(W1-1/p,p(Rn-1))=Bp,p-(1-1/p)(Rn-1) . In order to prove this result we follow [[31], Lemma 3.39].

First we claim that if φC0(Rn) is such that φ|xn=0=0 then u(φ)=0 . To that end, we let

φ±(x)=φ(x),ifxR±n={xRn:±xn>0},0,otherwise.

Graph

Then φ±W1,p(Rn) and therefore, by [[2], Proposition 9.18], φ±W01,p(R±n) . Thus, there exist sequences φj,±C0(R±n) such that φj,±φ± in W1,p(R±n) as j . Letting

χj(x)=φj,+(x),ifxR+n,φj,-(x),ifxR-n,

Graph

we see that χjC0(Rn) , χj=0 near {xn=0} , and χjφ in W1,p(Rn) . Hence, we have 0=u(χj)u(φ) , and therefore, u(φ)=0 , establishing the claim.

To proceed we need the following result, see [[33], [13], Theorem 1.5.1.1, p. 37]. The trace operator uu|xn=0 , which is defined on C0(Rn) , has a unique continuous extension as an operator,

γ:W1,p(Rn)W1-1/p,p(Rn-1),1<p<.

Graph

This operator has a right continuous inverse, the extension operator,

E:W1-1/p,p(Rn-1)W1,p(Rn)

Graph

so that γ(Eψ)=ψ for all ψW1-1/p,p(Rn-1) .

Now we define

2.19 v(φ)=u(Eφ),φC0(Rn-1).

Graph

We have

|v(φ)|uW-1,p(Rn)EφW1,p(Rn)CuW-1,p(Rn)φW1-1/p,p(Rn-1),

Graph

and therefore, v(W1-1/p,p(Rn-1)) . Note that when 1<p< ,

W1-1/p,p(Rn-1)=Bp,p1-1/p(Rn-1),(Bp,p1-1/p(Rn-1))=Bp,p-(1-1/p)(Rn-1),

Graph

see [[38], Section 2.5, p. 190, and Section 2.6.1, p. 198].

Finally, we claim that u-vδxn=0=0 . Indeed, letting φC0(Rn) and using (2.19) and our first claim, we get

(u-vδxn=0)(φ)=u(φ)-v(φ|xn=0)=u(φ-E(φ|xn=0))=0.

Graph

This completes the proof of Proposition 2.4.

Runge type approximation

Let Ω1Ω2Rn , n2 , be two bounded open sets with C boundaries such that Ω2\Ω1¯ . Suppose that Ω1Ω2=U¯ where UΩ1 is open with C boundary. Let G:C(Ω2¯)C(Ω2¯) , aw , be the solution operator to the Dirichlet problem,

-Δw=ainΩ2,w|Ω2=0.

Graph

The following result is an extension of [[8], Lemma 2.2] and [[24], Lemma 2.2], where the similar density results were obtained in the L2 and H1 topologies, respectively.

Lemma 2.5

The space

W:={Ga|Ω1:aC(Ω2¯),supp(a)Ω2\Ω1¯}

Graph

is dense in the space

S:={uC(Ω1¯):-Δu=0inΩ1,u|Ω1Ω2=0},

Graph

with respect to the W1,p(Ω1) -topology, for any 1<p< .

Proof

We shall follow the proof of [[24], Lemma 2.2] closely, adapting it to the Lp based Sobolev spaces. Let vW~-1,p(Ω1) , 1p+1p=1 , be such that

2.20 (v,Ga|Ω1)W~-1,p(Ω1),W1,p(Ω1)=0

Graph

for any aC(Ω2¯) , supp(a)Ω2\Ω1¯ . In view of the Hahn–Banach theorem, we have to prove that

(v,u)W~-1,p(Ω1),W1,p(Ω1)=0,

Graph

for any uS .

To that end, we first note that as GaC(Ω2¯) and Ga|Ω2=0 , we have GaW01,p(Ω2) . By [[2], Proposition 9.18], we can view Ga as an element of W1,p(Rn) via an extension by 0 to Rn\Ω2 . By the definition of W01,p(Ω2) , there exists a sequence φjC0(Ω2) such that φjGa in W1,p(Rn) . We have in view of (2.20) that

2.21 0=(v,Ga)W-1,p(Rn),W1,p(Rn)=limj(v,φj)W-1,p(Rn),W1,p(Rn)=limj(v,φj)W-1,p(Ω2),W01,p(Ω2)=(v,Ga)W-1,p(Ω2),W01,p(Ω2).

Graph

Next, Proposition 2.2 implies that there is a sequence vjC0(Ω1) such that vjv in W-1,p(Rn) . Consider the following Dirichlet problems,

2.22 -Δf=v|Ω2W-1,p(Ω2)inΩ2,f=0onΩ2,-Δfj=vjinΩ2,fj=0onΩ2.

Graph

By Theorem 2.3, the problems (2.22) have unique solutions fW01,p(Ω2) and fjC(Ω2¯)W01,p(Ω2) , respectively.

Using (2.21), (2.22), we get

2.23 0=(v,Ga)W-1,p(Ω2),W01,p(Ω2)=limj(vj,Ga)W-1,p(Ω2),W01,p(Ω2)=limj(-Δfj,Ga)W-1,p(Ω2),W01,p(Ω2)=limjΩ2(-Δfj)Ga¯dx=limjΩ2fja¯dx=Ω2fa¯dx.

Graph

Here we have used Green's formula, the fact that fj|Ω2=Ga|Ω2=0 , and that

f-fjW1,p(Ω2)Cv-vjW-1,p(Rn),

Graph

which is a consequence of Theorem 2.3.

It follows from (2.23) that f=0 in Ω2\Ω1¯ . This together with the fact that fW01,p(Ω2) , in view of [[2], Proposition 9.18], allows us to conclude that fW01,p(Ω1) . Thus, there exists a sequence f^jC0(Ω1) be such that f^jf in W1,p(Rn) , and therefore, -Δf^j-Δf in W-1,p(Rn) .

Let uS and let Ext(u)W1,p(Rn) be an extension of u. Using Green's formula, we get

2.24 (-Δf,Ext(u))W-1,p(Rn),W1,p(Rn)=limj((-Δf^j),Ext(u))W-1,p(Rn),W1,p(Rn)=limjΩ1(-Δf^j)u¯dx=0.

Graph

Let g=-Δf-vW-1,p(Rn) . We have that supp(g)Ω1 , in view of the fact that supp(v),supp(f)Ω1¯ , and (2.22). An application of Proposition 2.4 gives therefore

g=hδΩ1,hBp,p-(1-1/p)(Ω1).

Graph

It also follows from (2.22) that supp(g)Ω1Ω2=U¯ , and hence, supp(h)U¯ . Here UΩ1 is a bounded open set with C boundary, and therefore, there exists a sequence hjC0(U) such that hjh in Bp,p-(1-1/p)(Ω1) , see [[38], Section 4.3.2, p. 318]. Thus, we get

2.25 (g,Ext(u))W-1,p(Rn),W1,p(Rn)=(h,u|Ω1)Bp,p-(1-1/p)(Ω1),W1-1/p,p(Ω1)=limj(hj,u|Ω1)Bp,p-(1-1/p)(Ω1),Bp,p1-1/p(Ω1)=limjΩ1hju¯dS=0,

Graph

where the last equality follows from the fact that u|Ω1Ω2=0 . Combining (2.24) and (2.25), we see that

(v,u)W~-1,p(Ω1),W1,p(Ω1)=(-Δf,Ext(u))W-1,p(Rn),W1,p(Rn)-(g,Ext(u))W-1,p(Rn),W1,p(Rn)=0.

Graph

From local to global results. Completion of proof of Theorem 1.8

We follow [[8]]. Let Γ~=Ω\Γ . Assuming that f satisfies (1.3) and using Proposition 2.1, we would like to show that f vanishes inside Ω . To that end, let x0Γ and let us fix a point x1Ω . Let θ:[0,1]Ω¯ be a C1 curve joining x0 to x1 such that θ(0)=x0 , θ(0) is the interior normal to Ω at x0 and θ(t)Ω , for all t(0,1] . We set

Θε(t)={xΩ¯:d(x,θ([0,t]))ε}

Graph

and

I={t[0,1]:fvanishes a.e. onΘε(t)Ω}.

Graph

By Proposition 2.1, we have 0I if ε>0 is small enough. First as in [[8]], I is a closed subset of [0, 1]. If one proves that I is open then I=[0,1] due to the fact that [0, 1] is connected. This implies that x1supp(f) , and as x1 is an arbitrary point of Ω , we conclude that f=0 in Ω , and this will complete the proof of Theorem 1.8. Hence, we only need to prove that the set I is open in [0, 1].

To this end, let tI and ε>0 be small enough so that Θε(t)ΩΓ . Arguing as in [[8], [24]], we smooth out Ω\Θε(t) into an open subset Ω1 of Ω with smooth boundary such that

Ω1Ω\Θε(t),ΩΩ1Γ~,

Graph

and Ω1Ω=U¯ where UΩ1 is an open set with C boundary. By smoothing out the set ΩB(x0,ε) , with 0<εε sufficiently small, we enlarge the set Ω into an open set Ω2 with smooth boundary so that

Ω2ΩΩ1Ω=Ω1Ω2Γ~.

Graph

Let GΩ2 be the Green kernel associated to the open set Ω2 ,

-ΔyGΩ2(x,y)=δ(x-y),GΩ2(x,·)|Ω2=0.

Graph

We have GΩ2(x,y)C(Ω×Ω¯\{x=y}) , see [[40], Section 8.1]. Let us consider

v(x(1),,x(m+1))=Ω1f(y)(k=1mr=1,rkm(ω·yGΩ2(x(r),y))yGΩ2(x(k),y))·yGΩ2(x(m+1),y)dy,

Graph

where x(1),,x(m+1)Ω2\Ω1¯ . The function v is harmonic in all variables x(1),,x(m+1)Ω2\Ω1¯ . Since f=0 on Θε(t)Ω , we have

v(x(1),,x(m+1))=Ωf(y)(k=1mr=1,rkm(ω·yGΩ2(x(r),y))yGΩ2(x(k),y))·yGΩ2(x(m+1),y)dy,

Graph

where x(1),,x(m+1)Ω2\Ω1¯ . Now when x(l)Ω2\Ω¯ , the Green function GΩ2(x(l),·)C(Ω¯) is harmonic on Ω , and GΩ2(x(l),·)|Γ~=0 . By the orthogonality condition (1.3), we have v(x(1),,x(m+1))=0 when x(l)Ω2\Ω¯ , l=1,,m+1 .

As v(x(1),,x(m+1)) is harmonic in all variables x(1),,x(m+1)Ω2\Ω1¯ , and Ω2\Ω1¯ is connected, by unique continuation, we get that v(x(1),,x(m+1))=0 when x(1),,x(m+1)Ω2\Ω1¯ , i.e.

2.26 Ω1f(y)(k=1mr=1,rkm(ω·yGΩ2(x(r),y))yGΩ2(x(k),y))·yGΩ2(x(m+1),y)dy=0,x(1),,x(m+1)Ω2\Ω1¯.

Graph

Let alC(Ω2¯) , supp(al)Ω2\Ω1¯ , l=1,,m+1 . Multiplying (2.26) by a1(x(1))am+1(x(m+1)) , and integrating, we get

2.27 Ω1f(y)(k=1mr=1,rkmΩ2(ω·yGΩ2(x(r),y))ar(x(r))dx(r)Ω2yGΩ2(x(k),y)ak(x(k))dx(k))·Ω2yGΩ2(x(m+1),y)am+1(x(m+1))dx(m+1)dy=0.

Graph

Now it follows from the definition of W in Lemma 2.5 that any vW is given by

v(y)=Ω2GΩ2(x,y)a(x)dx,yΩ1,

Graph

where aC(Ω2¯) , supp(a)Ω2\Ω1¯ . This together with (2.27) gives that

2.28 Ω1f(y)(k=1mr=1,rkm(ω·v(r))v(k))·v(m+1)dy=0,

Graph

for all v(1),,v(m+1)W .

The (m+1) -linear form,

W1,m+1(Ω1)××W1,m+1(Ω1)C,(v(1),,v(m))Ω1f(y)(k=1mr=1,rkm(ω·v(r))v(k))·v(m+1)dy

Graph

is continuous in view of Hölder's inequality. An application of Lemma 2.5 with p=m+1 shows that (2.28) holds for all v(1),,v(m)C(Ω1¯) harmonic in Ω1 which vanish on Ω1Ω2 . Proposition 2.1 implies that f vanishes on a neighborhood of Ω1\(Ω1Ω2) , and therefore, I is an open set. The proof of Theorem 1.8 is complete.

Proof of Theorem 1.2

First it follows from (i) and (ii) that for each τC fixed, γ can be expanded into a power series

3.1 γ(x,τ,z)=1+k=1zkγ(x,τ,0)zkk!,zkγ(x,τ,0)C1,α(Ω¯),τ,zC,

Graph

converging in the Cα(Ω¯) topology. Furthermore, the map Cτzkγ(x,τ,0) is holomorphic with values in Cα(Ω¯) .

Let λΣ be arbitrary but fixed. Let ε=(ε1,,εm)Cm , m2 , and consider the Dirichlet problem (1.2) with

3.2 f=k=1mεkfk,fkC(Ω),supp(fk)Γ,k=1,,m.

Graph

Then for all |ε| sufficiently small, the problem (1.2) has a unique solution u(·;ε)C2,α(Ω¯) close to λ in C2,α(Ω¯) -topology, which depends holomorphically on εneigh(0,Cm) , with values in C2,α(Ω¯) .

We shall use an induction argument on m2 to prove that the equality

Λγ1Γ(λ+k=1mεkfk)=Λγ2Γ(λ+k=1mεkfk),

Graph

for all |ε| sufficiently small and all fkC(Ω) , supp(fk)Γ , k=1,,m , gives that zm-1γ1(x,λ,0)=zm-1γ1(x,λ,0) .

First let m=2 and we proceed to carry out a second order linearization of the partial Dirichlet-to-Neumann map. Let uj=uj(x;ε)C2,α(Ω¯) be the unique solution close to λ in C2,α(Ω¯) -topology of the Dirichlet problem,

3.3 Δuj+div(k=1zkγj(x,uj,0)(ω·uj)kk!uj)=0inΩ,uj=λ+ε1f1+ε2f2onΩ,

Graph

for j=1,2 . The solution uj is C with respect to ε for |ε| sufficiently small in view of Theorem B.1. Applying εl|ε=0 , l=1,2 , to (3.3), and using that uj(x,0)=λ , we get

3.4 Δvj(l)=0inΩ,vj(l)=flonΩ,

Graph

where vj(l)=εluj|ε=0 . It follows that v(l):=v1(l)=v2(l)C(Ω¯) .

Applying ε1ε2|ε=0 to (3.3) and letting wj=ε1ε2uj|ε=0 , we obtain that

3.5 Δwj+div(zγj(x,λ,0)((ω·v(1))v(2)+(ω·v(2))v(1)))=0inΩ,wj=0onΩ,

Graph

j=1,2 .

The fact that Λγ1Γ(λ+ε1f1+ε2f2)=Λγ1Γ(λ+ε1f1+ε2f2) for all small ε , and all f1,f2C(Ω) with supp(f1),supp(f2)Γ , gives that

3.6

Graph

An application of ε1ε2|ε=0 to (3.6) yields that

3.7 (νw1-νw2)|Γ+(zγ1(x,λ,0)-zγ2(x,λ,0))×((ω·v(1))νv(2)+(ω·v(2))νv(1))|Γ=0.

Graph

Multiplying the difference of two equations in (3.5) by v(3)C(Ω¯) harmonic in Ω , integrating over Ω , using Green's formula and (3.7), we obtain that

3.8 Ω(zγ1(x,λ,0)-zγ2(x,λ,0))((ω·v(1))v(2)+(ω·v(2))v(1))·v(3)dx=Ω\Γ(zγ1(x,λ,0)-zγ2(x,λ,0))((ω·v(1))νv(2)+(ω·v(2))νv(1))v(3)dS+Ω\Γ(νw1-νw2)v(3)dS=0,

Graph

provided that supp(v(3)|Ω)Γ . It follows from (3.8) that

3.9 Ω(zγ1(x,λ,0)-zγ2(x,λ,0))((ω·v(1))v(2)+(ω·v(2))v(1))·v(3)dx=0,

Graph

for all v(l)C(Ω¯) harmonic in Ω such that supp(v(l)|Ω)Γ , l=1,2,3 . An application of Theorem 1.8 with m=2 allows us to conclude that zγ1(·,λ,0)=zγ2(·,λ,0) in Ω . Now as λΣ is arbitrary and the functions Cτzγj(x,τ,0) , j=1,2 , are holomorphic, by the uniqueness properties of holomorphic functions, we have zγ1(·,·,0)=zγ2(·,·,0) in Ω¯×C .

Let m3 and assume that

3.10 zkγ1(·,·,0)=zkγ2(·,·,0)inΩ¯×C,

Graph

for all k=1,,m-2 . Let λΣ be arbitrary but fixed. To prove that zm-1γ1(·,λ,0)=zm-1γ2(·,λ,0) in Ω¯ , we carry out the mth order linearization of the partial Dirichlet-to-Neumann map. In doing so, we let uj=uj(x;ε)C2,α(Ω¯) be the unique solution close to λ in C2,α(Ω¯) -topology of the Dirichlet problem,

3.11 Δuj+div(k=1zkγj(x,uj,0)(ω·uj)kk!uj)=0inΩ,uj=λ+ε1f1++εmfmonΩ,

Graph

for j=1,2 . We shall next apply ε1εm|ε=0 to (3.11). To this end, we first note that ε1εm(k=mzkγj(x,uj,0)(ω·uj)kk!uj) is a sum of terms each of them containing positive powers of uj , which vanishes when ε=0 . The only term in ε1εm(zm-1γj(x,uj,0)(ω·uj)m-1(m-1)!uj) which does not contain a positive power of uj is

3.12 zm-1γj(x,uj,0)(k=1mr=1,rkm(ω·εruj)εkuj).

Graph

Finally, the expression ε1εm(k=1m-2zkγj(x,uj,0)(ω·uj)kk!uj)|ε=0 is independent of j=1,2 . Indeed, this follows from (3.10), the fact that this expression contains only the derivatives of uj of the form εl1,,εlssuj|ε=0 with s=1,,m-1 , εl1,,εls{ε1,,εm} , and the fact that

3.13 εl1,,εlssu1|ε=0=εl1,,εlssu2|ε=0,

Graph

for s=1,,m-1 , εl1,,εls{ε1,,εm} . The latter can be seen by induction on s, applying the operator εl1,,εlss|ε=0 to (3.11) and using (3.10) as well as the unique solvability of the Dirichlet problem for the Laplacian. Thus, an application ε1εm|ε=0 to (3.11) gives

3.14 Δwj+div(zm-1γj(x,λ,0)(k=1mr=1,rkm(ω·v(r))v(k)))=HminΩ,wj=0onΩ,

Graph

cf. (3.12). Here wj=ε1εmuj|ε=0 and

Hm(x,λ):=-div(ε1εm(k=1m-2zkγj(x,uj,0)(ω·uj)kk!uj)|ε=0).

Graph

The fact that Λγ1Γ(λ+ε1f1++εmfm)=Λγ1Γ(λ+ε1f1++εmfm) for all small ε and all fkC(Ω) with supp(fk),Γ , k=1,,m , yields (3.6). Applying of ε1εm|ε=0 to (3.6), using (3.10) and (3.13), we obtain that

3.15 (νw1-νw2)|Γ+(zm-1γ1(x,λ,0)-zm-1γ2(x,λ,0))(k=1mr=1,rkm(ω·v(r))νv(k))|Γ=0.

Graph

Using (3.14), (3.15), and proceeding as in the case m=2 , we get

3.16 Ω(zm-1γ1(x,λ,0)-zm-1γ1(x,λ,0))(k=1mr=1,rkm(ω·v(r))v(k))·v(m+1)dx=0,

Graph

for all v(l)C(Ω¯) harmonic in Ω such that supp(v(l)|Ω)Γ , l=1,,m+1 . Applying Theorem 1.8, we conclude that zm-1γ1(·,λ,0)=zm-1γ2(·,λ,0) in Ω¯ . Now as λΣ is arbitrary and the functions Cτzm-1γj(x,τ,0) , j=1,2 , are holomorphic, we have zm-1γ1(·,·,0)=zm-1γ2(·,·,0) in Ω¯×C . This completes the proof of Theorem 1.2.

Proof of Theorem 1.1

First it follows from (a) and (b) that γ can be expanded into the following power series,

4.1 γ(x,λ)=1+k=1λkγ(x,0)λkk!,λkγ(x,0)C1,α(Ω¯),λC,

Graph

converging in the C1,α(Ω¯) topology.

Let ε=(ε1,,εm)Cm , m2 , and consider the Dirichlet problem (1.1) with f given by (3.2). For all |ε| sufficiently small, the problem (1.1) has a unique small solution u(·;ε)C2,α(Ω¯) , which depends holomorphically on εneigh(0,Cm) .

As in the proof of Theorem 1.2, we use an induction argument on m2 to show that Λγ1Γ=Λγ2Γ implies that λm-1γ1(x,0)=λm-1γ1(x,0) .

First let m=2 and we perform a second order linearization of the partial Dirichlet-to-Neumann map. Let uj=uj(x;ε)C2,α(Ω¯) be the unique solution small solution of the Dirichlet problem,

4.2 Δuj+div(k=1λkγj(x,0)ujkk!uj)=0inΩ,uj=ε1f1+ε2f2onΩ,

Graph

for j=1,2 . Applying εl|ε=0 , l=1,2 , to (4.2), and using that uj(x,0)=0 , we see that

4.3 Δvj(l)=0inΩ,vj(l)=flonΩ,

Graph

where vj(l)=εluj|ε=0 . We have therefore v(l):=v1(l)=v2(l)C(Ω¯) .

Applying ε1ε2|ε=0 to (4.2) and setting wj=ε1ε2uj|ε=0 , we get

4.4 Δwj+div(λγj(x,0)(v(1)v(2)+v(2)v(1)))=0inΩ,wj=0onΩ,

Graph

j=1,2 . The fact that Λγ1Γ(ε1f1+ε2f2)=Λγ1Γ(ε1f1+ε2f2) for all small ε , and all f1,f2C(Ω) with supp(f1),supp(f2)Γ , implies that

4.5 (1+k=1λkγ1(x,0)u1kk!)νu1|Γ=(1+k=1zkγ2(x,0)u2kk!)νu2|Γ.

Graph

Applying ε1ε2|ε=0 to (4.5), we get

4.6 (νw1-νw2)|Γ+(λγ1(x,0)-λγ2(x,0))(v(1)νv(2)+v(2)νv(1))|Γ=0.

Graph

Multiplying the difference of two equations in (4.4) by v(3)C(Ω¯) harmonic in Ω , integrating over Ω , using Green's formula and (4.6), we obtain that

4.7 Ω(λγ1(x,0)-λγ2(x,0))(v(1)v(2)+v(2)v(1))·v(3)dx=Ω\Γ(λγ1(x,0)-λγ2(x,0))(v(1)νv(2)+v(2)νv(1))v(3)dS+Ω\Γ(νw1-νw2)v(3)dS=0,

Graph

provided that supp(v(3)|Ω)Γ . Thus, (4.7) gives that

Ω(λγ1(x,0)-λγ2(x,0))(v(1)v(2)+v(2)v(1))·v(3)dx=0,

Graph

for all v(l)C(Ω¯) harmonic in Ω such that supp(v(l)|Ω)Γ , l=1,2,3 . By Theorem 1.9 with m=2 , we get λγ1(·,0)=λγ2(·,0) in Ω¯ .

Let m3 and assume that λkγ1(·,0)=λkγ2(·,0)inΩ¯ , for all k=1,,m-2 . To prove that λm-1γ1(·,0)=λm-1γ2(·,·,0) in Ω¯ , we perform the mth order linearization of the partial Dirichlet-to-Neumann map. In doing so, we let uj=uj(x;ε)C2,α(Ω¯) be the unique small solution of the Dirichlet problem,

4.8 Δuj+div(k=1λkγj(x,0)ujkk!uj)=0inΩ,uj=ε1f1++εmfmonΩ,

Graph

for j=1,2 . Applying ε1εm|ε=0 to (4.8), and arguing as in Theorem 1.2, we obtain that

4.9 Δwj+div(λm-1γj(x,0)(k=1mr=1,rkmv(r)v(k)))=HminΩ,wj=0onΩ.

Graph

Here wj=ε1εmuj|ε=0 and

Hm(x):=-div(ε1εm(k=1m-2λkγj(x,0)ujkk!uj)|ε=0),

Graph

which is independent of j.

Now the equality Λγ1Γ(ε1f1++εmfm)=Λγ1Γ(ε1f1++εmfm) for all small ε and all fkC(Ω) with supp(fk),Γ , k=1,,m , implies (4.5). Applying of ε1εm|ε=0 to (4.5), we obtain that

4.10 (νw1-νw2)|Γ+(λm-1γ1(x,0)-λm-1γ2(x,0))(k=1mr=1,rkmv(r)νv(k))|Γ=0.

Graph

Proceeding as in the case m=2 , and using (4.9), (4.10), we get

Ω(λm-1γ1(x,0)-λm-1γ1(x,,0))(k=1mr=1,rkmv(r)v(k))·v(m+1)dx=0,

Graph

for all v(l)C(Ω¯) harmonic in Ω such that supp(v(l)|Ω)Γ , l=1,,m+1 . An application of Theorem 1.9 allows us to conclude that λm-1γ1(·,0)=λm-1γ2(·,0) in Ω¯ . This completes the proof of Theorem 1.1.

Acknowledgements

We are very grateful to the referees for helpful comments and suggestions. The work of Y.K is partially supported by the French National Research Agency ANR (project MultiOnde) grant ANR-17-CE40-0029. The research of K.K. is partially supported by the National Science Foundation (DMS 1815922, DMS 2109199). The research of G.U. is partially supported by NSF, a Walker Professorship at UW and a Si-Yuan Professorship at IAS, HKUST. Part of the work was supported by the NSF grant DMS-1440140 while K.K. and G.U. were in residence at MSRI in Berkeley, California, during Fall 2019 semester.

Data availability statement

Data sharing not applicable to this article as no datasets were generated or analysed during the current study.

Declarations

Conflict of interest

The authors have no conflicts of interest to declare that are relevant to this article.

Appendix A. Proof of Theorem 1.2 in the case of full data

Note that the result of Theorem 1.2 is new even in the case of full data, i.e. Γ=Ω , and the purpose of this appendix is to present an alternative simple proof in this case.

Using the linearization of the Dirichlet-to-Neumann map ΛγΩ , we shall see below that the proof of Theorem 1.2 in the full data case will be a consequence of the following density result.

Proposition A.1

Let ΩRn , n2 , be a bounded open set with C boundary, let ωSn-1 be fixed and let m=2,3,, be fixed. Let fL(Ω) be such that

A.1 Ωf(ω·v1)m-1v1·v2dx=0,

Graph

for all functions v1,v2C(Ω¯) harmonic in Ω . Then f=0 in Ω .

Proof

Let ξSn-1 and consider kSn-1 such that ξ·k=0 . Let h>0 . Setting

v1(x)=e1hx·(k+iξ),v2(x)=emhx·(-k+iξ),

Graph

so that v1,v2C(Rn) and harmonic. Substituting v1 and v2 into (A.1) and using that (k+iξ)·(-k+iξ)=-2 , we get

(ω·(k+iξ))m-1Ωf(x)e2mhix·ξdx=0,

Graph

and therefore, we have

Ωf(x)e2mhix·ξdx=0,

Graph

for all ξSn-1 , ξ·ω0 , and all h>0 . Hence, f=0 .

Let λΣ be arbitrary but fixed. We shall use an induction argument on m2 to prove that the equality

Λγ1Γ(λ+k=1mεkfk)=Λγ2Γ(λ+k=1mεkfk),

Graph

for all |ε| sufficiently small and all fkC(Ω) , supp(fk)Γ , k=1,,m , gives that zm-1γ1(x,λ,0)=zm-1γ1(x,λ,0) .

First when m=2 , taking v(1)=v(2) in (3.9) and using Proposition A.1 with m=2 , we get zγ1(·,λ,0)=zγ2(·,λ,0) in Ω . Now as λΣ is arbitrary, we have zγ1(·,·,0)=zγ2(·,·,0) in Ω¯×C .

Let m=3,4, . Let λΣ be arbitrary but fixed. Letting v(1)==v(m) in (3.16) and using Proposition A.1, we see that zm-1γ1(·,λ,0)=zm-1γ2(·,λ,0) in Ω . Again, as λΣ is arbitrary, we get zm-1γ1(·,·,0)=zm-1γ2(·,·,0) in Ω¯×C . This completes the proof of Theorem 1.2 in the full data case.

Appendix B. Well-posedness of the Dirichlet problem for a quasilinear conductivity equation

In this appendix we shall recall a standard argument for showing the well-posedness of the Dirichlet problem for a quasilinear conductivity equation.

Let ΩRn , n2 , be a bounded open set with C boundary. Let kN{0} and 0<α<1 and let Ck,α(Ω¯) be the standard Hölder space on Ω , see [[16], [24]]. We observe that Ck,α(Ω¯) is an algebra under pointwise multiplication, with

B.1 uvCk,α(Ω¯)C(uCk,α(Ω¯)vL(Ω)+uL(Ω)vCk,α(Ω¯)),u,vCk,α(Ω¯),

Graph

see [[16], Theorem A.7]. We write Cα(Ω¯)=C0,α(Ω¯) .

Let ωSn-1={ωRn,|ω|=1} , be fixed. Consider the Dirichlet problem for the following isotropic quasilinear conductivity equation,

B.2 div(γ(x,u,ω·u)u)=0inΩ,u=λ+fonΩ,

Graph

with λC . We assume that the function γ:Ω¯×C×CC satisfies the following conditions,

  • the map C×C(τ,z)γ(·,τ,z) is holomorphic with values in C1,α(Ω¯) with some 0<α<1 ,
  • γ(x,0,0)=1 .

It follows from (i) and (ii) that γ can be expand into a power series

B.3 γ(x,τ,z)=1+j+k1,j0,k0τjzkγ(x,0,0)τjzkj!k!,τjzkγ(x,0,0)C1,α(Ω¯),

Graph

converging in the C1,α(Ω¯) topology.

We have the following result.

Theorem B.1

Let λC be fixed. Then under the above assumptions, there exist δ>0 , C>0 such that for any fBδ(Ω):={fC2,α(Ω):fC2,α(Ω)<δ} , the problem (B.2) has a solution u=uλ,fC2,α(Ω¯) which satisfies

u-λC2,α(Ω¯)CfC2,α(Ω).

Graph

The solution u is unique within the class {uC2,α(Ω¯):u-λC2,α(Ω¯)<Cδ} and it is depends holomorphically on fBδ(Ω) . Furthermore, the map

Bδ(Ω)C1,α(Ω¯),fνu|Ω

Graph

is holomorphic.

Proof

Let λC be fixed, and let

B1=C2,α(Ω),B2=C2,α(Ω¯),B3=Cα(Ω¯)×C2,α(Ω).

Graph

Consider the map,

B.4 F:B1×B2B3,F(f,u)=(div(γ(x,u,ω·u)u),u|Ω-λ-f).

Graph

Following [[29]], we shall make use of the implicit function theorem for holomorphic maps between complex Banach spaces, see [[37], p. 144]. First we check that F enjoys the mapping property (B.4). To that end in view of the fact that C1,α(Ω¯) is an algebra under pointwise multiplication, we only need to show that γ(x,u,ω·u)C1,α(Ω¯) . In doing so, by Cauchy's estimates, we get

B.5 τjzkγ(x,0,0)C1,α(Ω¯)j!k!R1jR2ksup|τ|=R1,|z|=R2γ(·,τ,z)C1,α(Ω¯),R1,R2>0,

Graph

for all j0 , k0 , and j+k1 . With the help of (B.1) and (B.5), we obtain that

B.6 τjzkγ(x,0,0)uj(ω·u)kj!k!C1,α(Ω¯)Cj+kR1jR2kuC1,α(Ω¯)jω·uC1,α(Ω¯)ksup|τ|=R1,|z|=R2γ(·,τ,z)C1,α(Ω¯).

Graph

Taking R1=2CuC1,α(Ω¯) and R2=2Cω·uC1,α(Ω¯) , we see that the series

j+k1,j0,k0τjzkγ(x,0,0)uj(ω·u)kj!k!

Graph

converges in C1,α(Ω¯) . Hence, in view of (B.3), γ(x,u,ω·u)C1,α(Ω¯) .

Let us show that F in (B.4) is holomorphic. First F is locally bounded as it is continuous in (fu). Hence, we only need to check that F is weak holomorphic, see [[37], p. 133]. To that end, letting (f0,u0),(f1,u1)B1×B2 , we show that the map

Graph

is holomorphic in C with values in B3 . Clearly, we only have to check that the map μγ(x,u0(x)+μu1(x),ω·(u0(x)+μu1(x))) is holomorphic in C with values in C1,α(Ω¯) . This is a consequence of the fact that the series

j+k1,j0,k0τjzkγ(x,0,0)(u0+μu1)j(ω·(u0+μu1))kj!k!

Graph

converges in C1,α(Ω¯) , locally uniformly in λC , in view of (B.6).

We have F(0,λ)=0 and the partial differential uF(0,λ):B2B3 is given by

uF(0,λ)v=(Δv,v|Ω).

Graph

It follows from [[11], Theorem 6.15] that the map uF(0,λ):B2B3 is a linear isomorphism.

An application of the implicit function theorem, see [[37], p. 144], shows that there exists δ>0 and a unique holomorphic map S:Bδ(Ω)C2,α(Ω¯) such that S(0)=λ and F(f,S(f))=0 for all fBδ(Ω) . Letting u=S(f) and using that S is Lipschitz continuous and S(0)=λ , we have

u-λC2,α(Ω¯)CfC2,α(Ω).

Graph

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

References 1 Agranovich, M: Sobolev Spaces, their Generalizations and Elliptic Problems in Smooth and Lipschitz Domains, Springer Monographs in Mathematics. Springer, Cham (2015) 2 Brezis H. Functional Analysis, Sobolev Spaces and Partial Differential Equations. Universitext. 2011: New York; Springer. 10.1007/978-0-387-70914-7. 1220.46002 3 Browder, F: Functional analysis and partial differential equations. II. Math. Ann. 145, 81–226 (1961/62) 4 Cârstea, C, Feizmohammadi, A: A density property for tensor products of gradients of harmonic functions and applications, preprint, arXiv:2009.11217 5 Cârstea C, Feizmohammadi A. An inverse boundary value problem for certain anisotropic quasilinear elliptic equations. J. Differ. Equ. 2021; 284: 318-349. 4227095. 10.1016/j.jde.2021.02.044. 1484.35405 6 Cârstea, C, Feizmohammadi, A, Kian, Y, Krupchyk, K, Uhlmann, G: The Calderón inverse problem for isotropic quasilinear conductivities. Adv. Math. 391, Paper No. 107956 (2021) 7 Cârstea C, Nakamura G, Vashisth M. Reconstruction for the coefficients of a quasilinear elliptic partial differential equation. Appl. Math. Lett. 2019; 98: 121-127. 3964222. 10.1016/j.aml.2019.06.009. 1423.35437 8 Dos Santos Ferreira D, Kenig C, Sjöstrand J, Uhlmann G. On the linearized local Calderón problem. Math. Res. Lett. 2009; 16; 6: 955-970. 2576684. 10.4310/MRL.2009.v16.n6.a4. 1198.31003 9 Eskin G. Lectures on Linear Partial Differential Equations, Graduate Studies in Mathematics. 2011: Providence; American Mathematical Society Feizmohammadi A, Oksanen L. An inverse problem for a semi-linear elliptic equation in Riemannian geometries. J. Differ. Equ. 2020; 269; 6: 4683-4719. 4104456. 10.1016/j.jde.2020.03.037. 1448.58016 Gilbarg D, Trudinger N. Elliptic Partial Differential Equations of Second Order, Reprint of the Classics in Mathematics. 20011998: Berlin; Springer. 10.1007/978-3-642-61798-0. 1042.35002 Glowinski R, Rappaz J. Approximation of a nonlinear elliptic problem arising in a non-Newtonian fluid flow model in glaciology. ESAIM: Math. Modell. Numer. Anal. 2003; 37: 175-186. 1972657. 10.1051/m2an:2003012. 1046.76002 Grisvard P. Elliptic Problems in Nonsmooth Domains, Monographs and Studies in Mathematics, 24. 1985: Boston; Pitman (Advanced Publishing Program). 0695.35060 Grubb G. Distributions and Operators, Graduate Texts in Mathematics, 252. 2009: New York; Springer Hervas D, Sun Z. An inverse boundary value problem for quasilinear elliptic equations. Commun. Partial Differ. Equ. 2002; 27; 11–12: 2449-2490. 1944036. 10.1081/PDE-120016164. 1013.35084 Hörmander L. The boundary problems of physical geodesy. Arch. Rational Mech. Anal. 1976; 62; 1: 1-52. 602181. 10.1007/BF00251855. 0331.35020 Imanuvilov O, Uhlmann G, Yamamoto M. The Calderón problem with partial data in two dimensions. J. Am. Math. Soc. 2010; 23; 3: 655-691. 10.1090/S0894-0347-10-00656-9. 1201.35183 Imanuvilov O, Yamamoto M. Unique determination of potentials and semilinear terms of semilinear elliptic equations from partial Cauchy data. J. Inverse Ill-Posed Probl. 2013; 21; 1: 85-108. 3023500. 10.1515/jip-2012-0033. 1273.35318 Isakov V. Uniqueness of recovery of some quasilinear partial differential equations. Commun. Partial Differ. Equ. 2001; 26; 11–12: 1947-1973. 1876409. 10.1081/PDE-100107813. 0997.35018 Isakov V, Nachman A. Global uniqueness for a two-dimensional semilinear elliptic inverse problem. Trans. Am. Math. Soc. 1995; 347; 9: 3375-3390. 1311909. 10.1090/S0002-9947-1995-1311909-1. 0849.35148 Isakov V, Sylvester J. Global uniqueness for a semilinear elliptic inverse problem. Commun. Pure Appl. Math. 1994; 47; 10: 1403-1410. 1295934. 10.1002/cpa.3160471005. 0817.35126 Kang K, Nakamura G. Identification of nonlinearity in a conductivity equation via the Dirichlet-to-Neumann map. Inverse Problems. 2002; 18; 4: 1079-1088. 1929283. 10.1088/0266-5611/18/4/309. 1043.35146 Kapanadze D, Mishuris G, Pesetskaya E. Exact solution of a nonlinear heat conduction problem in a doubly periodic 2D composite material. Arch. Mech. 2015; 67; 2: 157-178. 3381129. 1329.80005arXiv:1403.7914 Krupchyk K, Uhlmann G. Partial data inverse problems for semilinear elliptic equations with gradient nonlinearities. Math. Res. Lett. 2020; 27: 1801-1824. 4216606. 10.4310/MRL.2020.v27.n6.a10. 1459.35399 Krupchyk K, Uhlmann G. A remark on partial data inverse problems for semilinear elliptic equations. Proc. Am. Math. Soc. 2020; 148; 2: 681-685. 4052205. 10.1090/proc/14844. 1431.35246 Krupchyk, K, Uhlmann, G: Inverse problems for nonlinear magnetic Schrödinger equations on conformally transversally anisotropic manifolds, preprint, arXiv:2009.05089 (to appear in Analysis and PDE) Kurylev Y, Lassas M, Uhlmann G. Inverse problems for Lorentzian manifolds and non-linear hyperbolic equations. Invent. Math. 2018; 212; 3: 781-857. 3802298. 10.1007/s00222-017-0780-y. 1396.35074 Lai, R.-Y, Zhou, T: Partial Data Inverse Problems for Nonlinear Magnetic Schrödinger Equations, preprint arXiv:2007.02475 (to appear in Mathematical Research Letters) Lassas M, Liimatainen T, Lin Y-H, Salo M. Inverse problems for elliptic equations with power type nonlinearities. J. Math. Pures Appl. 2021; 145: 44-82. 4188325. 10.1016/j.matpur.2020.11.006. 1460.35395 Lassas M, Liimatainen T, Lin Y-H, Salo M. Partial data inverse problems and simultaneous recovery of boundary and coefficients for semilinear elliptic equations. Revista Matemática Iberoamericana. 2021; 37: 1553-1580. 4269409. 10.4171/rmi/1242. 1473.35653 McLean W. Strongly Elliptic Systems and Boundary Integral Equations. 2000: Cambridge; Cambridge University Press. 0948.35001 Medková D. The Laplace Equation, Boundary Value Problems on Bounded and Unbounded Lipschitz Domains. 2018: Cham; Springer. 1457.35002 Mironescu, P: Note on Gagliardo's theorem "". Ann. Univ. Buchar. Math. Ser. 6(LXIV)(1), 99–103 (2015) Muñoz, C, Uhlmann, G: The Calderón problem for quasilinear elliptic equations. Ann. Inst. H. Poincaré Anal. Non Linéaire 37(5), 1143–1166 (2020) Nachman A. Global uniqueness for a two-dimensional inverse boundary value problem. Ann. Math. (2). 1996; 143; 1: 71-96. 1370758. 10.2307/2118653. 0857.35135 Ponte Castañeda P, Kailasam M. Nonlinear electrical conductivity in heterogeneous media. Proc. R. Soc. Lond. Ser. A. 1997; 453; 1959: 793-816. 1447153. 10.1098/rspa.1997.0044. 0886.73039 Pöschel J, Trubowitz E. Inverse Spectral Theory, Pure and Applied Mathematics, 130. 1987: Boston; Academic Press Inc. 0623.34001 Triebel H. Interpolation theory, function spaces, differential operators. 1978: Amsterdam; North-Holland. 0387.46033Mathematical Library, 18 Triebel H. Spaces of Besov-Hardy-Sobolev type on complete Riemannian manifolds. Ark. Mat. 1986; 24; 2: 299-337. 884191. 0664.46026 Schlag W. A Course in Complex Analysis and Riemann Surfaces, Graduate Studies in Mathematics, 154. 2014: Providence; American Mathematical Society. 10.1090/gsm/154 Shankar R. Recovering a quasilinear conductivity from boundary measurements. Inverse Problems. 2020; 27: 015014. 4197837. 1457.65172 Sun Z. On a quasilinear inverse boundary value problem. Math. Z. 1996; 221; 2: 293-305. 1376299. 10.1007/PL00022738. 0843.35137 Sun Z. Inverse boundary value problems for a class of semilinear elliptic equations. Adv. Appl. Math. 2004; 32; 4: 791-800. 2053845. 10.1016/j.aam.2003.06.001. 1059.35176 Sun Z. Conjectures in inverse boundary value problems for quasilinear elliptic equations. Cubo. 2005; 7; 3: 65-73. 2191048. 1100.35123 Sun Z. An inverse boundary-value problem for semilinear elliptic equations. Electron. J. Differ. Equ. 2010; 2010; 37: 5. 2602870. 1189.35380 Sun Z, Uhlmann G. Inverse problems in quasilinear anisotropic media. Am. J. Math. 1997; 119; 4: 771-797. 1465069. 10.1353/ajm.1997.0027. 0886.35176 Sylvester J, Uhlmann G. A global uniqueness theorem for an inverse boundary value problem. Ann. Math. (2). 1987; 125; 1: 153-169. 873380. 10.2307/1971291. 0625.35078

By Yavar Kian; Katya Krupchyk and Gunther Uhlmann

Reported by Author; Author; Author

Titel:
Partial data inverse problems for quasilinear conductivity equations.
Autor/in / Beteiligte Person: Kian, Yavar ; Krupchyk, Katya ; Uhlmann, Gunther
Link:
Zeitschrift: Mathematische Annalen, Jg. 385 (2023-04-01), Heft 3/4, S. 1611-1638
Veröffentlichung: 2023
Medientyp: academicJournal
ISSN: 0025-5831 (print)
DOI: 10.1007/s00208-022-02367-y
Schlagwort:
  • INVERSE problems
  • HARMONIC functions
  • SEMILINEAR elliptic equations
  • EQUATIONS
  • Subjects: INVERSE problems HARMONIC functions SEMILINEAR elliptic equations EQUATIONS
Sonstiges:
  • Nachgewiesen in: DACH Information
  • Sprachen: English
  • Document Type: Article
  • Author Affiliations: 1 = Aix Marseille University, Université de Toulon, CNRS, CPT, Marseille, France ; 2 = Department of Mathematics, University of California, 92697-3875, Irvine, CA, USA ; 3 = Department of Mathematics, University of Washington, 98195-4350, Seattle, WA, USA ; 4 = Institute for Advanced Study of the Hong Kong University of Science and Technology, Kowloon, Hong Kong
  • Full Text Word Count: 14306

Klicken Sie ein Format an und speichern Sie dann die Daten oder geben Sie eine Empfänger-Adresse ein und lassen Sie sich per Email zusenden.

oder
oder

Wählen Sie das für Sie passende Zitationsformat und kopieren Sie es dann in die Zwischenablage, lassen es sich per Mail zusenden oder speichern es als PDF-Datei.

oder
oder

Bitte prüfen Sie, ob die Zitation formal korrekt ist, bevor Sie sie in einer Arbeit verwenden. Benutzen Sie gegebenenfalls den "Exportieren"-Dialog, wenn Sie ein Literaturverwaltungsprogramm verwenden und die Zitat-Angaben selbst formatieren wollen.

xs 0 - 576
sm 576 - 768
md 768 - 992
lg 992 - 1200
xl 1200 - 1366
xxl 1366 -