## Computational & Applied Mathematics

*On-line version* ISSN 1807-0302

### Comput. Appl. Math. vol.25 no.2-3 Petrópolis 2006

**On an inverse boundary value problem ^{*}**

**Alberto P. Calderón**

Department of Mathematics, The University of Chicago

**ABSTRACT**

This paper is a reprint of the original work by A. P. Calderón published by the Brazilian Mathematical Society (SBM) in ATAS of SBM (Rio de Janeiro), pp. 65-73, 1980. The original paper had no abstract, so this reprint to be truthful to the original work is published with no abstract.

**Mathematical subject classification:** 35J25, 35Q60, 47F05, 86A20, 86A22.

**Key words:** inverse problems, boundary value problems, identification problem, elliptic equations.

In this note we shall discuss the following problem. Let *D* be a boundedDomain in ^{n}, *n* __>__ 2, with Lipschitzian boundary *dD*, and g be a real bounded measurable function in *D* with a positive lower bound. Consider the differential operator

*L*_{g} (*w*) = Ñ . (g Ñ *w*)

acting on function of *H*^{1}(*D*) and the quadratic form *Q*_{g}(f) , where the functions in *H*^{1}(* ^{n}*) , defined by

The problem is then to decide whether g is uniquely determined by *Q*_{g} and to calculate g in terms *Q*_{g}, if g is indeedDetermined by *Q*_{g}.

This problem originates in the following problem of electrical prospection. If *D* represents an inhomogeneous conducting body with electrical conductivity g, determine g by means of direct current steady state electrical measurements carried out on the surface of *D*, that is without penetrating *D*. In this physical situation *Q*_{g}(f) represents the power necessary to maintain an electrical potential g in ¶*D*.

In principle *Q*_{g} can be determined through measurements effected on *dD* and contains all information about g which can be thus obtained.

But let us return to our mathematical problem. Let us introduce the following norms in the space of functions g on *dD* and in the space of quadratic forms *Q*(f)

Then the mapping

F : g ®* Q*_{g}

is bounded and analytic in the subset of *L*^{¥}(*D*) consisting of functions which are real and have a positive lower bound. Our goal is then to determine whether F is injective, and invert F if this is the case. This we are yet unable to do and is, as far as we know, an open problem. However we shall show that *d*F|_{g = const} is indeed injective, that is, the linearized problem has an affirmative answer. If *d*F|_{g = const}, which is a linear operator, had a closer range, one could conclude that F itself is injective in a sufficiently small neighborhood of g = const. But the range of *d*F is not closed and the desired conclusion cannot be obtained in this fashion. Nevertheless, as we shall see below, if g is sufficiently close to a constant, it is nearly determined by *Q*_{g} and in some cases it can be calculated with an error much smaller than ||g - const||*L*^{¥}.

To show this let us first obtain an expression for the solution of the equation

*L*_{g} (*W*) = Ñ . (g Ñ *W*), *W*|*dD* = f. g = 1+d

Let *W* = *u* + *v*, where D*u* = *L*_{1} *u* = 0, *u*|_{dD} = f. Then

*L*_{g}(*W*) = *L*_{1+g} (*u*+*v*) = *L*_{1}*v* + *L*_{g}*v* + *L*_{g}*u* = 0

Since *u*|_{dD} = *W*|*dD* we have *v*|_{dD} = 0 and *v* Î (*D*), the closure in *H*^{1}(^{n}) of functions of *C*^{¥} with support in *D*. But *L*_{1}, as an operator from (D) into *H*^{-1}(D), has a bounded inverse *G*, and from the last expression we obtain

*v* + *GL*_{d}*v* = - *GL*_{d}*u*

and

Since for *W* Î (*D*), if A denotes the norm of *G*, the series above will converge for *A* < 1 and

From (1) it follows that f is analytic at g = 1. The same argument would show that f is analytic at any other point g.

Next let us calculate *d* f|_{g = 1}. We have

The contribution of the second term in the integrand of the last integral vanishes on account of the fact that D*u* = 0. Furthermore, from (1) one sees readily that the parts linear in d of the last two terms in the integrand vanish. Thus setting *d* g = d we obtain

To show that *d Q*_{g}(f) |_{g = 1} is injective, we merely have to show that if the last integral vanishes for all *u* with D*u* = 0 then d = 0. But if the integral vanishes for all such *u*, then we also have

whenever D*u*_{1} = D*u*_{2} = 0 in *D*. Now let *Z* be any vector in ^{n} and __ a__ another vector such that |

__| = |__

*a**Z*|,(

__.__

*a**Z*) = 0 . Then the functions

are harmonic, and substituting in (3) we obtain

whence it follows that d = 0.

Now let us return to *Q*_{g} (*W*) . We set again g = 1+d and introduce the bilinear form

and setting *W _{j}* =

*u*+

_{j}*v*,

_{j}*j*= 1,2, D

*u*= 0,

_{j}*u*|

_{j}_{dD}= f

*we obtain*

_{j} Now, substitution of the exponentials in (5) for *u*_{1} and *u*_{2} in the preceding expression (taking __ a__ to be a function of

*Z*such that |

*a*| = |

*Z*|, (

__.__

*a**Z*) = 0 ) yields

where (*Z*) is a Fourier transform of g extended to be zero outside *D*, the function

is known and, as follows readily from (2),

provided that *A*__ <__ 1-e, where *C* depends only on *D* and e, and *r* is the radius of the smallest sphere containing *D*. Now *R*(*Z*) is too large to permit estimating g(*x*). However, under favorable circumstances it is still possible to obtain satisfactory information about g. Choose a, 1 < a < 2, then for

we have |*R*(*Z*)| __<__ *C*. Let h be a function such that Î *C*^{¥}, supp Ì {|*Z*| __<__ 1} (0) = 1, and let h_{s} = s^{n}h(s*Z*). Then we have

and

where * denotes convolution and

where *C*_{1} depends only on *D*, a e e.

Thus if is sufficiently small, (9) gives an approximation for g*h_{s} with an error which is much smaller than . Clearly, if is small then s is large and g*h_{s} is itself in some sense, a good approximation to g.

Approximations to the function g itself may be obtained if one assumes that g, extended to be equal to 1 outside *D*, is in *C ^{m}*,

*m*>

*n*. In this case one obtains

(*Z*) = _{1} (*Z*) + (*Z*)

where *F*_{1} is known and *R*(*Z*) is the same as in (6). One then calculates d(*x*) by integrating over |*Z*| __<__ s with s as in (8) and estimates the error by using the decay of at ¥. Thus one obtains

g(*x*) = *F*_{2}(*x*) + s(*x*)

where *F*_{2}(*x*) is known and

where *M* is a bound for the derivatives of order *m* of g.

**REFERENCES**

We have been unable to find treatments of the problem discussed above in the literature, at least not in the general setting in which we are interested. Similar problems have been studied in the papers listed below.

[1] J.R. Cannon, Determination of certain parameters in heat conduction problems. *J. Math. Anal. Appl.*, **8** (1964), 188-201. [ Links ]

[2] J.R. Cannon, Determination of an unknown coefficient in a parabolic differential equation. *Duke Math. J.*, **30** (1963), 313-323. [ Links ]

[3] J.R. Cannon, Determination of the unknown coefficient* k*(*u*) in the equation *V*·*k*(*u*) Ñ*u* = 0 from overspecified boundary data. *J. Math. Anal. Appl.*, ** 18** (1967), 112-114. [ Links ]

[4] J.R.Cannon and Paul Du Chateau, Determining unknown coefficients in a non-linear heat conduction problem. *SIAM J. Appl. Math.*, **24**(3) (1973), 298-314. [ Links ]

[5] Jim Douglas, Jr and G.F. Jones, Determination of a coefficient in a parabolic differential equation, Part. II: numerical approximation. *J. Math. Mech.*, **11** (1962), 919-926. [ Links ]

[6] B.F. Jones, Various methods for finding unknown coefficients in parabolic differential equations. *Comm. Pure Appl. Math.*, **16** (1963), 33-44. [ Links ]

Received: 14/XI/06

#693/06.

* We thank the Brazilian Mathematical Society (SBM) for the permission to reproduce this article originally published at the "Seminar on Numerical Analysis and its Application to Continuum Physics'' in the ATAS of SBM (Rio de Janeiro), 1980, pp. 65-73.