CONDITIONAL WELL-POSEDNESS FOR AN ELLIPTIC INVERSE PROBLEM
|
|
- John Ryan
- 5 years ago
- Views:
Transcription
1 CONDITIONAL WELL-POSEDNESS FOR AN ELLIPTIC INVERSE PROBLEM IAN KNOWLES AND MARY A. LARUSSA Abstract. For the elliptic equation (p(x) v) + λq(x)v = f, x R n, the problem of determining when one or more of the coefficient functions p, q, and f depends continuously on the solution functions v = v p,q,f,λ, is considered. 1. Introduction and Scientific Context. According to Hadamard [5, 6], a problem, defined as the obtaining of a solution from some given data, is well-posed if three requirements are met: (1) a solution exists; (2) that solution is unique, and (3) depends continuously on the data. Otherwise the problem is ill-posed. In the context of computational inverse problems it is commonly believed, at least in the wider scientific community, that all inverse problems are ill-posed. This has led to the general perception in these communities that the methods, such as Tychonov regularization, that have been specifically developed for truly ill-posed problems, have to be applied to all inverse problems. One of the aims of this article is to show that in some situations there are alternatives that should be considered. Specifically, consider the elliptic differential equation Lv = (p(x) v) + λq(x)v = f(x), x, (1.1) where λ > 0, and is a bounded domain in R n with a C 2,1 -boundary, n 1, and are real, and q and p satisfy f C 0,1 (), q L (), and p C 0,1 (), (1.2) q(x) 0, p(x) ν > 0, x, (1.3) for some constant ν. It is known [11, 30,V] that, under these conditions, for any φ C 1,1 () the Dirichlet problem formed from (1.1) and the boundary condition v = φ on, (1.4) has a unique solution v = v p,q,f,λ in the Hölder space C 1,1 (). In the proofs that follow, we work with solutions u in a Sobolev space context, specifically the Sobolev space H 1 () = W 1,2 (), which contains C 1,1 (). Having the solution v C 1,1 () provides sufficient smoothness on the boundary for integration by parts, which is an important tool in the proofs that follow. Also, note that since p and f are both Lipschitz continuous (i.e. in C 0,1 ()), they are both bounded, and therefore in L 2 (). We are concerned here with determining precise conditions under which inverse problems related to (1.1) are well-posed: specifically, when are one or more of the coefficients p > 0, q 0, f continuously dependent on the solutions v p,q,f,λ for one or more values of λ? Department of Mathematics, University of Alabama at Birmingham, Birmingham, Alabama, 1
2 2 I. KNOWLES AND M. A. LARUSSA Interest in these particular elliptic inverse problems arises naturally from the study of underground aquifer systems (as well as in oil reservoir simulations) which are often modeled by the diffusion equation: q(x) w t = [p(x) w(x, t)] + R(x, t). (1.5) Here, w represents the (known) piezometric head (well pressure), p > 0 the hydraulic conductivity (or sometimes, for a two dimensional aquifer, the transmissivity), R the recharge, and the function q 0 the storativity of the aquifer (see, for example, [1, 2]). It is well known among hydro-geologists [1, Chapter 8] that the inability to obtain reliable values for the coefficients in (1.5) is a serious impediment to the confident use of such models. An effective new method for the simultaneous recovery of p, q, and R (including its time dependence) from a known w is presented in [8, 9, 14]. This method relies upon first transforming the known solution w of the parabolic equation (1.5) to solutions of the elliptic equation (1.1) via a finite Laplace transform, and then recovering the coefficients by solving the above elliptic inverse problem. During the course of this work we were struck by the fact that the computations appeared to converge steadily, oftentimes over thousands of iterations, to the correct functions, which were typically quite discontinuous. As this was clearly not the behavior one normally associates with ill-posed problems, we conjectured that these were actually well-posed problems, at least under certain conditions. As existence is evident from the formulation, and uniqueness for these elliptic inverse problems is treated in [7], our main goal here is to present conditions under which these inverse recoveries exhibit continuous dependence on the data. We assume throughout that uniqueness holds for each of the inverse problems considered in the paper. 2. Preliminary Theory. For later use we need a new trace estimate for the conormal derivatives of H 1 () functions. Such an estimate is known for H 2 () functions (see for example [4, Theorem ]), but when u H 1 (), even the definition of the co-normal derivative of u on is initially problematical. For an open bounded region R n we define the boundary space H 1/2 ( ) to be H 1 ()/H 1 0(), and denote the dual space H 1/2 ( ) by H 1/2 ( ). It is known [12, A.6] that if we define W to be the subspace of all w in H 1 () such that (I )w = 0 distributionally, then we have, with respect to the usual inner product on H 1 (), that H 1 () = H 1 0() W. From this one can also see that the trace map γ : W H 1/2 ( ) defined by where γw = {w} + H 1 0(), γw H 1/2 ( ) = inf{ v H 1 () : v = w + v 0 γw} = w H 1 () is a linear isometry of W onto H 1/2 ( ), and thus an isometric isomorphism. For later use, note also that, for any u H 1 (), u = w + u 0, where w W and u 0 H 1 0(), and we have the standard trace inequality γu H 1/2 ( ) = γw H 1/2 ( ) = w H1 () u H1 (). (2.1)
3 CONDITIONAL WELL-POSEDNESS 3 We henceforth regard u as γu. Now, for u H 1 () and p L () satisfying (1.3), we can define the continuous linear functional f = (I p )u H 1 0() = H 1 () by f(φ) = p u φ + uφ (2.2) for all φ H 1 0(). By the Hahn-Banach theorem f can be extended to a linear functional F that is continuous on H 1 () H 1 0() with F H1 () = f H 1 (). (2.3) For any F continuously extending f, define the functional g H 1/2 ( ) by g(γφ) = F (φ) + p u φ + uφ, (2.4) for all φ H 1 (). It follows from (2.4) and the Lax-Milgram lemma that, as g H 1/2 ( ) and F H 1 (), the weak Neumann problem (I p )u = F (2.5) p ν u = g, (2.6) has u H 1 () as the unique solution. So, each choice of F H 1 () extending f gives rise to a corresponding co-normal derivative g defined by (2.4), and vice versa. The key point here is that once F (extending f) and g are chosen so that (2.4) holds, the F and g uniquely associate with u. Notice that for all F extending a given f, we have that their norms in H 1 () satisfy F H 1 () = sup φ H 1 () φ H 1 () =1 F (φ) sup f(φ) = f H 1(). (2.7) φ H 0 1() φ H 1 0 () =1 Also, by the Riesz representation theorem, if we define an equivalent inner product (and associated norm) on H 1 () (and H 1 0()) by (v, w) 1,p = p v w + vw then, as H 1 0() is a Hilbert space, there is a unique w H 1 0() such that for all φ H 1 0(), f(φ) = (w, φ) 1,p (2.8) and f H 1 () = w 1,p. By the same reasoning, there is also a unique v H 1 () such that for all φ H 1 () and F (φ) = (v, φ) 1,p (2.9) v 1,p = F H1 () = f H 1 () = w 1,p. (2.10)
4 4 I. KNOWLES AND M. A. LARUSSA Setting φ = w in (2.9) and (2.8) we have (v, w) 1,p = w 2 1,p. Assuming without loss that v 0 (and hence w 0) and using (2.10) we have ( ) v w v, = 1, 1,p w 1,p from which it follows that v = w. This means that among all the continuous extensions F of a given f, the Hahn-Banach extension F is uniquely specified by the condition F H1 () = f H 1 (). Notice also that from (2.4) that the unique g associated with F (and u) via (2.4) satisfies for all φ H 1 (), so that 1,p g (γφ) = (u v, φ) 1,p = (u w, φ) 1,p (2.11) g 2 H 1/2 ( ) = u w 2 1,p. (2.12) Observe that for any F H 1 () and g H 1/2 ( ) the weak Neumann problem (I p )u = F (2.13) p ν u = g, (2.14) has a unique solution u H 1 (), by the Lax-Milgram lemma. We partition H 1 () H 1/2 ( ) into equivalence classes in which (F 1, g 1 ) equivalent to (F 2, g 2 ) means that they associate with the same u H 1 (). Motivated by (2.3), (2.7) and (2.12), we define the norm of the equivalence class [(F, g)] generated by u to be [(F, g)] = f H 1 () + u w 2 1,p, (2.15) and denote the set of all such equivalence classes by B. Consider a Cauchy sequence {[(F n, g n )]} in B, where the pair (F n, g n ) associates with u n H 1 (). Let f n be defined by (2.2) with u n replacing u, let w n be defined by (2.8) with f n replacing f. As {[(F n, g n )]} is Cauchy, it follows from (2.15) that the sequence {f n } in H 1 () is Cauchy, hence convergent to f H 1 (). Consequently, from (2.8) the sequence {w n } is Cauchy, hence convergent to w H 1 0(). As the sequence {u n w n } in H 1 () is also Cauchy from (2.15) again, and u n = (u n w n ) + w n, we see that the sequence {u n } is Cauchy, and hence convergent to u H 1 (). Define F H 1 () by F = w and g H 1/2 ( ) by (2.11). Then [(F, g )] is the equivalence class generated by u. It follows that {[(F n, g n )]} {[(F, g )]} in B, and that B is a Hilbert space. So, for each u H 1 () we associate the equivalence class [(F, g)] B, and each such class contains a distinguished pair (F, g ), where F is the unique Hahn- Banach extension associated with f = (I p )u H 1 (), and g is defined by (2.11). For any u H 1 () we define the co-normal derivative of u on by p ν u = g H 1/2 ( ) and (I p )u = F H 1 (). Finally, note that for smooth functions u, such as u C 1,1 () for example, and assuming p C 0,1 (), the co-normal derivative p ν u and the function (I p )u exist independently, and it follows from (2.13) and (2.14) that there is a unique F extending f in this case, and the aforementioned equivalence class associated with u becomes a singleton set. We now have
5 CONDITIONAL WELL-POSEDNESS 5 Theorem 2.1. If is Lipschitz and p ν w denotes the co-normal derivative of w H 1 () on and p C 0,1 () satisfies (1.3) then the mapping w [((I p )w, p ν w)] is an isomorphism of H 1 () onto the Hilbert space B, where B is the Hilbert space of equivalence classes defined above. In particular there exists a constant C, independent of w H 1 (), such that p ν w H 1/2 ( ) C w H1 (). (2.16) Proof. From the discussion above, we see that this mapping is defined on H 1 (). By the Lax-Milgram lemma, for any F H 1 () and any g H 1/2 ( ) the weak Neumann problem (I p )v = F p ν v = g has a unique weak solution v H 1 (), so that the mapping is onto. As is continuous, by [3, V.4.7] the restrictions to of functions in C0 (R n ) are dense in H 1 (). Let u H 1 (). Then there is a sequence {u n }, of restrictions to of C0 (R n ) functions, with limit u in H 1 (). Using Green s theorem, we have that u 2 n + p u n 2 = p ν u n.γu n + (I p )u n.u n = g n(γu n ) + F n(u n ), (2.17) where (Fn, gn) is the unique representative in the singleton equivalence class generated by u n, i.e. in this case gn(γφ) = p ν u n.γφ and Following (2.2) define f n H 1 () by Fn(φ) = (I p )u n.φ f n (φ) = (u n, φ) 1,p = (w n, φ) 1,p (2.18) for some unique w n H 1 0() and all φ H 1 0(). Notice that w n 1,p = f n H 1 (), for all n. As F n is the Hahn-Banach extension of f n, from the preceding discussion we have that F n = w n, together with g n(γφ) = (u n w n, φ) 1,p for all φ H 1 (). From (2.18) we have for all m, n that f n f m H 1 () = sup φ H 1 0 () φ H 1 () =1 sup φ H 1 () φ H 1 () =1 = u n u m 1,p. (f n f m )(φ) (u n u m, φ) 1,p
6 6 I. KNOWLES AND M. A. LARUSSA Consequently, as the sequence {u n } converges strongly to u in H 1 (), the sequence {f n } is Cauchy, and hence is convergent to an f H 1 (). As w n w m 1,p = f n f m H 1 () for all m, n, we also have that the sequence {w n } converges to a function w H 1 0(). Finally we see that the sequence {gn} converges to g H 1/2 ( ) where g (γφ) = (u w, φ) 1,p for all φ H 1 (), and {Fn} (= {w n }) converges to w = F where F H1 () = w 1,p = f H 1 (). Consequently, on letting n in (2.17) one obtains u 2 + p u 2 = g (γu) + F (u). From (2.1) so that, if ν = min(ν, 1), Thus γu H 1/2 ( ) u H 1 (), ν u 2 H 1 () g H 1/2 ( ). γu H 1/2 ( ) + F H 1 (). u H 1 () g H 1/2 ( ). u H 1 () + F H 1 (). u H 1 (). ν u H 1 () g H 1/2 ( ) + F H 1 () = [(F, g )] B So the map in question, in addition to being onto, is closed and injective, and has a continuous inverse; its continuity follows by the closed graph theorem. We conclude this section with a lemma that will find frequent use in the sequel. Lemma 2.2. Let H be a Hilbert space and let T : D H H be a densely defined positive linear selfadjoint operator in H. Assume also that T 1 : H D is compact. If {g n } converges to zero weakly in H, and {T v n } is bounded in H, then {(g n, v n ) H } converges to zero. Proof. Let y n = T 1 g n for all n. By [13, Theorem 6.3] and the compactness of T 1 and the weak convergence of {g n }, the sequence {y n } converges strongly to zero in H. Then (g n, v n ) H = (T y n, v n ) H = (y n, T v n ) H. From the foregoing and the Cauchy-Schwarz inequality, this goes to zero as n. 3. The One-coefficient Case. Consider first that case in which p, f, and λ > 0 in (1.1) are given, and the coefficient q is to be recovered from the known solution u p,q,f,λ H 1 (). The recovery method discussed earlier makes use of a gradient descent technique in which a sequence of iterates {q n } is produced with the property that the corresponding sequence of solutions {u p,qn,f,λ} converges to u p,q,f,λ strongly in H 1 () (c.f. [10, Theorems 5.1 and 5.3]). In the next theorem we connect the convergence of the solutions {u p,qn,f,λ} with the convergence of their inversions {q n }. Analogous statements may also be made about each of the subsequent theorems in the paper. Theorem 3.1. Assume that λ > 0, p, and f satisfying (1.2) and (1.3) are given, and for some constant M q 0, 0 q n M, for all n. (3.1)
7 CONDITIONAL WELL-POSEDNESS 7 Assume also that u = u p,q,f,λ is a given solution of (1.1) for which u p,q,f,λ c > 0 on (3.2) for some constant c. Let u p,qn,f,λ denote solutions of (1.1), with q replaced by q n, that have the same boundary data as u. Then the sequence {q n } converges weakly to q in L 2 () if and only if the sequence {u p,qn,f,λ} converges to u p,q,f,λ in H 1 (). Proof. Let u n = u p,qn,f,λ, and assume first that u n u in H 1 (). Subtracting the equations satisfied by u p,q,f,λ and u p,qn,f,λ gives λ(q n q)u = p (u u n ) + λq n (u u n ), (3.3) where, by (1.2), all the terms above are functions in L (). Let φ C 0 () be chosen arbitrarily. Noting that ψ = φ /(λu) H 1 0() and applying Green s theorem, we have (q n q)φ = = p (u u n )ψ + λ p (u u n ). ψ + λ q n (u u n )ψ q n (u u n )ψ, which approaches zero as n. Let φ L 2 () and ɛ > 0 be given, and choose φ C 0 () so that φ φ L2 () < ɛ 4M.m() 1/2, where m() denotes the Lebesgue measure of the bounded open set. Then choose N so that for n > N, (q q n )φ < ɛ 2. Consequently, for n > N (q q n )φ (q q n )φ + 2M. φ φ < ɛ. As φ was chosen arbitrarily, we have the desired weak convergence in L 2 (). Conversely, assume that q n q weakly in L 2 (). After multiplying (3.3) by u u n H 1 0() and using Green s Theorem, we have that p (u u n ) 2 + λq n (u u n ) 2 = λ u.(q n q).(u u n ). (3.4) Using the assumptions that p ν, λ > 0, and q n 0 we have p (u u n ) 2 + λq n (u u n ) 2 ν (u u n ) ν (u u n ) νλ 0 u u n ν. min(1, λ 0). u u n 2 H 1 (), (3.5)
8 8 I. KNOWLES AND M. A. LARUSSA where λ 0 is the smallest positive eigenvalue of the negative Laplacian on. Now, from (3.4), (3.5), and q n q 2M, and the inequality 2ab ɛa 2 + (1/ɛ)b 2, for any ɛ > 0 we have u u n 2 H K 1 u(q 1 n q)(u u n ) (3.6) () Choose ɛ so that K 2 ɛ = 1 2. It follows that for all n. Hence, K 2 (ɛ u u n 2 L ɛ u 2 L 2) (3.7) K 2 ɛ u u n 2 H 1 () + K 2 ɛ u 2 L 2 (). (3.8) u u n L 2 () u u n H 1 () K 3, u n L 2 () u L 2 () + u u n L 2 () (3.9) u L 2 () + K 3, (3.10) for all n, and the sequence {u n } is bounded in L 2 (). Finally, the linear operator L = p + λq defined on the dense domain D = {v H 1 0() : Lv L 2 ()} in L 2 () is selfadjoint and positive, and has a compact inverse operator with domain L 2 (), by [3, Theorem 1.4, p. 306]. From (3.3) we have that the sequence {g n }, where g n = p (u u n ) + λq n (u u n ) = λ(q n q)u, converges weakly to zero in L 2 (), and from a calculation, L(u u n ) = λ(q n q)u n, so that, by (3.10), {L(u u n )} is bounded in L 2 (). Finally, from (3.5) we have u u n 2 H 1 () K 4. p (u u n ) 2 + λq n u u n 2 = K 4.(g n, u u n ) L2 (), which approaches zero as n by Lemma 2.2. It follows that u n u in H 1 (), as promised. Remarks. 1. The inverse problem may be written as q = F (u), where the solution u = u p,q,f,λ, A p,f,λ,m = {v H 1 () : v satisfies (1.1) for some q L (), 0 q M }, and F : A p,f,λ,m L () L 2 (). While the set A p,f,λ,m is not a linear space, it can be given a metric topology via the metric d defined by d(v, w) = v w H 1 (), v, w A p,f,λ,m. If the range of F is considered to be a subset of the space L 2 (), then we can infer from the above theorem that the inverse problem of obtaining q from u is well-posed
9 CONDITIONAL WELL-POSEDNESS 9 Fig Recovery of q from u p,q,λ Z Z X Y X Y (a) True q (b) Computed q in the sense that the continuity required of part (3) of the definition of well-posedness is taken in the sense of the above metric topology on A p,f,λ,m and the weak topology of L 2 () on the range of F. 2. A typical recovery of a discontinuous coefficient q from a known solution u (using the method of [8]) is illustrated in figure 3. Notice that while the location and height of the discontinuity is well recovered here, the pointwise errors ( jaggies ) are quite noticeable. This is in line with the weak convergence proven above. In fact if, for any set one chooses φ L 2 () to be the characteristic function of that set, then we have that (q n q)(x) dx = (q n q)(x)φ(x) dx 0 as n. This indicates that observed errors in the recovery, (q q n )(x), are small on average over any subset, rather than necessarily small in a pointwise sense at a given x. Notice also that while the q n approximate q only in a weak sense, the corresponding functions u n are strongly convergent to u in H 1 (). This shows that the modeling with the inversely recovered q is effective in the sense that the the data, u, is faithfully reproduced when one solves equation (1.1) using the recovered q. Similar comments apply to the remaining theorems in the paper. 3. The condition (3.2) is satisfied naturally in practice. Our solutions u p,q,f,λ (x, λ) = 1 0 w(x, t)e λt dt, x, are finite Laplace transforms of piezometric head values w(x, t). Here, one normally knows that w(x, t) β > 0 uniformly in x and t, and one can easily infer (3.2) from this. We consider next the situation in which p, q, and λ > 0 and u p,q,f,λ are given and the continuous recovery of f is desired. Theorem 3.2. Assume that λ > 0, p and q satisfy (1.2) and (1.3), and that there exists a constant K such that, f L 2 () K, and f n L 2 () K, for all n. (3.11) Assume also that u = u p,q,f,λ is a given solution of (1.1). Let u p,q,fn,λ denote solutions of (1.1), with f replaced by f n, that have the same boundary data as u. Then the
10 10 I. KNOWLES AND M. A. LARUSSA sequence {f n } converges weakly to f in L 2 () if and only if the sequence {u p,q,fn,λ} converges to u p,q,f,λ in H 1 (). Proof. Let u n = u p,q,fn,λ, and assume first that u n u in H 1 (). Subtracting the equations satisfied by u p,q,f,λ and u p,q,fn,λ we have f f n = p (u u n ) + λq(u u n ). (3.12) For ψ C 0 () we have (f f n )ψ = p (u u n )ψ + λ q(u u n )ψ = p (u u n ). ψ + λ q(u u n )ψ, which approaches zero as n. Thus, we have the desired weak convergence in L 2 () by a similar argument to that used in the previous theorem. Conversely, assume that f n f weakly in L 2 (). From (3.12) we have that the sequence {g n }, where g n = L(u u n ) = p (u u n ) + λq(u u n ) = f f n, converges weakly to zero in L 2 (). The sequence {L(u u n )} is also bounded in L 2 (). Thus, by Lemma 2.2 it follows as before that u n u in H 1 (). If q, f, and λ > 0 and u p,q,f,λ are given then we have the following theorem on the continuous recovery of p. Theorem 3.3. Assume that λ > 0, q, and f satisfying (1.2) and (1.3) are given. Assume also that u = u p,q,f,λ is a given solution of (1.1) with the property that for each ψ L 2 () there is a φ H 1 () solving the first order equation φ u p,q,f,λ = ψ. (3.13) Let u pn,q,f,λ denote solutions of (1.1), with p replaced by p n, that have the same boundary data as u, and assume that 0 < ν p n K, for all n, and (p p n ) = 0, for all n. (3.14) If the sequence {u pn,q,f,λ} converges to u p,q,f,λ in H 1 () then the sequence {p n } converges weakly to p in L 2 (). Proof. Let u n = u pn,q,f,λ, and assume that u n u in H 1 (). Subtracting the equations satisfied by u p,q,f,λ and u pn,q,f,λ we have (p p n ) u = p n (u u n ) + λq(u u n ). (3.15) For an arbitrary ψ L 2 () let φ H 1 () satisfy (3.13). We have, via Green s theorem, (3.14), and (3.15), (p p n )ψ = (p p n ) u φ = φ (p p n ) u = φ p n (u u n ) λq.(u u n ).φ = φp n ν (u u n ) p n (u u n ) φ + λq.(u n u).φ = φp ν (u u n ) p n (u u n ) φ + λq.(u n u).φ.
11 CONDITIONAL WELL-POSEDNESS 11 The second and third terms on the right above approach zero as n from the convergence to zero in H 1 () of the sequence {u u n } and the Cauchy-Schwarz inequality. As for the first term we have, by the trace inequalities (2.1) and (2.16), φ.p ν (u u n ) φ H ( ). p 1/2 ν (u u n ) H 1/2 ( ) C. φ H 1 (). u u n H 1 (), and this term also tends to zero. As ψ was chosen arbitrarily, we have the desired weak convergence in L 2 (). As a partial converse of the previous theorem, we also have Theorem 3.4. Assume that λ > 0, q, and f satisfying (1.2) and (1.3) are given. Assume also that u = u p,q,f,λ is a given solution of (1.1), and let u pn,q,f,λ denote solutions of (1.1), with p replaced by p n, that have the same boundary data as u, and assume that for some constant K, 0 < ν p n K and p n K for all n, and that (3.14) holds. If (p p n ) u p,q,f,λ φ 0 (3.16) for all φ H 1 (), then {u pn,q,f,λ} converges to u p,q,f,λ in H 1 (). Proof. Assume that (3.16) holds. From (3.15) we have, as before by first considering φ C 0 () and then using a density argument, that the sequence {g n } defined by g n = p n (u u n ) + λq(u u n ) = (p p n ) u, converges weakly to zero in L 2 (). Now L(u u n ) = p (u u n ) + λq(u u n ) = p u n + p u n p u + λq(u u n ), and from the differential equation satisfied by u n u n = ( p n u n + λqu n f)/p n. So, from the assumptions on p n and p n and the boundedness of {u n } in H 1 () (proved as before), the sequence {L(u u n )} is bounded in L 2 (). By Lemma 2.2 it follows that (g n, u u n ) L2 () 0 and thus that u n u in H 1 () as required. Remarks. 1. By the method of characteristics, the first order linear partial differential equation (3.13) will have solutions φ when the flow lines of the vector field u p,q,f,λ are such that any point in exits to the boundary along a flow line. Given that the function u p,q,f,λ forms part of the data of the inverse problem, one can in principle at least check this approximately via a computer plot of the flow field of u p,q,f,λ. 2. If we forego the condition (3.13), combining the remaining conditions of the last two theorems one may also prove, using the same method, that u n u 0 in H 1 0() if and only if for all φ H 1 () (p p n ) u φ 0.
12 12 I. KNOWLES AND M. A. LARUSSA One may consider this last statement as a kind of weighted weak convergence of the sequence {p n }. 3. One should note that uniqueness for the inverse problem is guaranteed if one specifies the coefficient p on the boundary [7]. Because of this, the descent methods mentioned above all enforce the same boundary data on the iterates that converge to p by using Neuberger gradients in the descent processes. The condition (3.14) in the theorem above is thus automatically true for these inverse methods. 4. The Two-coefficient Case Recovering q and f with p known. Consider first the case in which p is known and q and f are to be recovered from known solutions u 1 = u p,q,f,λ1 and u 2 = u p,q,f,λ2, where λ 1 λ 2. Theorem 4.1. With p and λ 2 > λ 1 > 0 satisfying (1.2) and (1.3) given, assume that for some constants K and M, f n L2 () K, 0 q n M, for all n. (4.1) Assume also that u i = u p,q,f,λi, i = 1, 2, are given solutions of (1.1) for which λ 1 u p,q,f,λ1 λ 2 u p,q,f,λ2 c > 0, (4.2) for some constant c. Let u p,qn,f n,λ i denote solutions of (1.1), with q replaced by q n and f replaced by f n, that have the same boundary data as u i. Then the sequences {q n } and {f n } converge weakly to q and f, respectively, in L 2 () if and only if for i = 1, 2 the sequences {u p,qn,f n,λ i } converge to u p,q,f,λi in H 1 (). Proof. Assume that for i = 1, 2 the sequences {u p,qn,f n,λ i } converge to u p,q,f,λi in H 1 (). We have and Subtracting these we get p u p,q,f,λ1 + λ 1 qu p,q,f,λ1 = f (4.3) p u p,q,f,λ2 + λ 2 qu p,q,f,λ2 = f (4.4) p u p,qn,f n,λ 1 + λ 1 q n u p,qn,f n,λ 1 = f n (4.5) p u p,qn,f n,λ 2 + λ 2 q n u p,qn,f n,λ 2 = f n. (4.6) p (u p,q,f,λ1 u p,q,f,λ2 ) + q(λ 1 u p,q,f,λ1 λ 2 u p,q,f,λ2 ) = 0 (4.7) p (u p,qn,f n,λ 1 u p,qn,f n,λ 2 ) + q n (λ 1 u p,qn,f n,λ 1 λ 2 u p,qn,f n,λ 2 ) = 0. (4.8) If we set v i,n = u p,q,f,λi u p,qn,f n,λ i, subtracting the last two equations then gives p [v 1,n v 2,n ] + q n [λ 1 v 1,n λ 2 v 2,n ] = (q n q)(λ 1 u p,q,f,λ1 λ 2 u p,q,f,λ2 ) (4.9) As v i,n 0 in H 1 () for i = 1, 2 the proof that the sequence {q n } converges weakly to q in L 2 () now follows by the corresponding argument in the one coefficient case, using (4.2). Subtracting (4.3) and (4.5) gives p v 1,n + λ 1 [(q q n )u p,q,f,λ1 + q n v 1,n ] = f f n, (4.10)
13 CONDITIONAL WELL-POSEDNESS 13 and we have also that {f n } converges weakly to f in L 2 (). Conversely, assume that the sequences {q n } and {f n } converge weakly to q and f, respectively, in L 2 (). On rearranging (4.10) we find that From (4.11) we have that p v 1,n + λ 1 q n v 1,n = λ 1 (q n q)u p,q,f,λ1 + (f f n ). (4.11) g n = p v 1,n + λ 1 q n v 1,n (4.12) converges weakly to zero in L 2 (). Note that from (4.11) the functions Lv 1,n = p v 1,n + λ 1 q.v 1,n = λ 1 (q n q)u p,qn,f n,λ 1 + (f f n ) form a bounded set in L 2 (). It follows from Lemma 2.2 that {u p,qn,f n,λ 1 } converges to u p,q,f,λ1 in H 1 (). The convergence of the sequence {u p,qn,f n,λ 2 } is handled similarly. Remark. We may write λ 1 u p,q,f,λ1 λ 2 u p,q,f,λ2 = 1 0 (λ 1 e λ1t λ 2 e λ2t )w(x, t) dt, where the w(x, t) values are piezometric head measurements. equation λ 1 e λ1t λ 2 e λ2t = 0 has a unique root t 0 = λ2 ln( λ 1 ) λ 1 ( λ2 λ 1 1) If we note that the then if we choose λ 1, λ 2 so that t 0 > 1 then the factor λ 1 e λ1t λ 2 e λ2t 0 for t [0, 1]. As the head measurements w(x, t) are usually uniformly bounded below by a positive number, one can see that the restriction (4.2) is quite natural here Recovering p and f with q known. Consider next the case in which q is known and p and f are to be recovered from known solutions u 1 = u p,q,f,λ1 and u 2 = u p,q,f,λ2, λ 1 λ 2. Theorem 4.2. With q and λ 2 > λ 1 > 0 satisfying (1.2) and (1.3) given, assume that for some constants K and M, 0 < ν p n K, f n L 2 () K, for all n. (4.13) Assume also that each p n satisfies (1.3) and that for all n (p p n ) = 0, (4.14) and u i = u p,q,f,λi, i = 1, 2, are given solutions of (1.1) with u p,q,f,λi c i > 0, for some constants c i, i = 1, 2, and for any ψ L 2 (), there is a solution φ H 1 () of the first order equation φ (u p,q,f,λ1 u p,q,f,λ2 ) = ψ. (4.15) Let u pn,q,f n,λ i denote solutions of (1.1), with p replaced by p n and f replaced by f n, that have the same boundary data as u p,q,f,λi. If, for i = 1, 2 the sequences
14 14 I. KNOWLES AND M. A. LARUSSA {u pn,q,f n,λ i } converge, respectively, to u p,q,f,λi in H 1 (), then the sequences {p n } and {f n } converge weakly to p and f, respectively, in L 2 (). Proof. Assume that for i = 1, 2 the sequences {u pn,q,f n,λ i } converge to u p,q,f,λi in H 1 (). We have and Subtracting these we get p u p,q,f,λ1 + λ 1 qu p,q,f,λ1 = f (4.16) p u p,q,f,λ2 + λ 2 qu p,q,f,λ2 = f (4.17) p n u pn,q,f n,λ 1 + λ 1 qu pn,q,f n,λ 1 = f n (4.18) p n u pn,q,f n,λ 2 + λ 2 qu pn,q,f n,λ 2 = f n. (4.19) p (u p,q,f,λ1 u p,q,f,λ2 ) + q(λ 1 u p,q,f,λ1 λ 2 u p,q,f,λ2 ) = 0 (4.20) p n (u pn,q,f n,λ 1 u pn,q,f n,λ 2 ) + q(λ 1 u pn,q,f n,λ 1 λ 2 u pn,q,f n,λ 2 ) = 0. (4.21) If we set v i,n = u p,q,f,λi u pn,q,f n,λ i, i = 1, 2, subtracting the last two equations then gives (p n p) (u p,q,f,λ1 u p,q,f,λ2 ) = p n (v 1,n v 2,n )+q[λ 1 v 1,n λ 2 v 2,n ] (4.22) Consequently, multiplying by φ H 1 () and applying Green s theorem, using (4.14) and the trace estimate (2.16) on the boundary term on the right, gives (p n p) (u p,q,f,λ1 u p,q,f,λ2 ) φ 0 for all φ in H 1 (). It follows from (4.15) that {p n } converges weakly in L 2 (). If we now subtract (4.16) and (4.18) we see that f f n = (p p n ) u p,q,f,λ1 p n v 1,n + λ 1 qv 1,n. It follows that the sequence {f n } converges weakly to f in L 2 (). We also have the following partial converse. Theorem 4.3. With q and λ 2 > λ 1 > 0 given, assume that for some constants K and M, 0 < ν p K, f L2 () K, 0 q M and 0 < ν p n K, p n K, f n L 2 () K, for all n. (4.23) Assume also that each p n satisfies (1.3) and that for all n (4.14) holds and u i = u p,q,f,λi, i = 1, 2, are given solutions of (1.1) with u p,q,f,λi c i > 0, for some constants c i, i = 1, 2. Let u pn,q,f n,λ i denote solutions of (1.1), with p replaced by p n and f replaced by f n, that have the same boundary data as u p,q,f,λi. If the sequence {f n } converges weakly to f in L 2 () and for i = 1, 2 (p p n ) u p,q,f,λi φ 0 (4.24) for all φ H 1 () then for i = 1, 2 the sequences {u pn,q,f n,λ i } converge, respectively, to u p,q,f,λi in H 1 ().
15 CONDITIONAL WELL-POSEDNESS 15 Proof. Assume that the sequence {f n } converges weakly to f in L 2 () and for i = 1, 2 and every φ H 1 () (p p n ) u p,q,f,λi φ 0. A further subtraction gives g n = p n v 1,n + λ 1 qv 1,n = f f n + (p p n ) u p,q,f,λ1, (4.25) so that the sequence {g n } converges weakly to zero in L 2 (). Also, Lv 1,n = p v 1,n + λ 1 qv 1,n = f λ 1 qu p,q,f,λ1 + p u pn,q,f n,λ 1 + λ 1 qv 1,n = f λ 1 qu p,q,f,λ1 + p u pn,q,f n,λ 1 + p u pn,q,f n,λ 1 + λ 1 qv 1,n, and from the equation for u pn,q,f n,λ 1 we have u pn,q,f n,λ 1 = ( p n u pn,q,f n,λ 1 + λ 1 qu pn,q,f n,λ 1 f n )/p n. So, from the boundedness of {u pn,q,f n,λ 1 } in H 1 () (proven as before) and the assumptions on p n we see that the sequence {Lv 1,n } is bounded in L 2 (). It follows from Lemma 2.2 that {u pn,q,f n,λ 1 } converges to u p,q,f,λ1 in H 1 (); the convergence of {u pn,q,f n,λ 2 } to u p,q,f,λ2 is similar Recovering p and q with f known. Consider next the case in which f is known and p and q are to be recovered from known solutions u 1 = u p,q,f,λ1 and u 2 = u p,q,f,λ2, λ 1 λ 2. Theorem 4.4. With f and λ 2 > λ 1 > 0 given, assume that for some constants K and M, 0 < ν p K, 0 q M, 0 < ν p n K, 0 q n M, for all n. (4.26) Assume also that each p n satisfies (1.3) and that for all n (p p n ) = 0, (4.27) and u i = u p,q,f,λi, i = 1, 2, are given solutions of (1.1) with u p,q,f,λi c i > 0 and with the property that for any ψ L 2 () the first order equation (λ 2 u p,q,f,λ2 u p,q,f,λ1 λ 1 u p,q,f,λ1 u p,q,f,λ2 ) φ + (λ 2 λ 1 )φ( u p,q,f,λ2 u p,q,f,λ1 ) = ψ (4.28) has a solution φ H 1 (). Let u pn,q n,f,λ i denote solutions of (1.1), with p replaced by p n and q replaced by q n, that have the same boundary data as u p,q,f,λi. If for i = 1, 2 the sequences {u pn,q n,f,λ i } converge to u p,q,f,λi in H 1 () then the sequences {p n } and {q n } converge weakly to p and q, respectively, in L 2 (). Proof. Assume that for i = 1, 2 the sequences {u pn,q n,f,λ i } converge to u p,q,f,λi in H 1 (). We have p u p,q,f,λ1 + λ 1 qu p,q,f,λ1 = f (4.29) p u p,q,f,λ2 + λ 2 qu p,q,f,λ2 = f (4.30)
16 16 I. KNOWLES AND M. A. LARUSSA and p n u pn,q n,f,λ 1 + λ 1 q n u pn,q n,f,λ 1 = f (4.31) p n u pn,q n,f,λ 2 + λ 2 q n u pn,q n,f,λ 2 = f. (4.32) Combining these we get λ 2 u p,q,f,λ2 p u p,q,f,λ1 + λ 1 u p,q,f,λ1 p u p,q,f,λ2 = (λ 2 u p,q,f,λ2 λ 1 u p,q,f,λ1 ).f (4.33) λ 2 u pn,q n,f,λ 2 p n u pn,q n,f,λ 1 + λ 1 u pn,q n,f,λ 1 p n u pn,q n,f,λ 2 = (λ 2 u pn,q n,f,λ 2 λ 1 u pn,q n,f,λ 1 ).f (4.34) If we set v i,n = u p,q,f,λi u pn,q n,f,λ i, i = 1, 2, combining the last two equations then gives Here λ 2 u p,q,f,λ2 (p p n ) u p,q,f,λ1 + λ 1 u p,q,f,λ1 (p p n ) u p,q,f,λ2 = λ 2 u p,q,f,λ2 p n u p,q,f,λ1 λ 2 u pn,q n,f,λ 2 p n u pn,q n,f,λ 1 λ 1 u p,q,f,λ1 p n u p,q,f,λ2 + λ 1 u pn,q n,f,λ 1 p n u pn,q n,f,λ 2 + [(λ 2 u p,q,f,λ2 λ 1 u p,q,f,λ1 ) (λ 2 u pn,q n,f,λ 2 λ 1 u pn,q n,f,λ 1 )].f (4.35) λ 2 u p,q,f,λ2 p n u p,q,f,λ1 λ 2 u pn,q n,f,λ 2 p n u pn,q n,f,λ 1 = λ 2 (v 2,n p n u p,q,f,λ1 + u pn,q n,f,λ 2 p n v 1,n ), (4.36) and λ 1 u p,q,f,λ1 p n u p,q,f,λ2 λ 1 u pn,q n,f,λ 1 p n u pn,q n,f,λ 2 = λ 1 (v 1,n p n u p,q,f,λ2 + u pn,q n,f,λ 1 p n v 2,n ), (4.37) [(λ 2 u p,q,f,λ2 λ 1 u p,q,f,λ1 ) (λ 2 u pn,q n,f,λ 2 λ 1 u pn,q n,f,λ 1 )].f = (λ 2 v 2,n λ 1 v 1,n ).f (4.38) So, it now follows that (p p n )[ u p,q,f,λ1 (λ 2 u p,q,f,λ2 φ) u p,q,f,λ2 (λ 1 u p,q,f,λ1 φ)] 0 and hence that (p p n )[(λ 2 u p,q,f,λ2 u p,q,f,λ1 λ 1 u p,q,f,λ1 u p,q,f,λ2 ) φ + (λ 2 λ 1 )φ u p,q,f,λ2 u p,q,f,λ1 )] 0, for all φ in H 1 (). It follows from (4.28) that the sequence {p n } converges weakly to p in L 2 ().
17 CONDITIONAL WELL-POSEDNESS 17 Finally, subtracting (4.29) and (4.31) we have λ 1 (q n q)u p,q,f,λ1 = p n v 1,n λ 1 q n v 1,n + (p p n ) u p,q,f,λ1, (4.39) The weak convergence of {q n } to q now follows from the weak convergence of {p n } proved above by following the method used at the same stage in the proof of Theorem 3.1. We also have the following partial converse. Theorem 4.5. With f and λ 2 > λ 1 > 0 given, assume that and for some constants K and M, 0 < ν p K, 0 q M, 0 < ν p n K, p n K, 0 q n M, for all n. (4.40) Assume also that each p n satisfies (1.3) and that for all n (4.27) holds and u i = u p,q,f,λi, i = 1, 2, are given solutions of (1.1) with u p,q,f,λi c i > 0. Let u pn,q n,f,λ i denote solutions of (1.1), with p replaced by p n and q replaced by q n, that have the same boundary data as u p,q,f,λi. If {q n } converges weakly to q in L 2 () and for i = 1, 2 (p p n ) u p,q,f,λi φ 0 (4.41) for all φ H 1 (), then for i = 1, 2 the sequences {u pn,q n,f,λ i } converge to u p,q,f,λi in H 1 (). Proof. Assume that the sequence {q n } converges weakly to q in L 2 (), and that (4.41) holds. Rearranging (4.39) gives g n = p n v 1,n + λ 1 q n v 1,n = λ 1 (q q n )u p,q,f,λ1 + (p p n ) u p,q,f,λ1, (4.42) so that the sequence {g n } converges weakly to zero in L 2 (). Also, Lv 1,n = p v 1,n + λ 1 qv 1,n = f λ 1 qu p,q,f,λ1 + p u pn,q n,f,λ 1 + λ 1 qv 1,n = f λ 1 qu p,q,f,λ1 + p u pn,q n,f,λ 1 + p u pn,q n,f,λ 1 + λ 1 qv 1,n, where, from the differential equation for u pn,q n,f,λ 1, u pn,q n,f,λ 1 = ( p u pn,q n,f,λ 1 + λ 1 q n u pn,q n,f,λ 1 f)/p n. From the boundedness of {u pn,q,f n,λ 1 } in H 1 () (proven as before) and the assumptions on p n and q n, the sequence {Lv 1,n } is bounded in L 2 (). So, by Lemma 2.2 again, {u pn,q,f n,λ 1 } converges to u p,q,f,λ1 in H 1 (); the convergence of {u pn,q,f n,λ 2 } to u p,q,f,λ2 is similar. 5. The Three-coefficient Case. In practice, the most relevant case is the one in which p, q, and f are recovered simultaneously. Theorem 5.1. Assume that for some positive numbers λ i, 1 i 3, with λ 3 > λ 2 > λ 1, solutions u i = u p,q,f,λi of the equation (1.1) are given with the property that for any ψ L 2 () the first order equation φ [w 13 (u p,q,f,λ1 u p,q,f,λ2 ) w 12 (u p,q,f,λ1 u p,q,f,λ3 )] + [ w 13 (u p,q,f,λ1 u p,q,f,λ2 ) w 12 (u p,q,f,λ1 u p,q,f,λ3 )].φ = ψ, (5.1)
18 18 I. KNOWLES AND M. A. LARUSSA where for i = 2, 3, w 1i = λ 1 u p,q,f,λ1 λ i u p,q,f,λi, has a solution φ H 1 (). Assume that w 1,i c > 0 for some constant c and i = 2, 3, and that for some constants K and M, 0 < ν p K, f L2 () K, 0 q M and 0 < ν p n K, f n L2 () K, 0 q n M for all n (5.2) Assume also that each p n satisfies (1.3) and that for all n (p p n ) = 0. Let u pn,q n,f n,λ i denote solutions of (1.1), with p replaced by p n, q replaced by q n, and f replaced by f n that have the same boundary data as u p,q,f,λi. If, for 1 i 3 the sequences {u pn,q n,f n,λ i } converge to u p,q,f,λi in H 1 () then the sequences {p n }, {q n }, and {f n } converge weakly to p, q, and f, respectively, in L 2 (). Proof. Assume that for 1 i 3 the sequences {u pn,q n,f,λ i } converge to u p,q,f,λi in H 1 (). We have and Subtracting gives and p u p,q,f,λ1 + λ 1 qu p,q,f,λ1 = f (5.3) p u p,q,f,λ2 + λ 2 qu p,q,f,λ2 = f (5.4) p u p,q,f,λ3 + λ 3 qu p,q,f,λ3 = f (5.5) p n u pn,q n,f n,λ 1 + λ 1 q n u pn,q n,f n,λ 1 = f n (5.6) p n u pn,q n,f n,λ 2 + λ 2 q n u pn,q n,f n,λ 2 = f n (5.7) p n u pn,q n,f n,λ 3 + λ 3 q n u pn,q n,f n,λ 3 = f n. (5.8) p (u p,q,f,λ1 u p,q,f,λ2 ) + q(λ 1 u p,q,f,λ1 λ 2 u p,q,f,λ2 ) = 0 (5.9) p (u p,q,f,λ1 u p,q,f,λ3 ) + q(λ 1 u p,q,f,λ1 λ 3 u p,q,f,λ3 ) = 0 (5.10) p n (u pn,q n,f n,λ 1 u pn,q n,f n,λ 2 ) + q n (λ 1 u pn,q n,f n,λ 1 λ 2 u pn,q n,f n,λ 2 ) = 0 (5.11) p n (u pn,q n,f n,λ 1 u pn,q n,f n,λ 3 ) + q n (λ 1 u pn,q n,f n,λ 1 λ 3 u pn,q n,f n,λ 3 ) = 0. (5.12) If we set v i,n = u p,q,f,λi u pn,q n,f n,λ i, 1 i 3, subtracting again then gives and (p p n ) (u p,q,f,λ1 u p,q,f,λ2 ) p n (v 1,n v 2,n ) = (q n q)(λ 1 u p,q,f,λ1 λ 2 u p,q,f,λ2 ) + q n [λ 2 v 2,n λ 1 v 1,n ] (5.13) (p p n ) (u p,q,f,λ1 u p,q,f,λ3 ) p n (v 1,n v 3,n ) = (q n q)(λ 1 u p,q,f,λ1 λ 3 u p,q,f,λ3 ) + q n [λ 3 v 3,n λ 1 v 1,n ]. (5.14) Finally, eliminating the terms in q q n in the last two equations gives w 13 (p p n ) (u p,q,f,λ1 u p,q,f,λ2 ) + w 12 (p p n ) (u p,q,f,λ1 u p,q,f,λ3 ) = [w 13 p n (v 1,n v 2,n ) w 12 p n (v 1,n v 1,3 )] +q n [w 12 (λ 1 v 1,n λ 3 v 3,n ) w 13 (λ 1 v 1,n λ 2 v 2,n )].
19 It now follows that CONDITIONAL WELL-POSEDNESS 19 (p p n ){ φ [w 13 (u p,q,f,λ1 u p,q,f,λ2 ) w 12 (u p,q,f,λ1 u p,q,f,λ3 )]} + [ w 13 (u p,q,f,λ1 u p,q,f,λ2 ) w 12 (u p,q,f,λ1 u p,q,f,λ3 )].φ 0 (5.15) for all φ H 1 (). From (5.1) we have that the sequence {p n } converges weakly to p in L 2 (). The weak convergence of the sequence {q n } to q is then a consequence of (5.13). Finally, subtracting (5.3) and (5.6) gives f f n = (p p n ) u p,q,f,λ1 +λ 1 (q q n )u p,q,f,λ1 (p n ) v 1,n +λ 1 q n v 1,n, (5.16) from which follows the weak convergence of the sequence {f n } to f. The following partial converse also holds. Theorem 5.2. Assume that for some positive numbers λ i, 1 i 3, with λ 3 > λ 2 > λ 1, solutions u i = u p,q,f,λi of the equation (1.1) are given. Assume that w 1,i c > 0 for some constant c and i = 2, 3, and that for some constants K and M, 0 < ν p K, f L2 () K, 0 q M and 0 < ν p n K, p n K, f n L2 () K, 0 q n M, for all n. (5.17) Assume also that each p n satisfies (1.3) and that for all n (p p n ) = 0. Let u pn,q n,f n,λ i denote solutions of (1.1), with p replaced by p n, q replaced by q n, and f replaced by f n that have the same boundary data as u p,q,f,λi. If {q n } converges weakly to q and {f n } converges weakly to f in L 2 (), and for i = 1, 2, 3 (p p n ) u p,q,f,λi φ 0 (5.18) for all φ H 1 (), then for i = 1, 2, 3 the sequence {u pn,q n,f,λ i } converges to u p,q,f,λi in H 1 (). Proof. Set g n = p n v 1,n + λ 1 q n v 1,n = λ 1 (q n q)u p,q,f,λ1 + (p p n ) u p,q,f,λ1 + f f n, (5.19) so that the sequence {g n } converges weakly to zero in L 2 (). Also, Lv 1,n = p v 1,n + λ 1 qv 1,n = f λ 1 qu pn,q n,f n,λ 1 + p u pn,q n,f n,λ 1 + λ 1 qv 1,n = p u pn,q n,f n,λ 1 + p u pn,q n,f n,λ 1 λ 1 qu pn,q n,f n,λ 1 + λ 1 qv 1,n + f and the sequence {Lv 1,n } is bounded in L 2 (). So by Lemma 2.2 {u pn,q n,f n,λ 1 } converges to u p,q,f,λ1 in H 1 (); the convergence of {u pn,q,f n,λ i } to u p,q,f,λi, i = 2, 3, is again similar.
20 20 I. KNOWLES AND M. A. LARUSSA Acknowledgment. The authors are grateful to the referees for their careful reading of the original manuscript and catching several errors in the process, and for other helpful suggestions on improving the readability of the paper. REFERENCES [1] Mary P. Anderson and William W. Woessner. Applied Groundwater Modeling. Academic Press, New York, [2] J. Bear. Dynamics of Fluids in Porous Media. American Elsevier, New York, [3] D. E. Edmunds and W. D. Evans. Spectral Theory and Differential Operators. Oxford Mathematical Monographs. Clarendon Press, Oxford, [4] P. Grisvard. Elliptic problems in nonsmooth domains, volume 24 of Monographs and Studies in Mathematics. Pitman (Advanced Publishing Program), Boston, MA, [5] J. Hadamard. Bull. Univ. Princeton, 13:49, [6] J. Hadamard. Lectures on the Cauchy Problem in Linear Partial Differential Equations. Yale University Press, New Haven, [7] Ian Knowles. Uniqueness for an elliptic inverse problem. SIAM J. Appl. Math., 59(4): , Available online at [8] Ian Knowles, Tuan A. Le, and Aimin Yan. On the recovery of multiple flow parameters from transient head data. J. Comp. Appl. Math., 169:1 15, [9] Ian Knowles, Michael Teubner, Aimin Yan, Paul Rasser, and Jong Wook Lee. Inverse groundwater modelling in the Willunga Basin, South Australia. Hydrogeology Journal, 15: , [10] Ian Knowles and Robert Wallace. A variational method for numerical differentiation. Numerische Mathematik, 70:91 110, [11] Carlo Miranda. Partial Differential Equations of Elliptic Type. Springer-Verlag, Berlin, second edition, [12] Michael Struwe. Variational methods, volume 34 of Ergebnisse der Mathematik und ihrer Grenzgebiete (3) [Results in Mathematics and Related Areas (3)]. Springer-Verlag, Berlin, second edition, Applications to nonlinear partial differential equations and Hamiltonian systems. [13] Joachim Weidmann. Linear operators in Hilbert spaces. Springer-Verlag, New York, [14] Aimin Yan. An Inverse Groundwater Model. PhD thesis, University of Alabama at Birmingham, 2004.
CONDITIONAL WELL-POSEDNESS FOR AN ELLIPTIC INVERSE PROBLEM
CONDITIONAL WELL-POSEDNESS FOR AN ELLIPTIC INVERSE PROBLEM IAN KNOWLES AND MARY A. LARUSSA Abstract. For the elliptic equation (p(x) v) + λq(x)v = f, x R n, the problem of determining when one or more
More informationFact Sheet Functional Analysis
Fact Sheet Functional Analysis Literature: Hackbusch, W.: Theorie und Numerik elliptischer Differentialgleichungen. Teubner, 986. Knabner, P., Angermann, L.: Numerik partieller Differentialgleichungen.
More informationThe Dirichlet s P rinciple. In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation:
Oct. 1 The Dirichlet s P rinciple In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation: 1. Dirichlet s Principle. u = in, u = g on. ( 1 ) If we multiply
More informationCHAPTER VIII HILBERT SPACES
CHAPTER VIII HILBERT SPACES DEFINITION Let X and Y be two complex vector spaces. A map T : X Y is called a conjugate-linear transformation if it is a reallinear transformation from X into Y, and if T (λx)
More informationTheory of PDE Homework 2
Theory of PDE Homework 2 Adrienne Sands April 18, 2017 In the following exercises we assume the coefficients of the various PDE are smooth and satisfy the uniform ellipticity condition. R n is always an
More informationEigenvalues and Eigenfunctions of the Laplacian
The Waterloo Mathematics Review 23 Eigenvalues and Eigenfunctions of the Laplacian Mihai Nica University of Waterloo mcnica@uwaterloo.ca Abstract: The problem of determining the eigenvalues and eigenvectors
More informationVariational Formulations
Chapter 2 Variational Formulations In this chapter we will derive a variational (or weak) formulation of the elliptic boundary value problem (1.4). We will discuss all fundamental theoretical results that
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More informationElliptic Partial Differential Equations of Second Order
David Gilbarg Neil S.Trudinger Elliptic Partial Differential Equations of Second Order Reprint of the 1998 Edition Springer Chapter 1. Introduction 1 Part I. Linear Equations Chapter 2. Laplace's Equation
More informationYour first day at work MATH 806 (Fall 2015)
Your first day at work MATH 806 (Fall 2015) 1. Let X be a set (with no particular algebraic structure). A function d : X X R is called a metric on X (and then X is called a metric space) when d satisfies
More informationFunctional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability...
Functional Analysis Franck Sueur 2018-2019 Contents 1 Metric spaces 1 1.1 Definitions........................................ 1 1.2 Completeness...................................... 3 1.3 Compactness......................................
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More information4. Solvability of elliptic PDEs
4. Solvability of elliptic PDEs 4.1 Weak formulation Let us first consider the Dirichlet problem for the Laplacian with homogeneous boundary conditions in a bounded open domain R n with C 1 boundary, u
More informationOn m-accretive Schrödinger operators in L p -spaces on manifolds of bounded geometry
On m-accretive Schrödinger operators in L p -spaces on manifolds of bounded geometry Ognjen Milatovic Department of Mathematics and Statistics University of North Florida Jacksonville, FL 32224 USA. Abstract
More informationChapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems
Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Elliptic boundary value problems often occur as the Euler equations of variational problems the latter
More informationTHEOREMS, ETC., FOR MATH 515
THEOREMS, ETC., FOR MATH 515 Proposition 1 (=comment on page 17). If A is an algebra, then any finite union or finite intersection of sets in A is also in A. Proposition 2 (=Proposition 1.1). For every
More informationEXISTENCE OF WEAK SOLUTIONS FOR A NONUNIFORMLY ELLIPTIC NONLINEAR SYSTEM IN R N. 1. Introduction We study the nonuniformly elliptic, nonlinear system
Electronic Journal of Differential Equations, Vol. 20082008), No. 119, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu login: ftp) EXISTENCE
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationTitle: Localized self-adjointness of Schrödinger-type operators on Riemannian manifolds. Proposed running head: Schrödinger-type operators on
Title: Localized self-adjointness of Schrödinger-type operators on Riemannian manifolds. Proposed running head: Schrödinger-type operators on manifolds. Author: Ognjen Milatovic Department Address: Department
More information2.3 Variational form of boundary value problems
2.3. VARIATIONAL FORM OF BOUNDARY VALUE PROBLEMS 21 2.3 Variational form of boundary value problems Let X be a separable Hilbert space with an inner product (, ) and norm. We identify X with its dual X.
More informationTHE FORM SUM AND THE FRIEDRICHS EXTENSION OF SCHRÖDINGER-TYPE OPERATORS ON RIEMANNIAN MANIFOLDS
THE FORM SUM AND THE FRIEDRICHS EXTENSION OF SCHRÖDINGER-TYPE OPERATORS ON RIEMANNIAN MANIFOLDS OGNJEN MILATOVIC Abstract. We consider H V = M +V, where (M, g) is a Riemannian manifold (not necessarily
More informationCOMBINED EFFECTS FOR A STATIONARY PROBLEM WITH INDEFINITE NONLINEARITIES AND LACK OF COMPACTNESS
Dynamic Systems and Applications 22 (203) 37-384 COMBINED EFFECTS FOR A STATIONARY PROBLEM WITH INDEFINITE NONLINEARITIES AND LACK OF COMPACTNESS VICENŢIU D. RĂDULESCU Simion Stoilow Mathematics Institute
More informationOn uniqueness in the inverse conductivity problem with local data
On uniqueness in the inverse conductivity problem with local data Victor Isakov June 21, 2006 1 Introduction The inverse condictivity problem with many boundary measurements consists of recovery of conductivity
More informationNONLINEAR FREDHOLM ALTERNATIVE FOR THE p-laplacian UNDER NONHOMOGENEOUS NEUMANN BOUNDARY CONDITION
Electronic Journal of Differential Equations, Vol. 2016 (2016), No. 210, pp. 1 7. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu NONLINEAR FREDHOLM ALTERNATIVE FOR THE p-laplacian
More informationON THE EXISTENCE OF TRANSMISSION EIGENVALUES. Andreas Kirsch1
Manuscript submitted to AIMS Journals Volume 3, Number 2, May 29 Website: http://aimsciences.org pp. 1 XX ON THE EXISTENCE OF TRANSMISSION EIGENVALUES Andreas Kirsch1 University of Karlsruhe epartment
More information08a. Operators on Hilbert spaces. 1. Boundedness, continuity, operator norms
(February 24, 2017) 08a. Operators on Hilbert spaces Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ [This document is http://www.math.umn.edu/ garrett/m/real/notes 2016-17/08a-ops
More informationDue Giorni di Algebra Lineare Numerica (2GALN) Febbraio 2016, Como. Iterative regularization in variable exponent Lebesgue spaces
Due Giorni di Algebra Lineare Numerica (2GALN) 16 17 Febbraio 2016, Como Iterative regularization in variable exponent Lebesgue spaces Claudio Estatico 1 Joint work with: Brigida Bonino 1, Fabio Di Benedetto
More informationNUMERICAL RADIUS OF A HOLOMORPHIC MAPPING
Geometric Complex Analysis edited by Junjiro Noguchi et al. World Scientific, Singapore, 1995 pp.1 7 NUMERICAL RADIUS OF A HOLOMORPHIC MAPPING YUN SUNG CHOI Department of Mathematics Pohang University
More informationCompact operators on Banach spaces
Compact operators on Banach spaces Jordan Bell jordan.bell@gmail.com Department of Mathematics, University of Toronto November 12, 2017 1 Introduction In this note I prove several things about compact
More informationContents: 1. Minimization. 2. The theorem of Lions-Stampacchia for variational inequalities. 3. Γ -Convergence. 4. Duality mapping.
Minimization Contents: 1. Minimization. 2. The theorem of Lions-Stampacchia for variational inequalities. 3. Γ -Convergence. 4. Duality mapping. 1 Minimization A Topological Result. Let S be a topological
More informationDIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION
Meshless Methods in Science and Engineering - An International Conference Porto, 22 DIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION Robert Schaback Institut für Numerische und Angewandte Mathematik (NAM)
More informationAn introduction to some aspects of functional analysis
An introduction to some aspects of functional analysis Stephen Semmes Rice University Abstract These informal notes deal with some very basic objects in functional analysis, including norms and seminorms
More information4 Sobolev spaces, trace theorem and normal derivative
4 Sobolev spaces, trace theorem and normal derivative Throughout, n will be a sufficiently smooth, bounded domain. We use the standard Sobolev spaces H 0 ( n ) := L 2 ( n ), H 0 () := L 2 (), H k ( n ),
More informationOn the simplest expression of the perturbed Moore Penrose metric generalized inverse
Annals of the University of Bucharest (mathematical series) 4 (LXII) (2013), 433 446 On the simplest expression of the perturbed Moore Penrose metric generalized inverse Jianbing Cao and Yifeng Xue Communicated
More informationFUNCTIONAL ANALYSIS LECTURE NOTES: WEAK AND WEAK* CONVERGENCE
FUNCTIONAL ANALYSIS LECTURE NOTES: WEAK AND WEAK* CONVERGENCE CHRISTOPHER HEIL 1. Weak and Weak* Convergence of Vectors Definition 1.1. Let X be a normed linear space, and let x n, x X. a. We say that
More informationI teach myself... Hilbert spaces
I teach myself... Hilbert spaces by F.J.Sayas, for MATH 806 November 4, 2015 This document will be growing with the semester. Every in red is for you to justify. Even if we start with the basic definition
More informationand finally, any second order divergence form elliptic operator
Supporting Information: Mathematical proofs Preliminaries Let be an arbitrary bounded open set in R n and let L be any elliptic differential operator associated to a symmetric positive bilinear form B
More informationTHE STOKES SYSTEM R.E. SHOWALTER
THE STOKES SYSTEM R.E. SHOWALTER Contents 1. Stokes System 1 Stokes System 2 2. The Weak Solution of the Stokes System 3 3. The Strong Solution 4 4. The Normal Trace 6 5. The Mixed Problem 7 6. The Navier-Stokes
More informationThe Dirichlet-to-Neumann operator
Lecture 8 The Dirichlet-to-Neumann operator The Dirichlet-to-Neumann operator plays an important role in the theory of inverse problems. In fact, from measurements of electrical currents at the surface
More informationLinear Analysis Lecture 5
Linear Analysis Lecture 5 Inner Products and V Let dim V < with inner product,. Choose a basis B and let v, w V have coordinates in F n given by x 1. x n and y 1. y n, respectively. Let A F n n be the
More informationSELF-ADJOINTNESS OF SCHRÖDINGER-TYPE OPERATORS WITH SINGULAR POTENTIALS ON MANIFOLDS OF BOUNDED GEOMETRY
Electronic Journal of Differential Equations, Vol. 2003(2003), No.??, pp. 1 8. ISSN: 1072-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) SELF-ADJOINTNESS
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationNonlinear Analysis 72 (2010) Contents lists available at ScienceDirect. Nonlinear Analysis. journal homepage:
Nonlinear Analysis 72 (2010) 4298 4303 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na Local C 1 ()-minimizers versus local W 1,p ()-minimizers
More informationTraces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains
Traces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains Sergey E. Mikhailov Brunel University West London, Department of Mathematics, Uxbridge, UB8 3PH, UK J. Math. Analysis
More informationLecture 17. Higher boundary regularity. April 15 th, We extend our results to include the boundary. Let u C 2 (Ω) C 0 ( Ω) be a solution of
Lecture 7 April 5 th, 004 Higher boundary regularity We extend our results to include the boundary. Higher a priori regularity upto the boundary. Let u C () C 0 ( ) be a solution of Lu = f on, u = ϕ on.
More informationMAT 578 FUNCTIONAL ANALYSIS EXERCISES
MAT 578 FUNCTIONAL ANALYSIS EXERCISES JOHN QUIGG Exercise 1. Prove that if A is bounded in a topological vector space, then for every neighborhood V of 0 there exists c > 0 such that tv A for all t > c.
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationFunctional Analysis F3/F4/NVP (2005) Homework assignment 3
Functional Analysis F3/F4/NVP (005 Homework assignment 3 All students should solve the following problems: 1. Section 4.8: Problem 8.. Section 4.9: Problem 4. 3. Let T : l l be the operator defined by
More informationChapter 2 Finite Element Spaces for Linear Saddle Point Problems
Chapter 2 Finite Element Spaces for Linear Saddle Point Problems Remark 2.1. Motivation. This chapter deals with the first difficulty inherent to the incompressible Navier Stokes equations, see Remark
More informationApplied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books.
Applied Analysis APPM 44: Final exam 1:3pm 4:pm, Dec. 14, 29. Closed books. Problem 1: 2p Set I = [, 1]. Prove that there is a continuous function u on I such that 1 ux 1 x sin ut 2 dt = cosx, x I. Define
More informationProof by bootstrapping
Proof by bootstrapping Jordan Bell jordan.bell@gmail.com Department of Mathematics, University of Toronto May 4, 2015 The Oxford English Dictionary defines to bootstrap as the following: To make use of
More informationINF-SUP CONDITION FOR OPERATOR EQUATIONS
INF-SUP CONDITION FOR OPERATOR EQUATIONS LONG CHEN We study the well-posedness of the operator equation (1) T u = f. where T is a linear and bounded operator between two linear vector spaces. We give equivalent
More informationReal Analysis Notes. Thomas Goller
Real Analysis Notes Thomas Goller September 4, 2011 Contents 1 Abstract Measure Spaces 2 1.1 Basic Definitions........................... 2 1.2 Measurable Functions........................ 2 1.3 Integration..............................
More informationAn ill-posed boundary value problem for the Helmholtz equation on Lipschitz domains
c de Gruyter 2009 J. Inv. Ill-Posed Problems 17 (2009), 703 711 DOI 10.1515 / JIIP.2009.041 An ill-posed boundary value problem for the Helmholtz equation on Lipschitz domains W. Arendt and T. Regińska
More informationON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM
ON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM OLEG ZUBELEVICH DEPARTMENT OF MATHEMATICS THE BUDGET AND TREASURY ACADEMY OF THE MINISTRY OF FINANCE OF THE RUSSIAN FEDERATION 7, ZLATOUSTINSKY MALIY PER.,
More informationProblem Set 6: Solutions Math 201A: Fall a n x n,
Problem Set 6: Solutions Math 201A: Fall 2016 Problem 1. Is (x n ) n=0 a Schauder basis of C([0, 1])? No. If f(x) = a n x n, n=0 where the series converges uniformly on [0, 1], then f has a power series
More informationA metric space X is a non-empty set endowed with a metric ρ (x, y):
Chapter 1 Preliminaries References: Troianiello, G.M., 1987, Elliptic differential equations and obstacle problems, Plenum Press, New York. Friedman, A., 1982, Variational principles and free-boundary
More informationPartial Differential Equations and Sobolev Spaces MAT-INF4300 autumn Snorre H. Christiansen November 10, 2016
Partial Differential Equations and Sobolev Spaces MAT-INF4300 autumn 2016 Snorre H. Christiansen November 10, 2016 1 2 Contents 1 Introduction 5 2 Notations for differential operators 7 3 Generalities
More informationTHE RECONSTRUCTION OF GROUNDWATER PARAMETERS FROM HEAD DATA IN AN UNCONFINED AQUIFER
THE RECONSTRUCTION OF GROUNDWATER PARAMETERS FROM HEAD DATA IN AN UNCONFINED AQUIFER IAN KNOWLES AND AIMIN AN Abstract. The estimation of groundwater flow parameters from head measurements and other ancillary
More informationMAT 570 REAL ANALYSIS LECTURE NOTES. Contents. 1. Sets Functions Countability Axiom of choice Equivalence relations 9
MAT 570 REAL ANALYSIS LECTURE NOTES PROFESSOR: JOHN QUIGG SEMESTER: FALL 204 Contents. Sets 2 2. Functions 5 3. Countability 7 4. Axiom of choice 8 5. Equivalence relations 9 6. Real numbers 9 7. Extended
More informationOn Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma. Ben Schweizer 1
On Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma Ben Schweizer 1 January 16, 2017 Abstract: We study connections between four different types of results that
More informationTechnische Universität Graz
Technische Universität Graz Stability of the Laplace single layer boundary integral operator in Sobolev spaces O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2016/2 Technische
More informationTakens embedding theorem for infinite-dimensional dynamical systems
Takens embedding theorem for infinite-dimensional dynamical systems James C. Robinson Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. E-mail: jcr@maths.warwick.ac.uk Abstract. Takens
More informationNotions such as convergent sequence and Cauchy sequence make sense for any metric space. Convergent Sequences are Cauchy
Banach Spaces These notes provide an introduction to Banach spaces, which are complete normed vector spaces. For the purposes of these notes, all vector spaces are assumed to be over the real numbers.
More informationt y n (s) ds. t y(s) ds, x(t) = x(0) +
1 Appendix Definition (Closed Linear Operator) (1) The graph G(T ) of a linear operator T on the domain D(T ) X into Y is the set (x, T x) : x D(T )} in the product space X Y. Then T is closed if its graph
More informationHilbert Spaces. Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space.
Hilbert Spaces Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space. Vector Space. Vector space, ν, over the field of complex numbers,
More informationWentzell Boundary Conditions in the Nonsymmetric Case
Math. Model. Nat. Phenom. Vol. 3, No. 1, 2008, pp. 143-147 Wentzell Boundary Conditions in the Nonsymmetric Case A. Favini a1, G. R. Goldstein b, J. A. Goldstein b and S. Romanelli c a Dipartimento di
More informationAnalysis Preliminary Exam Workshop: Hilbert Spaces
Analysis Preliminary Exam Workshop: Hilbert Spaces 1. Hilbert spaces A Hilbert space H is a complete real or complex inner product space. Consider complex Hilbert spaces for definiteness. If (, ) : H H
More informationMath Tune-Up Louisiana State University August, Lectures on Partial Differential Equations and Hilbert Space
Math Tune-Up Louisiana State University August, 2008 Lectures on Partial Differential Equations and Hilbert Space 1. A linear partial differential equation of physics We begin by considering the simplest
More informationThe Navier-Stokes problem in velocity-pressure formulation :convergence and Optimal Control
The Navier-Stokes problem in velocity-pressure formulation :convergence and Optimal Control A.Younes 1 A. Jarray 2 1 Faculté des Sciences de Tunis, Tunisie. e-mail :younesanis@yahoo.fr 2 Faculté des Sciences
More informationOn Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities
On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
More informationDomain Perturbation for Linear and Semi-Linear Boundary Value Problems
CHAPTER 1 Domain Perturbation for Linear and Semi-Linear Boundary Value Problems Daniel Daners School of Mathematics and Statistics, The University of Sydney, NSW 2006, Australia E-mail: D.Daners@maths.usyd.edu.au
More informationTHE L 2 -HODGE THEORY AND REPRESENTATION ON R n
THE L 2 -HODGE THEORY AND REPRESENTATION ON R n BAISHENG YAN Abstract. We present an elementary L 2 -Hodge theory on whole R n based on the minimization principle of the calculus of variations and some
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationADJOINTS, ABSOLUTE VALUES AND POLAR DECOMPOSITIONS
J. OPERATOR THEORY 44(2000), 243 254 c Copyright by Theta, 2000 ADJOINTS, ABSOLUTE VALUES AND POLAR DECOMPOSITIONS DOUGLAS BRIDGES, FRED RICHMAN and PETER SCHUSTER Communicated by William B. Arveson Abstract.
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More information3 (Due ). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?
MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due ). Show that the open disk x 2 + y 2 < 1 is a countable union of planar elementary sets. Show that the closed disk x 2 + y 2 1 is a countable
More informationOn Ekeland s variational principle
J. Fixed Point Theory Appl. 10 (2011) 191 195 DOI 10.1007/s11784-011-0048-x Published online March 31, 2011 Springer Basel AG 2011 Journal of Fixed Point Theory and Applications On Ekeland s variational
More informationPROBLEMS. (b) (Polarization Identity) Show that in any inner product space
1 Professor Carl Cowen Math 54600 Fall 09 PROBLEMS 1. (Geometry in Inner Product Spaces) (a) (Parallelogram Law) Show that in any inner product space x + y 2 + x y 2 = 2( x 2 + y 2 ). (b) (Polarization
More informationFriedrich symmetric systems
viii CHAPTER 8 Friedrich symmetric systems In this chapter, we describe a theory due to Friedrich [13] for positive symmetric systems, which gives the existence and uniqueness of weak solutions of boundary
More informationSEMILINEAR ELLIPTIC EQUATIONS WITH DEPENDENCE ON THE GRADIENT
Electronic Journal of Differential Equations, Vol. 2012 (2012), No. 139, pp. 1 9. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu SEMILINEAR ELLIPTIC
More informationPartial Differential Equations, 2nd Edition, L.C.Evans The Calculus of Variations
Partial Differential Equations, 2nd Edition, L.C.Evans Chapter 8 The Calculus of Variations Yung-Hsiang Huang 2018.03.25 Notation: denotes a bounded smooth, open subset of R n. All given functions are
More informationWELL POSEDNESS OF PROBLEMS I
Finite Element Method 85 WELL POSEDNESS OF PROBLEMS I Consider the following generic problem Lu = f, where L : X Y, u X, f Y and X, Y are two Banach spaces We say that the above problem is well-posed (according
More informationSobolev spaces. May 18
Sobolev spaces May 18 2015 1 Weak derivatives The purpose of these notes is to give a very basic introduction to Sobolev spaces. More extensive treatments can e.g. be found in the classical references
More informationl(y j ) = 0 for all y j (1)
Problem 1. The closed linear span of a subset {y j } of a normed vector space is defined as the intersection of all closed subspaces containing all y j and thus the smallest such subspace. 1 Show that
More informationChapter 1. Introduction
Chapter 1 Introduction Functional analysis can be seen as a natural extension of the real analysis to more general spaces. As an example we can think at the Heine - Borel theorem (closed and bounded is
More informationOn the distributional divergence of vector fields vanishing at infinity
Proceedings of the Royal Society of Edinburgh, 141A, 65 76, 2011 On the distributional divergence of vector fields vanishing at infinity Thierry De Pauw Institut de Recherches en Mathématiques et Physique,
More informationEquations paraboliques: comportement qualitatif
Université de Metz Master 2 Recherche de Mathématiques 2ème semestre Equations paraboliques: comportement qualitatif par Ralph Chill Laboratoire de Mathématiques et Applications de Metz Année 25/6 1 Contents
More informationMTH 503: Functional Analysis
MTH 53: Functional Analysis Semester 1, 215-216 Dr. Prahlad Vaidyanathan Contents I. Normed Linear Spaces 4 1. Review of Linear Algebra........................... 4 2. Definition and Examples...........................
More informationSPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS
SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties
More informationSTOKES PROBLEM WITH SEVERAL TYPES OF BOUNDARY CONDITIONS IN AN EXTERIOR DOMAIN
Electronic Journal of Differential Equations, Vol. 2013 2013, No. 196, pp. 1 28. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu STOKES PROBLEM
More informationA Posteriori Existence in Adaptive Computations
Report no. 06/11 A Posteriori Existence in Adaptive Computations Christoph Ortner This short note demonstrates that it is not necessary to assume the existence of exact solutions in an a posteriori error
More informationIntroduction to Spectral Theory
P.D. Hislop I.M. Sigal Introduction to Spectral Theory With Applications to Schrodinger Operators Springer Introduction and Overview 1 1 The Spectrum of Linear Operators and Hilbert Spaces 9 1.1 TheSpectrum
More informationReal Analysis, 2nd Edition, G.B.Folland Elements of Functional Analysis
Real Analysis, 2nd Edition, G.B.Folland Chapter 5 Elements of Functional Analysis Yung-Hsiang Huang 5.1 Normed Vector Spaces 1. Note for any x, y X and a, b K, x+y x + y and by ax b y x + b a x. 2. It
More information3. (a) What is a simple function? What is an integrable function? How is f dµ defined? Define it first
Math 632/6321: Theory of Functions of a Real Variable Sample Preinary Exam Questions 1. Let (, M, µ) be a measure space. (a) Prove that if µ() < and if 1 p < q
More informationCONTINUITY OF CP-SEMIGROUPS IN THE POINT-STRONG OPERATOR TOPOLOGY
J. OPERATOR THEORY 64:1(21), 149 154 Copyright by THETA, 21 CONTINUITY OF CP-SEMIGROUPS IN THE POINT-STRONG OPERATOR TOPOLOGY DANIEL MARKIEWICZ and ORR MOSHE SHALIT Communicated by William Arveson ABSTRACT.
More informationSemi-strongly asymptotically non-expansive mappings and their applications on xed point theory
Hacettepe Journal of Mathematics and Statistics Volume 46 (4) (2017), 613 620 Semi-strongly asymptotically non-expansive mappings and their applications on xed point theory Chris Lennard and Veysel Nezir
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationA Brief Introduction to Functional Analysis
A Brief Introduction to Functional Analysis Sungwook Lee Department of Mathematics University of Southern Mississippi sunglee@usm.edu July 5, 2007 Definition 1. An algebra A is a vector space over C with
More informationParameter Identification
Lecture Notes Parameter Identification Winter School Inverse Problems 25 Martin Burger 1 Contents 1 Introduction 3 2 Examples of Parameter Identification Problems 5 2.1 Differentiation of Data...............................
More information