arxiv: v1 [math.ap] 23 Jan 2019

Transcription

1 A PATIAL DATA INVESE POBLEM FO THE CONVECTION-DIFFUSION EUATION SUMAN KUMA SAHOO AND MANMOHAN VASHISTH arxiv: v1 [math.ap] 23 Jan 2019 Abstract. In this article, we study the unique determination of convection term and the time dependent density coefficient appearing in a convection-diffusion equation from partial Dirichlet to Neumann map measured on boundary. Keywords: Inverse problems, parabolic equation, Carleman estimates, partial boundary data. Mathematics Subject Classifications (2010): 3530, 35K Introduction Let Ω n, n 3 be a bounded simply connected open set with C 2 boundary. For T > 0, let := (0,T) Ω and denote its lateral boundary by Σ := (0,T) Ω. Consider the following initial boundary value problem ( t ) n j=1 ( j +A j (x)) 2 +q(t,x) u(t,x) = 0, (t,x) u(0,x) = 0, x Ω u(t,x) = f(t,x), (t,x) Σ. We assume that A j W 3, (Ω) for all 1 j n and q W 1, (0,T;L (Ω)). Let us denote by and by A(x) := (A 1 (x),a 2 (x),,a n (x)) L A,q := t n ( j +A j (x)) 2 +q(t,x). j=1 Before going to the main context of the article, let us briefly mention the well-posedness of the forward problem. Following [21], for r,s we define the time-dependent Sobolev spaces by H r,s () := L 2 (0,T;H r (Ω)) H s (0,T;L 2 (Ω)); H r,s (Σ) := L 2 (0,T;H r ( Ω)) H s (0,T;L 2 ( Ω)); and H r,s 0 (Σ) := {f Hr,s (Σ) : f(0,x) = 0, x Σ}. (1.1) Now as shown in ([21], Theorem 1.43) that for f H 3 2,3 4 0 (Σ), there exists a unique solution u H 2,1 () of (1.1) with ν u H 1 2,1 4(Σ). Using this we define the Dirichlet to Neumann map Λ A,q (DN in short) by Λ A,q (f) = ν u Σ (1.2) and then it follows from well-posedness (see [21]) that Λ A,q : H 3 2,1 4 0 (Σ) H 1 2,3 4(Σ) is continuous. In this article, we are interested in recovery of A(x) and q(t,x) in (1.1) from the information of the DN map Λ A,q measured on a subset of Σ. It is well-known [see [51]] that one cannot determine the coefficient A uniquely from the DN map Λ A,q measured on lateral boundary Σ because there is a gauge invariance associated with A(x). So one can only hope to recover A uniquely upto a potential term however the coefficient q(t, x) can be determined uniquely. 1

2 2 SAHOO AND VASHISTH The initial boundary value problem (1.1) is known as a convection-diffusion equation with constant diffusion. The coefficients A and q are called convection term and density coefficient respectively. The convection-diffusion equations appear in chemical engineering, heat transfer and probabilistic study of diffusion process etc. Determination of the coefficients from boundary measurements appearing in parabolic partial differential equations have been studied by several authors. Isakov in [31] considered the problem of determining the time-independent coefficient q for the case when A = 0 in (1.1), from the DN map and he proved the uniqueness result by showing the density of the product of solutions (inspired by the work of [52]). Avdonin and Seidman in [1] studied the problem of determining the time-independent density coefficient q(x) appearing in (1.1) by using the boundary control method pioneered by Belishev, Kurylev, Lassas and others see [2, 4, 35] and references therein. In [21], Choulli proved the stability estimate analogous to the uniqueness problem considered in [31]. In [23] problem of determining the first order coefficients appearing in a parabolic equations in one dimension from the data measured at final time is studied. Cheng and Yamamoto in [14] proved the unique determination of convection term A(x) (when q = 0 in (1.1)) from a single boundary measurement in two dimension. Gaitan and Kian [28] using the global Carleman estimate used for hyperbolic equations [see [10]], proved the stable determination of the time-dependent coefficient q(t, x) in a bounded waveguide. Choulli and Kian in [19] proved the stability estimate for determining the time-dependent coefficient q from partial DN map. For more works related to parabolic inverse problems, we refer to [6, 16, 17, 18, 20, 21, 23, 28, 32, 33, 45] and the references therein. We also mention the work of [5, 7, 9, 41] related to dynamical Schrödinger equation and for steady state convection-diffusion equation we refer to [13, 15]. ecently Caro and Kian in [12] established the unique determination of convection term together with non-linearity appearing in the equation from the knowledge of DN map measured on Σ. Inspired by the work of [19], we consider the problem of determining the full first order space derivative perturbation of heat operator from the partial DN map. We have proved our uniqueness result by using the geometric optics solutions constructed using a Carleman estimate in a Sobolev space of negative order and inverting the ray transform of a vector field which is known only in a very small neighbourhood of the fix direction ω 0. For elliptic and hyperbolic inverse problems these kind of techniques have been used by several authors. elated to our work, we refer to [11, 24] for the elliptic case and to [30, 34, 36, 37, 38, 42] for the hyperbolic case. The article is organized as follows. In 2 we give the statement of the main result. 3 contains the boundary Carleman estimate. In 4 we construct the geometric optics solutions using a Carleman estimate in a Sobolev space of negative order. In 5 we derive an integral identity and 6 contains the proof of the main Theorem Statement of the main result We begin this section by fixing some notation which will be used to state the main result of this article. Following [11], fix an ω 0 S n 1, and define the ω 0 -shadowed and ω 0 -illuminated faces by Ω +,ω0 := {x Ω : ν(x) ω 0 0}, Ω,ω0 := {x Ω : ν(x) ω 0 0} of Ω, where ν(x) is outward unit normal to Ω at x Ω. Corresponding to Ω ±,ω0, we denote the lateral boundary parts by Σ ±,ω0 := (0,T) Ω ±,ω0. We denote by F = (0,T) F and G = (0,T) G where F and G are small enough open neighbourhoods of Ω +,ω0 and Ω,ω0 respectively in Ω. Since Ω is bounded so we can choose a smallest > 0 such that Ω B(0,), where B(0,) n is a ball of radius with center at origin. Now we define the admissible set A of vector fields A(x) appearing in (1.1) by A := { A W 3, (Ω) : A 1 9 We prove the uniqueness result for A A. More precisely, we prove the following: }. (2.1)

3 PATIAL DATA INVESE POBLEM 3 Theorem 2.1. Let ( A (1),q 1 ) and ( A (2),q 2 ) be two sets of coefficients such that A (i) A and q i W 1, (0,T;L (Ω)) for i = 1,2. Let u i be the solutions to (1.1) when (A,q) = ( A (i),q i ) and ΛA (i),q i for i = 1,2 be the DN map defined by (1.2) corresponding to u i. Now if then there exists a function Φ W 4, (Ω) such that and provided A (1) (x) = A (2) (x) for all x Ω. Λ A (1),q 1 (f) G = Λ A (2),q 2 (f) G, for all f H 2,1 0 (Σ), A (1) (x) A (2) (x) = x Φ(x); for all x Ω q 1 (t,x) = q 2 (t,x), (t,x) emark 2.2. The assumption 2.1 on the convection term A is required to prove the boundary and interior Carleman estimates for the operator L A,q. In Theorem 2.1, if in addition one assumes that x A 1 (x) = x A 2 (x) in Ω, then we can show that A 1 (x) = A 2 (x). This kind of situations has been considered in [[7, 22]] for determination of vector field term appearing in the dynamical Schrödinger equation and also in [12] for non-linear parabolic equation. With this extra assumption on A is, we can prove the uniqueness for (A i,q i ) even for the case when both A i s and q is are time-dependent. emark 2.3. While finalising the article we became aware of the work [12], where they have proved uniqueness results for the determination of non-linear terms together with convection and density coefficients in a gauge class from the knowledge of DN map measured on whole of Σ. In contrast to [12], we have proved the uniqueness result for the linear parabolic PDE s from the knowledge of partial DN map measured on a suitable subset of Σ. 3. Boundary Carleman estimate In this section, we prove a Carleman estimate involving the boundary terms for the operator L A,q. We will use this estimate to control the boundary terms appearing in the integral identity given by (6.4), where no information is given. Theorem 3.1. Let ϕ(t,x) = λ 2 t+λω x, where ω S n 1 is fixed. Let u C 2 () such that u(0,x) = 0, for x Ω and u(t,x) = 0, for (t,x) Σ. If A A and q W 1, ((0,T;L (Ω)), then there exists C > 0 depending only on Ω,T,q and A such that λ 2 e 2ϕ u(t,x) 2 dxdt+ e 2ϕ x u(t,x) 2 dxdt+ e 2ϕ(T,x) u(t,x) 2 dx Ω +λ e 2ϕ ν u(t,x) 2 ω ν(x) ds x dt C e 2ϕ L A,q u(t,x) 2 dxdt Σ +,ω (3.1) +Cλ e 2ϕ ν u(t,x) 2 ω ν(x) ds x dt Σ,ω holds for λ large. Proof. Let and denote by L ϕ := e ϕ L A,q e ϕ q(t,x) := q(t,x) x A(x) A(x) 2.

4 4 SAHOO AND VASHISTH Then where L ϕ v(t,x) = e ϕ ( t 2A(x) x + q(t,x))(e ϕ v(t,x)) = ( t v 2 x ϕ x + ( t ϕ x ϕ 2 ϕ ) 2A(x) x 2A(x) x φ ) v(t,x) := P 1 v +P 2 v +P 3 v P 1 := + t ϕ ϕ 2 ϕ = P 2 := t 2 x ϕ x = t 2λω x P 3 := 2A(x) x 2A(x) x ϕ+ q(t,x) = 2A(x) x 2λω A(x)+ q(t,x). Now let I := L ϕ v(t,x) 2 dxdt 1 2 where (P 1 +P 2 )v(t,x) 2 dxdt I 1 := 1 2 (P 1 +P 2 )v(t,x) 2 dxdt I 2 := P 3 v(t,x) 2 dxdt. Next we estimate each of I j for j = 1,2. Now I 1 is I 1 = 1 P 1 v(t,x) P 2 v(t,x) 2 dxdt+ 2 2 (3.2) P 3 v(t,x) 2 dxdt := I 1 I 2 (3.3) P 1 v(t,x) P 2 v(t,x)dxdt. We consider each term separately on the right hand side of the above equation. Using integration by parts and the fact that v Σ = 0, we have x v(t,x) 2 dxdt = v v dxdt 1 v 2 + s v(t,x) 2 dxdt 2s 2 holds for any s > 0. Thus, we have P 1 v(t,x) 2 dxdt 2s v(t,x) 2 dxdt s 2 v(t,x) 2 dxdt. (3.4) Using the arguments from [19], we have P 2 v(t,x) 2 dxdt 1+4λ v(t,x) 2 dxdt (3.5) where > 0 is the radius of the smallest ball B(0,) such that Ω B(0,). Now consider 2 P 1 v(t,x) P 2 v(t,x)dxdt = 2 v(t,x) t v(t,x)dxdt+2λ v(t,x)ω x v(t,x)dxdt = Ω Σ v(t,x) 2 dx+λ ω ν(x) ν v(t,x) 2 ds x dt. (3.6)

5 PATIAL DATA INVESE POBLEM 5 Combining Equations (3.4),(3.5) and (3.6), we have I 1 s v(t,x) 2 dxdt s2 v(t,x) 2 dxdt+ 1+4λ Ω v(t,x) 2 dx+ λ 2 Σ ω ν(x) ν v(t,x) 2 ds x dt. Next we estimate I 2. I 2 ( 2A(x) x 2λω A(x)+ q(t,x))v(t,x) 2 dxdt 2 q 2 v(t,x) 2 dxdt+8λ 2 A 2 Using (3.7) and (3.8) in (3.3), we get I s x v(t,x) 2 dxdt s2 v(t,x) 2 dxdt+ 1+4λ λ 2 Σ 8 A 2 ω ν(x) ν v(t,x) 2 ds x dt 2 q 2 x v(t,x) 2 dxdt. ( 1+4λ 2 ) 32 2 s2 2 2 q 2 8λ 2 A 2 + ( s 8 A 2 ) x v(t,x) 2 dxdt+ 1 2 Ω v(t,x) 2 dxdt+8 A 2 v(t,x) 2 dxdt v(t,x) 2 dxdt+ 1 2 v(t,x) 2 dxdt 8λ 2 A 2 v(t,x) 2 dxdt x v(t,x) 2 dx++ λ 2 Σ v(t,x) 2 dxdt. Ω x v(t,x) 2 dx v(t,x) 2 dxdt ω ν(x) ν v(t,x) 2 ds x dt. Now since A 1 9, therefore taking λ large enough and using the Poincaré inequality, we have L ϕ v(t,x) 2 dxdt C (λ 2 v(t,x) 2 dxdt+ x v(t,x) 2 dxdt+ 1 v(t,x) 2 dx 2 +λ Σ ω ν(x) ν v(t,x) 2 ds x dt holds for large λ, where C > 0 depending only on, A and q. Now after substituting v(t,x) = e ϕ(t,x) u(t,x) in (3.9), we get Ω (3.7) (3.8) ) (3.9) λ 2 e 2ϕ u(t,x) 2 dxdt+ e 2ϕ x u(t,x) 2 dxdt+ e 2ϕ(T,x) u(t,x) 2 dx Ω +λ e 2ϕ ν u(t,x) 2 ω ν(x) ds x dt C e 2ϕ L A,q u(t,x) 2 dxdt Σ +,ω +Cλ Σ,ω e 2ϕ ν u(t,x) 2 ω ν(x) ds x dt

6 6 SAHOO AND VASHISTH which is what we wanted to prove. From Equation (3.9), we have λ 2 v 2 L 2 ) + xv 2 L 2 () C L ϕv 2 L 2 () ; for v C c (). (3.10) 4. Construction of geometric optics solutions In this section, we construct the exponentially growing solution to and exponentially decaying solution to where L A,q given by L A,q u(t,x) = 0; (t,x) L A,qu(t,x) = 0; (t,x) L A,q := t n ( j A j ) 2 + q j=1 is a formal L 2 adjoint of the operator L A,q. We construct these solutions by using a Carleman estimate in a Sobolev space of negative order as used in [24] for elliptic case and in [37, 42] for hyperbolic case. Before going further, following [37] we will give some definition and notation, which will be used later. For m, define space L 2 (0,T;Hλ m(n )) by { L 2 (0,T;Hλ m (n )) := u(t, ) S ( n ) : ( λ 2 + ξ 2) } m/2û(t,ξ) L 2 ( n ) with the norm T u 2 L 2 (0,T;Hλ m(n )) := 0 n ( λ 2 + ξ 2) m û(t,ξ) 2 dξdt where S ( n ) denote the space of all tempered distribution on n and û(t,ξ) is the Fourier transform with respect to space variable x n. We define by { (λ D x,λ m u = Fx ξ 2) m/2 Fx u} here F x and Fx 1 denote the Fourier transform and inverse Fourier transform respectively with respect to space variable x n. With this we define the symbol class Sλ m(n ) of order m by { Sλ m (n ) := C λ C ( n n ) : x α β ξ c ( λ(x,ξ) C α,β λ 2 + ξ 2) } m β for all multi-index α,β N n. With these notations and definitions we state the main theorem of this section as follows: Theorem 4.1. Let ϕ,a and q be as in Theorem 3.1. (1) (Exponentially growing solutions) For λ large, we can find v H 1( 0,T;H 1 (Ω) ) satisfying { L A,q v(t,x) = 0; (t,x), v(0,x) = 0, x Ω of the form where for χ C c v g (t,x) = e ϕ (B g (t,x)+ g (t,x,λ)), (4.1) ((0,T)) arbitrary, we have ( B g (t,x) = χ(t)e i(tτ+x ξ) exp 0 ) ω A(x+sω)ds, (4.2)

7 PATIAL DATA INVESE POBLEM 7 and g (t,x,λ) satisfies the following g (0,x,λ) = 0, for x Ω and g H1 (0,T;Hλ 1(n )) C. (4.3) (2) (Exponentially decaying solutions) Let L A,q be as before. Then for λ large, we can find v H 1( 0,T;H 1 (Ω) ) satisfying { L A,qv(t,x) = 0; (t,x), v(t,x) = 0, x Ω of the form where for χ C c v d (t,x) = e ϕ (B d (t,x)+ d (t,x,λ)), (4.4) ((0,T)) arbitrary, we have ( B d (t,x) = χ(t)exp and d (t,x,λ) satisfies the following 0 ) ω A(x+sω)ds, (4.5) d (T,x,λ) = 0, for x Ω and d H1 (0,T;Hλ 1(n )) C. (4.6) The proof of the above theorem is based on a Carleman estimate in a Sobolev space of negative order. To prove the Carleman estimate stated in Proposition 4.2, we follow the arguments similar to the one used in [24, 37, 42] for elliptic and hyperbolic cases. Proposition 4.2. Let ϕ,a and q be as in Theorem 3.1. Then for λ large enough, we have (1) (Interior Carleman estimate for L A,q ) There exists a constant C > 0 independent of λ and v such that Proof. v L 2 (0,T;L 2 ( n )) C L ϕv L2 (0,T;H 1 λ (n )), (4.7) holds for all v C 1 ([0,T];C c (Ω)) satisfying v(t,x) = 0 for x Ω. (2) (Interior Carleman estimate for L A,q ) There exists a constant C > 0 independent of λ and v such that v L 2 (0,T;L 2 ( n )) C L ϕ v L2 (0,T;H 1 λ (n )), (4.8) holds for all v C 1 ([0,T];C c (Ω)) satisfying v(0,x) = 0 for x Ω. (1) (Proof for (4.7)) Let and denote by Then where L ϕ := e ϕ L A,qe ϕ q (t,x) := q(t,x)+ x A(x) A(x) 2. L ϕv = e ϕ ( t +2A(x) x + q (t,x))e ϕ v(t,x) ] = [ t +2A(x) x + q (t,x) 2λω A(x)+2λω x v(t, x) := P 1v +P 2v +P 3v P 1 := P 2 := t +2λω x P 3 := 2A(x) x 2λω A(x)+ q (t,x). (4.9)

8 8 SAHOO AND VASHISTH Now from (3.2), we have P 1 = P 1 and P 2 = P 2. Hence using the arguments similar to Theorem 3.1, we have v(t,x) 2 dxdt+λ 2 v(t,x) 2 dxdt C L ϕ v(t,x) 2 dxdt for some constant C > 0 independent of λ and v. The above estimate can be written in compact form as v L2 (0,T;H 1 λ (n )) C L ϕ v L 2 () (4.10) for some constant C independent of λ and v. Next using the pseudodifferential operators techniques, we shift the index by 1 in the above estimate. Let us denote by Ω a bounded open subset of n such that Ω Ω. Fix w C 1 ([0,T];C c (Ω)) satisfying w(t,x) = 0 and consider the following D x,λ 1 (P 1 +P 2) D x,λ w. Using the composition of pseudodifferential operators [29, Theorem ], we have Using (4.11) and (3.7), we have D x,λ 1 (P 1 +P 2) D x,λ w = (P 1 +P 2)w. (4.11) (P 1 +P 2 ) D x,λ w L2 (0,T;H 1 λ (n )) = D x,λ 1 (P 1 +P 2 ) D x,λ w L 2 (0,T;L 2 ( n )) = (P 1 +P 2)w L 2 (0,T;L 2 ( n )) (4.12) s w L 2 (0,T;L 2 ( n )) +λ w L 2 (0,T;L 2 ( n )) holds for λ large. Now consider ( P3 D x,λ w L2 (0,T;H 1 λ (n )) 2 (λω A(x) D x,λ w L2 (0,T;H 1 λ (n )) + A(x) x D x,λ w L2 (0,T;H 1 λ (n )) + q w L2 (0,T;H 1 λ (n )) Using the boundedness of the coefficients, we have ( P3 D x,λ w L2 (0,T;H 1 λ (n )) 2 λ A w L 2 (0,T;L 2 ( n )) + A w L 2 (0,T;L 2 ( n )) ) + q w L 2 (0,T;L 2 ( n )). Hence using the inequality as used in (3.3), we get L ϕ D,λ w L 2 (0,T;H 1 λ (n )) C w L 2 (0,T;H 1 λ (n )). Now let χ C c ( Ω) such that χ = 1 in Ω 1, where Ω Ω 1 Ω. Fix w = χ D,λ 1 v in the above equation and using the fact ). (1 χ) D,λ 1 v L 2 (0,T;H m λ (n )) C λ 2 v L 2 (0,T;L 2 ( n ))

9 and PATIAL DATA INVESE POBLEM 9 v L 2 (0,T;L 2 ( n )) = D,λ 1 v L 2 (0,T;H 1 λ (n )) w L 2 (0,T;H 1 λ (n )) + (1 χ) D,λ 1 v L 2 (0,T;H 1 λ (n )) w L 2 (0,T;H 1 λ (n )) + C λ 2 v L 2 (0,T;L 2 ( n )) we get L ϕ v L 2 (0,T;H 1 λ (n )) L ϕ D,λ w L 2 (0,T;H 1 λ (n )) C λ 2 v L 2 (0,T;L 2 ( n )) holds for large λ. Thus finally, we have w L 2 (0,T;H 1 λ (n )) C λ 2 v L 2 (0,T;L 2 ( n )) C v L 2 (0,T;L 2 ( n )) v L 2 (0,T;L 2 ( n )) C L ϕ v L 2 (0,T;H 1 λ (n )), holds for v C 1 ([0,T];Cc (Ω)) such that v(t,x) = 0 and λ large. (2) (Proof for (4.8)) follows by exactly the same argument as that for (4.7). Proposition 4.3. Let ϕ, A and q be as in Theorem 3.1. (1) (Existence of solution to L A,q ) For λ > 0 large enough and v H 1 (0,T;L 2 (Ω)), there exists a solution u H 1( 0,T;H 1 λ (Ω)) of { L ϕ u(t,x) = v(t,x); (t,x), u(0,x) = 0; x Ω and it satisfies u H1 (0,T;H 1 λ (Ω)) C v H 1 (0,T;L 2 (Ω)), (4.13) where C > 0 is a constant independent of λ. (2) (Existence of solution to L A,q ) For λ > 0 large enough and v L2 (), there exists a solution u L 2( 0,T;Hλ 1(Ω)) of { L ϕu(t,x) = v(t,x); (t,x), u(t,x) = 0, x Ω and it satisfies where C > 0 is a constant independent of λ. u H1 (0,T;H 1 λ (Ω)) C v H 1 (0,T;L 2 (Ω)), (4.14) Proof. We will give the proof for the existence of the solution to L A,q while for L A,q follows by using the similar arguments. The proof is based on the standard functional analysis arguments. Consider thespaces := {L ϕ u : u C1 ([0,T];Cc (Ω)) andu(t,x) = 0} asasubspaceofl2( 0,T;H 1 λ (n ) ). Define the linear operator T on S by T(L ϕz) = z(t,x)v(t,x)dtdx, for v H 1 (0,T;L 2 (Ω)). Now using the Carleman estimates (4.7), we have T(L ϕz) z L 2 () v L 2 () C v L 2 () L ϕz L2 (0,T;H 1 λ (n ))

10 10 SAHOO AND VASHISTH holds for all z C 1 ([0,T];Cc (Ω)) with z(t,x) = 0. Therefore using Hahn-Banach theorem, we can extend the linear operator T to L 2( 0,T;H 1 λ (n ) ). We denote the extended map as T and it satisfies T C v L 2 (). Since T is bounded linear functional on L 2( 0,T;H 1 λ (n ) ) and L 2( 0,T;H 1 λ (n ) ) is a Hilbert space. Therefore by iesz representation theorem, there exists a unique u L 2 (0,T;Hλ 1(n )) such that T(f) = f,u L 2 (0,T;L 2 ( n )) for all f L 2( 0,T;H 1 λ (n ) ), (4.15) with u L 2 (0,T;Hλ 1(n )) C v L 2 (). Now for z C 1 ([0,T];Cc (Ω)) satisfying z(t,x) = 0, choose f = L ϕz, in the above equation, we have L ϕ u = v. Now differentiating L ϕ u = v with respect to t, we see L ϕ ( t u) = t v ( t q)u L 2 (). Hence using the previous argument, we can have t u L 2 (0,T;Hλ 1 (Ω)), thus finally we get u H 1 (0,T;Hλ 1 (Ω)) and L ϕ u(t,x) = v(t,x); and u H 1 (0,T;H 1 λ (Ω)) C v H 1 (0,T;L 2 (Ω)). Next we will show that u(0,x) = 0, to prove this, we choose f = L ϕz where z C 1 ([0,T];Cc (Ω)) and z(t,x) = 0. Using this choice of f in (4.15), we have L ϕ z(t,x)u(t,x)dxdt = z(t, x)v(t, x)dxdt. Now using integration by parts and the fact that L ϕ u = v, we have u(0,x)z(0,x)dx = 0. Ω The above identity holds for any z C 1 ([0,T];C c (Ω)) satisfying z(t,x) = 0. Therefore, we conclude that u(0,x) = 0 for x Ω. This completes the first part of the Proposition Proof of the Theorem 4.1. Using the expressions v g and B g from (4.1) and (4.2) respectively and we have the equation for g is L A,q v g (t,x) = 0, L ϕ g (t,x,λ) = L A,q B g (t,x); where L A,q B g (t,x) H 1 (). Now using Proposition 4.3, we have there exists g H 1( 0,T;H 1 λ (Ω)) solution to { L ϕ g (t,x,λ) = L A,q B g (t,x); (t,x), g (0,x,λ) = 0; x Ω and it satisfies the following estimate g H1 (0,T;H 1 λ (Ω)) C where C is a constant independent of λ. This completes the construction of the solution for the operator L A,q u = 0 and existence of the solution for L A,qv = 0, follows in a similar way.

11 PATIAL DATA INVESE POBLEM Integral identity In this section, we derive an integral identity which will beused to prove our uniqueness result. We prove this identity using the geometric optics solutions constructed in 4 and solution to an adjoint problem. Let u i be the solutions to the following initial boundary value problems with vector field coefficient A (i) and scalar potential q i for i = 1,2. L A (i),q i u i (t,x) = 0; (t,x) u i (0,x) = 0; x Ω (5.1) u i (t,x) = f(t,x), (t,x) Σ. Let us denote u(t,x) := (u 1 u 2 )(t,x) ( A(x) := A (1) A (2)) (x) := (A 1 (x),,a n (x)) q i (t,x) := x A (i) (x) A (i) (x) 2 +q i (t,x) (5.2) q(t,x) := q 1 (t,x) q 2 (t,x). Then u is the solution to the following initial boundary value problem: L A (1),q 1 u(t,x) = 2A(x) x u 2 (t,x)+ qu 2 (t,x); (t,x) u(0,x) = 0, x Ω u(t,x) = 0, (t,x) Σ. Let v(t,x) of the form given by (4.4) be the solution to following equation (5.3) L A (1),q 1 v(t,x) = 0 in. (5.4) Also let u 2 of the form given by (4.1) be solution to the following equation L A (2),q 2 u 2 (t,x) = 0, in. (5.5) By the well-posedness result from ([21], Theorem 1.43) we have u H 1 () and ν u L 2 (Σ). Now multiplying (5.3) by v(t,x) and integrating over, we have (2A(x) x u 2 (t,x)+ q(t,x)u 2 (t,x))v(t,x)dxdt = L A (1),q 1 u(t,x)v(t,x)dxdt = u(t,x)l A (1),q 1 v(t,x)dxdt ν u(t,x)v(t,x)ds x dt+ Ω u(t, x)v(t, x)dx where in deriving the above identity we have used the following: u Σ = 0, ν u G = 0, u t=0 = 0 and A (1) is compactly supported in. Now using the fact that L v(t,x) = 0 in and v(t,x) = 0 in Ω, we get, A (1),q 1 (2A(x) x u 2 (t,x)+ q(t,x)u 2 (t,x))v(t,x)dxdt = ν u(t,x)v(t,x)ds x dt (5.6) Lemma 5.1. Let u i for i = 1,2 solutions to (5.1) with u 2 of the form (4.1). Let u = u 1 u 2, and v be of the form (4.4). Then ν u(t,x)v(t,x)ds x dt Cλ 1/2 (5.7) for all ω S n 1 such that ω ω 0 ǫ.

12 12 SAHOO AND VASHISTH Proof. Using the expression of v from (4.4), in the right-hand side of (5.6), we have ν u(t,x)v(t,x)ds x dt 2 ν u(t,x)e ϕ(t,x)( B d (t,x)+ d (t,x) ) ds x dt C ( ) 1+ d 2 L 2 (Σ) e 2ϕ(t,x) ν u(t,x) 2 ds x dt 2 C ( ) 1+ d H 1 (Ω) e 2ϕ(t,x) ν u(t,x) 2 ds x dt. Using (4.6), we have ν u(t,x)v(t,x)ds x dt 2 C e 2ϕ(t,x) ν u(t,x) 2 ds x dt. For ε > 0, define Ω +,ε,ω := {x Ω : ν(x) ω > ε} and Σ +,ε,ω := (0,T) Ω +,ε,ω. Now since Σ +,ε,ω for all ω such that ω ω 0 ε and e 2ϕ(t,x) ν u(t,x) 2 ds x dt e 2ϕ(t,x) ν u(t,x) 2 1 ds x dt = εe 2ϕ(t,x) ν u(t,x) 2 ds x dt ε Σ +,ε,ω Σ +,ε,ω 1 ν ϕe 2ϕ(t,x) ν u(t,x) 2 ds x dt, for all ω near ω 0. λε Σ +,ε,ω Now using the Carleman estimate (3.1) and (5.3), we get ν u(t,x)v(t,x)ds x dt 2 Cλ 1 e 2ϕ(t,x) (2A u 2 (t,x)+ qu 2 (t,x)) 2 dxdt. (5.8) Using expression for u 2 from (4.1) and (4.3), we have e 2ϕ(t,x) (2A u 2 (t,x)+ qu 2 (t,x)) 2 dxdt Cλ 2. Hence finally we get ν u(t,x)v(t,x)ds x dt Cλ 1/2, holds for all ω such that ω ω 0 ǫ. This completes the proof of lemma. In this section, we prove the uniqueness results. 6. Proof of theorem 2.1

13 PATIAL DATA INVESE POBLEM ecovery of the convection term A(x). Since from (5.6), we have (2A(x) x u 2 (t,x)+ q(t,x)u 2 (t,x))v(t,x)dxdt = ν u(t,x)v(t,x)ds x dt Now using Equation (5.7), we have (2A(x) x u 2 (t,x)+ q(t,x)u 2 (t,x))v(t,x)dxdt Cλ 1/2. After diving the above equation by λ and taking λ, we have lim 1 (2A(x) x u 2 (t,x)+ q(t,x)u 2 (t,x))v(t,x)dxdt = 0. (6.1) λ λ Next using the expression for u 2 and v from (4.1) and (4.4) respectively, we get ω A(x)B g (t,x)b d (t,x)ddxdt = 0; for all ω S n 1 such that ω ω 0 ǫ. This after using the expressions for B g (t,x) and B d (t,x) from Equations (4.2) and (4.5) respectively, we get ω A(x)χ 2 (t)e iξ x iτt exp ω A(x+sω)ds dxdt = 0; for all ω S n 1 such that ω ω 0 ǫ. 0 (6.2) Since the above identity holds for all χ Cc (0,T), thus we get ω A(x)e iξ x exp ω A(x+sω)ds dxdt = 0 (6.3) n 0 where ξ ω = 0 for all ω with ω ω 0 ǫ. Now decompose n = ω ω and using this in the above equation, we have iξ k ω A(k +τω)exp ω A(k +τω +sω)ds dτdk = 0, for all ω with ω ω 0 ǫ ω e 0 here dk denotes the Lebesgue measure on ω. After substituting τ +s = s, we get iξ k ω A(k +τω)exp ω A(k + sω)d s dτdk = 0, for all ω with ω ω 0 ǫ. (6.4) ω e Now since iξ k ω e = iξ k e ω τ ω A(k +τω)exp τ exp ω A(k + sω)d s dτdk τ iξ k ω A(k + sω)ds dτdk. = e 1 exp τ ω ω A(k +sω)ds dk.

14 14 SAHOO AND VASHISTH Combining this with (6.4), we get ω A(k +sω)ds = 0, for all k ω, with ω ω 0 ǫ. (6.5) To prove the uniqueness for A(x), we will use the following lemma. Lemma 6.1. Let n 3 and F = (F 1,F 2,,F n ) be a real-valued vector field with F j Cc (Ω) for all 1 j n. Suppose IF(x,ω) := ω F(x+sω)ds = 0 for all ω S n 1 near a fixed ω 0 S n 1 and for all x n. Then there exists a Φ Cc (Ω) such that F(x) = x Φ(x). Proof. The proof uses the arguments similar to the one used in [50, 47, 42], for the case of light ray transforms. Let η := (η 1,η 2,,η n ) n be arbitrary vector and denote ω := (ω 1,ω 2,,ω n ) S n 1. Then we have n (η )IF(x,ω) = ω i η j j F i (x+sω)ds. (6.6) i,j=1 Also by fundamental theorem of calculus, we have d (η F)(x+sω)ds = 0. ds This gives us n i,j=1 ω i η j i F j (x+sω)ds = 0. (6.7) Subtracting (6.7) from (6.6), we get n ω i η j h ij (x+sω)ds = 0, for all x n and ω S n 1 near a fixed ω 0 S n 1. i,j=1 where h ij is a n n matrix with entries Now since h ij (x) = ( j F i i F j )(x), for 1 i,j n. n ω i η j ĥ ij (ξ) = i,j=1 n i,j=1 n ω i η j h ij (x)e iξ x dx. Next decomposing n = ω ω and using (6.7), we get n ω i η j ĥ ij (ξ) = 0; for all η n, ξ ω and ω near ω 0. (6.8) i,j=1 Let {e j : 1 j n} be the standard basis for n where e j is defined as follows; e j := (0,0,,0,}{{} 1,0,,0) jth

15 PATIAL DATA INVESE POBLEM 15 and for simplicity we fix ω 0 = e 1. Now let ξ 0 = e 2 be a fixed vector in n. First we show that ĥij(ξ 0 ) = 0 for all 1 i,j n, then later we will prove that ĥij(ξ) = 0 for all 1 i,j n and ξ near ξ 0. Consider a small perturbation ω 0 (a) of vector ω 0 = e 1 by ω 0 (a) := cosae 1 +sinae k where 3 k n. Then we have ω 0 (a) is near ω 0 for a near 0 and ξ 0 ω 0 (a) = 0. Hence using these choices of ω 0 (a) and η = e j in(6.8), we have This gives us cosaĥ1j(ξ 0 )+sinaĥkj(ξ 0 ) = 0; for all 1 j n, 3 k n and anear 0. ĥ 1j (ξ 0 ) = 0, ĥkj(ξ 0 ) = 0; for all 1 j n,and 3 k n. After using the fact that ĥij = ĥji for all 1 i,j n, we get ĥ ij (ξ 0 ) = 0; for all 1 i,j n. Next we show that ĥij(ξ) = 0 for all ξ ω with ω near ω 0. Now for ω near ω 0, we choose ξ ω as follows: ξ 1 = sinφ 1 cosφ 2 ξ 2 = cosφ 1 ξ 3 = sinφ 1 sinφ 2 cosφ 3. ξ n 1 = sinφ 1 sinφ 2 sinφ n 2 cosθ ξ n = sinφ 1 sinφ 2 sinφ n 2 sinθ. Let A be an orthogonal matrix such that Aξ = e 2, where A is given by ξ 1 ξ 2 ξ 3 ξ φ 1 φ 1 φ 1 n φ 1 ξ 1 ξ 2 ξ 3 ξ n a 31 a 32 a 33 a 3n. A = a k1 a k2 a k3 a. (6.9) kn a k+11 a k+11 a k+13 a k+1n a n1 a n2 a n3 a nn Now choose ω = ( ) ξ 1, ξ2, ξ3,, ξn S n 1 φ 1 φ 1 φ 1 φ 1 then ω is near ω 0 = e 1 when φ i s and θ are close to 0. Next choose ω 0 (a) = cosae 1 +sinae l with l 2, then ω 0 (a) is close to e 1 when a is close to zero. Now define ω(a) by cosacosφ 1 cosφ 2 +a l1 sina ω(a) := A T cosasinφ 1 +a l2 sina ω 0 (a) =.. cosacosφ 1 sinφ 2 sinφ n 2 cosθ+a ln sina

16 16 SAHOO AND VASHISTH Then, we have ω(a) is close ω for a near 0 and ω is close to ω 0 when φ i and θ are close to zero. Also ω(a) ξ = A T ω 0 (a) ξ = ω 0 (a) Aξ = ω 0 (a) e 2 = 0 because of the choice of ω 0 (a). Therefore using these choices of ω(a) and choosing η = e j in (6.8), we have cosa ( n i=1 ξ i φ 1 ĥ ij (ξ) ) ( n ) +sina a li ĥ ij (ξ) = 0, for all 1 j n and a near 0. Now since sina and cosa are linearly independent. We have n ξ i ĥ ij (ξ) = 0, for all 1 j n φ 1 i=1 i=1 n a li ĥ ij (ξ) = 0, for all 1 j n and l 2. i=1 (6.10) Above equations can be written as ξ 1 ξ 2 ξ 3 ξ φ 1 φ 1 φ 1 n ĥ 1j (ξ) φ 1 a 31 a 32 a 33 a 3n ĥ 2j (ξ) ĥ 3j (ξ) a k1 a k2 a k3 a kn. = 0. (6.11) a k+1,1 a k+1,2 a k+1,3 a k+1,n a n1 a n2 a n3 a nn. ĥ nj (ξ) Now let us define matrix B and a n-vector h j as follows: ξ 1 ξ 2 ξ 3 ξ φ 1 φ 1 φ 1 n ĥ 1j (ξ) φ 1 a 31 a 32 a 33 a 3n ĥ 2j (ξ) ĥ 3j (ξ) B = a k1 a k2 a k3 a kn and h j (ξ) =. a k+1,1 a k+1,2 a k+1,3 a k+1,n a n1 a n2 a n3 a nn ĥ nj (ξ) Using these, we have equation (6.11), can be written as Bh j (ξ) = 0; for all 1 j n. (6.12) Note that the matrix B is obtained from A by removing the second row and it is (n 1) n matrix. From the definition of A it is clear that rank of A, rank(a) = n, so the rank of B, rank(b) = n 1. i.e. there exists at-least one non-zero minor of order n 1 of the matrix B. Without loss of generality assume B is non-zero minor of order n 1, where B is given by ξ 2 ξ 3 ξ φ 1 φ 1 n φ 1 a 32 a 33 a 3n. B..... = a k2 a k3 a kn. a k+1,2 a k+1,3 a k+1,n a n2 a n3 a nn

17 PATIAL DATA INVESE POBLEM 17 Now using the fact ĥ11(ξ) = 0 in (6.12), we have B h 1(ξ) = 0 (6.13) T where h 1 (ĥ21 = (ξ),ĥ31(ξ),,ĥn1(ξ)). Since B has full rank, therefore h 1 (ξ) = 0. Also using the fact ĥ ij (ξ) = ĥji(ξ) and h 1 (ξ) = 0, we have B h j (ξ) = 0; for 2 j n (6.14) where h j is an (n 1) vector obtained after deleting first entry from h j = (ĥ1j,ĥ2j,,ĥnj) T. Now using (6.13) and (6.14) in (6.12), we get which gives us h j (ξ) = 0, for 1 j n ĥ ij (ξ) = 0, for all 1 i,j n and for ξ near e 2. Therefore by using the Paley-Wiener theorem, we have ĥ ij (ξ) = 0, for all i,j = 0,1,2,,n and for all ξ n. From the Fourier inversion formula, we get h ij (x) = 0 for all x n. Hence finally by using the definition of h ij (x) and the Poincaré lemma, we have there exists a Φ Cc ( n ) such that F(x) = x Φ. This completes the proof of Lemma 6.1. Using the Lemma 6.1 in Equation (6.5), we have there exists Φ W 4, (Ω) such that A(x) = Φ(x). This completes the proof of uniqueness for the convection term A(x) Uniqueness for the density coefficient q(t,x). Since from the Subsection 6.1, we have A (2) (x) A (1) (x) = Φ(x) for some Φ W 4, (Ω). Therefore after replacing the pair (A (1),q 1 ) by (A (3),q 3 ) where A (3) = A (1) + Φ and q 3 = q 1, we get A (3) = A (2). Using this pair of coefficients in (5.6), we get q(t,x)u 2 (t,x)v(t,x)dxdt = ν u(t,x)v(t,x)ds x dt where q(t,x) := q 1 (t,x) q 2 (t,x). Now using the expressions for u 2 and v from (4.1) and (4.4) respectively and taking λ, we get q(t,x)e i(τt+x ξ) dxdt = 0, for all τ and ξ ω, where ω S n 1 is near ω 0. Now since q L () is zero outside therefore by using the Paley-Wiener theorem we have q 1 (t,x) = q 2 (t,x) for all (t,x). This completes the proof of the Theorem 2.1. Acknowledgments MV thanks Ibtissem Ben Aïcha and Gen Nakamura for the discussions on this problem. Both the authors are thankful to Venky Krishnan for the useful discussions.

18 18 SAHOO AND VASHISTH eferences [1] Avdonin, S. and Seidman, T.I. Identication of q(x) in u t = u qu, from boundary observations, SIAM J. Control Optim. 33 (1995) [2] Avdonin S A and Belishev M Dynamical inverse problem for the Schrödinger equation (BC-method) ed N N Uraltseva Proc. St. Petersburg Mathematical Society vol 10 (Transl. from the ussian by Tamara ozhkovskaya) (Providence, I: American Mathematical Society) (Translations. Series 2. American Mathematical Society 214, 114 (2005); translation from Tr. St-Peterbg. Mat. Obshch. 10, 317 (2004)) [3] Baudouin, L. and Puel, J.-P. Uniqueness and stability in an inverse problem for the Schrödinger equation, Inverse Problems 23 (2007), no. 3, [4] Belishev, Mikhail I. ecent progress in the boundary control method. Inverse Problems, 23(5):1-67, [5] Bellassoued, Mourad and Ferreira, David Dos Santos. Stable determination of coefficients in the dynamical anisotropic Schrödinger equation from the Dirichlet-to-Neumann map, Inverse Problems, 26, 2010, no. 12, , 30, [6] Bellassoued, Mourad; Kian, Yavar and Soccorsi, Eric. An inverse stability result for non-compactly supported potentials by one arbitrary lateral Neumann observation. J. Differential Equations260 (2016), no. 10, [7] Bellassoued, Mourad; Kian, Yavar and Soccorsi, Eric. An inverse problem for the magnetic Schrödinger equation in infinite cylindrical domains, Publ. es. Inst. Math. Sci. 54 (2018), no. 4, [8] Bellassoued, Mourad and Yamamoto, Masahiro. Determination of a coefficient in the wave equation with a single measurement. Appl. Anal. 87 (2008), no. 8, [9] Ben Aïcha, Ibtissem. Stability estimate for an inverse problem for the Schrödinger equation in a magnetic field with time-dependent coefficient, J. Math. Phys. 58 (2017), no. 7, , 21 pp. [10] Bukhgeĭm, Alexnder L. and Klibanov, Michael V. Uniqueness in the large of a class of multidimensional inverse problems, Dokl. Akad. Nauk SSS, 260, 1981, no. 2, , [11] Bukhgeĭm, Alexander L. and Uhlmann, Gunther. ecovering a potential from partial Cauchy data. Comm. Partial Differential Equations, 27(3-4): , [12] Caro, Pedro and Kian, Yavar. Determination of convection terms and quasi-linearities appearing in diffusion equations, [13] Cheng, J. and Yamamoto, M. The global uniqueness for determining two convection coefficients from Dirichlet to Neumann map in two dimensions, Inverse Problems, 16 (2000). [14] Cheng, J. and Yamamoto, M. Identification of convection term in a parabolic equation with a single measurement, Nonlinear Analysis, 50 (2002), [15] Cheng, J. and Yamamoto, M. Determination of Two Convection Coefficients from Dirichlet to Neumann Map in the Two-Dimensional Case, SIAM J. Math. Anal., 35 (2004), [16] Choulli, Mourad. An abstract inverse problem, J. Appl. Math. Stoc. Ana. 4 (2) (1991) [17] Choulli, Mourad. An abstract inverse problem and application, J. Math. Anal. Appl. 160 (1) (1991), [18] Choulli, Mourad and Kian, Yavar. Stability of the determination of a time-dependent coefficient in parabolic equations, Math. Control elat. Fields, 3, 2013, no. 2, , [19] Choulli, Mourad and Kian, Yavar. Logarithmic stability in determining the time-dependent zero order coefficient in a parabolic equation from a partial Dirichlet-to-Neumann map. Application to the determination of a nonlinear term, J. Math. Pures Appl. (9) 114 (2018), [20] Choulli, Mourad; Kian, Yavar and Soccorsi, Eric. Stable determination of time-dependent scalar potential from boundary measurements in a periodic quantum waveguide, SIAM J. Math. Anal., 47, 2015, no. 6, , [21] Choulli, Mourad. Une introduction aux problmes inverses elliptiques et paraboliques, Mathmatiques et Applications, Vol. 65, Springer-Verlag, Berlin, [22] Cristofol, Michel and Soccorsi, Eric. Stability estimate in an inverse problem for non-autonomous magnetic Schrödinger equations, Appl. Anal., 90, 2011, no. 10, , [23] Zui-Cha Deng, Jian-Ning Yu and Liu Yang. Identifying the coefficient of first-order in parabolic equation from final measurement data, Mathematics and Computers in Simulation 77 (2008) [24] Dos Santos Ferreira, David; Kenig, Carlos E.; Sjöstrand, Johannes and Uhlmann, Gunther, Determining a magnetic Schrödinger operator from partial Cauchy data, Comm. Math. Phys., 271, 2007, no. 2, , [25] Eskin, G., Inverse problems for the Schrödinger equations with time-dependent electromagnetic potentials and the Aharonov-Bohm effect, J. Math. Phys., 49, 2008, no.2, , 18, [26] Eskin, G., Inverse hyperbolic problems with time-dependent coefficients, Comm. Partial Differential Equations, 32, 2007, no , , [27] Eskin, G., Inverse problems for general second order hyperbolic equations with time-dependent coefficients, Bull. Math. Sci., 7, 2017, no. 2, , [28] Gaitan, Patricia and Kian, Yavar. A stability result for a time-dependent potential in a cylindrical domain, Inverse Problems, 29, 2013, no. 6, , 18, [29] Hörmander, L. The Analysis of linear partial differential operators, Vol III, Springer-Verlag, Berlin, Heidelberg, 1983.

19 PATIAL DATA INVESE POBLEM 19 [30] Hu, Guang-Hui and Kian, Yavar. Determination of singular time-dependent coefficients for wave equations from full and partial data, Inverse Probl. Imaging 12 (2018), no. 3, [31] Isakov, Victor. Completeness of products of solutions and some inverse problems for PDE, Journal of Differential Equations, 92, 1991, no. 2, , [32] Isakov, Victor. On uniqueness in inverse problems for semilinear parabolic equations, Arch. ational Mech. Anal., 124, 1993, no. 1, 1 12, [33] Isakov, Victor. Inverse Problems for Partial Differential Equations, Volume 127, Springer-Verlag, Berlin, Heidelberg, [34] Kian, Yavar. Unique determination of a time-dependent potential for wave equations from partial data, Ann. Inst. H. Poincaré Anal. Non Linéaire, 34, 2017, no.4, , [35] Alexander Katchalov, Yaroslav Kurylev, and Matti Lassas. Inverse boundary spectral problems, volume 123 of Chapman & Hall/CC Monographs and Surveys in Pure and Applied Mathematics. Chapman & Hall/CC, Boca aton, FL, [36] Kian, Yavar. Stability in the determination of a time-dependent coefficient for wave equations from partial data, J. Math. Anal. Appl., 436, 2016, no. 1, , X. [37] Kian, Yavar. ecovery of time-dependent damping coefficients and potentials appearing in wave equations from partial data, SIAM J. Math. Anal., 48, 2016, no. 6, , [38] Kian, Yavar and Oksanen, Lauri. ecovery of time-dependent coefficient on iemanian manifold for hyperbolic equations, [39] Kian, Yavar; Phan, uang Sang and Soccorsi, Eric. A Carleman estimate for infinite cyclindrical quantum domains and the application to inverse problems, Inverse Problems, 30, 2014, no. 5, , 16, [40] Kian, Yavar; Phan, uang Sang and Soccorsi, Eric. Hölder stable determination of a quantum scalar potential in unbounded cylindrical domains, J. Math. Anal. Appl., 426, 2015, no. 1, , X. [41] Kian, Yavar and Soccorsi, Eric. Hölder stably determining the time-dependent electromagnetic potential of the Schrödinger equation, arxiv: [42] Krishnan, Venkateswaran P. and Vashisth, Manmohan. An inverse problem for the relativistic Schrdinger equation with partial boundary data, Applicable Analysis, doi.org/ / [43] Lions, J.-L. and Magenes, E., Problèmes aux limites non homogènes et applications. Vol. 3, Travaux et echerches Mathématiques, No. 20, Dunod, Paris, 1970, xvii+328. [44] Lions, J.-L. and Magenes, E., Problèmes aux limites non homogènes et applications. Vol. 2, Travaux et echerches Mathématiques, No. 18, Dunod, Paris, 1968, xvi+251. [45] Nakamura, Gen and Sasayama, Satoshi. Inverse boundary value problem for the heat equation with discontinuous coefficients. J. Inverse Ill-Posed Probl. 21 (2013), no. 2, [46] Pohjola, V. A. uniqueness result for an inverse problem of the steady state convection-diffusion equation, SIAM J. Math. Anal., 47 (2015), [47] abieniaharatbar, Siamak. Support theorem for the light-ray transform of vector fields on Minkowski spaces. Inverse Problems & Imaging, 12(2): , [48] akesh and Symes, William W. Uniqueness for an inverse problem for the wave equation. Comm. Partial Differential Equations, 13(1):87 96, [49] Stefanov, Plamen D., Uniqueness of the multi-dimensional inverse scattering problem for time dependent potentials, Math. Z., 201, 1989, no. 4, , [50] Stefanov, Plamen D.. Support theorems for the light ray transform on analytic Lorentzian manifolds. Proc. Amer. Math. Soc., 145(3): , [51] Sun, Ziqi. An inverse boundary value problem for the Schrödinger operator with vector potentials, Trans. Amer. Math. Soc, Vol No. 2, (1992), [52] Sylvester, John and Uhlmann, Gunther. A global uniqueness theorem for an inverse boundary value problem. Ann. of Math. (2), 125(1): , TIF Centre for Applicable Mathematics, Bangalore , India. suman@tifrbng.res.in Beijing Computational Science esearch Center, Beijing , China. mvashisth@csrc.ac.cn