Decoupling of modes for the elastic wave equation in media of limited smoothness
|
|
- Gervais Barton
- 5 years ago
- Views:
Transcription
1 Decoupling of modes for the elastic wave equation in media of limited smoothness Valeriy Brytik 1 Maarten V. de Hoop 1 Hart F. Smith Gunther Uhlmann September, 1 1 Center for Computational and Applied Mathemematics, and Geo-Mathematical Imaging Group, Purdue University, West Lafayette, Indiana, USA. Department of Mathematics, University of Washington, Seattle, Washington, USA. Abstract We establish a decoupling result for the P and S waves of linear, isotropic elasticity, in the setting of twicedifferentiable Lamé parameters. Precisely, we show that the P S components of the wave propagation operator are regularizing of order one on L data, by establishing the diagonalization of the elastic system modulo a L -bounded operator. Effecting the diagonalization in the setting of twice-differentiable coefficients depends upon the symbol of the conjugation operator having a particular structure. 1 Introduction We consider the linear, isotropic elastic wave equation on R n, and Cauchy initial value problem on a time interval [ t, t ], t w i Ax, Dw i =, w i t= = f i, 1.1 t w i t= = g i, where w = w 1,..., w n denotes the displacement vector, and Ax, Dw i = n k µxk w i + i w k + i λxk w k. 1. k=1 The elastic wave equation arises in many physical situations, including the propagation of seismic waves through the Earth. The identification and separation of P and S polarized constituents so-called phases underlies many applications of seismic data processing and imaging, which traditionally has been based on the assumption of smoothness of the medium and asymptotic considerations. As an example, we mention the notion of receiver functions, which were introduced and developed by Vinnik 1977 and Langston Throughout this paper, we assume that the Lamé parameter functions µ and λ belong to C 1,1 R n, in that the first derivatives µ and λ are globally bounded, Lipschitz functions on R n. The principal symbol σ Ax, ξ of Ax, D takes the form σ A = µx ξ I + µx + λx ξ ξ = µx + λx ξ σ Π P + µx ξ σ Π S 1.3 1
2 where Π P and Π S are the projection operators corresponding to the Fourier multipliers Π P w = 1 div w, Π S w = I Π P w σ Π P il = ξ iξ l ξ, σ Π S il = ξ δ il ξ i ξ l ξ. The symbols of Π P and Π S are thus, for each x, ξ, projections onto the eigenspaces of Ax, ξ, respectively of dimension 1 and n 1. This yields that Ax, D Π P w = µx + λx Π P w + l.o.t., Ax, D Π S w = µx Π S w + l.o.t., 1.4 where l.o.t. denotes an operator of first order in w on some finite range of Sobolev spaces, as discussed in the next section. We call Π P w and Π S w, respectively, the P and S modes of the elastic wave w. If the Lamé parameters are constant, then the wave propagation operator commutes with Π P and Π S, so that wave evolution preserves modes which are of pure P or pure S form. In the case of nonconstant µ and λ, there will in general be a nonvanishing coupling of the P and S modes in the evolution. The main result of this paper is that the P S terms in the wave propagator are one order smoother than the P P and S S terms. Precisely, let the solution to?? be written in the form wt, = ctf + stg, and decompose ct = c P P t + c P S t + c SP t + c SS t, where c P P t = Π P ct Π P, etc. We similarly decompose st. The operators c P P t and c SS t are bounded on a range of Sobolev spaces H s R n, in particular for s 1, and s P P t, s SS t map H s 1 R n to H s R n for the same range. All maps are strongly continuous in t. We prove that c P S t : L R n H 1 R n, s P S t : H 1 R n H 1 R n. Additionally, we show that c P P t agrees with the propagator for A P P = Π P AΠ P up to a similarly regularizing term; see Corollary?? for a precise statement. Establishing regularization for the P S components of the propagator requires diagonalizing Ax, D to second order; that is, one order beyond the principal symbol consideration of??, since treating a first order derivative of w as a driving force term will produce terms of the same order as w in the solution. This diagonalization is carried out through conjugation by I + K, where K is a compound pseudodifferential operator with finite regularity symbol of order 1. In the case of smooth µ and λ, a similar diagonalization was carried out, in a more general setting, to all orders in Taylor 1975 using the pseudodifferential calculus. He used this diagonalization to study the reflection of singularities for the system of isotropic elasticity, as well as for other systems with constant multiplicities. Stolk and de Hoop have used Taylor s diagonalization in the study of the linearized seismic inverse scattering problem. Dencker 198 analyzed the propagation of polarization in Euclidean space for linear systems of principal type with smooth coefficients, which includes isotropic elasticity. Wang 1998 applied Dencker s propagation of polarization result to the linearized seismic inverse problem; in particular he discusses the issue of polarization in reflections caused by discontinuities of a perturbation of the medium. In the case of C 1,1 coefficients, the construction of K is more delicate, and carrying through the composition calculus to first order depends on the symbol of K having a special form. We believe that this is the first article studying the interaction of P and S waves for the elastic system in non-smooth media, and it remains open under what conditions this can be carried out in the C 1,1 setting for more general diagonalizable
3 systems with constant multiplicity eigenvalues. We anticipate that our result will have applications in the study of the linearized seismic inverse problem in media of limited smoothness. There has been considerable activity in studying unique continuation in the stationary case for non-smooth Lamé parameters, see Lin et. al. 1 and the references therein for the most recent results. Notation. All pseudodifferential operators that arise in this paper will have symbols that are finite sums px, ξ = N j=1 p m jx, ξ, where p m j is homogeneous in ξ of degree m j. We use σ m p to denote the principal symbol p m x, ξ of p, where m is the order of p. We will use L R n and H s R n to denote the L and L -based Sobolev spaces, in the case of complex valued functions as well as n-dimensional vector valued functions, such as the solution w. Throughout the paper, R will denote an operator that is bounded on L R n ; the form of R may change in each occurrence. The Conjugation Operator In this section we construct the conjugation operator K, which diagonalizes Ax, D up to a L -bounded operator. We decompose A = A P P + A P S + A SP + A SS, where A P P = Π P Ax, D Π P, A P S = Π P Ax, D Π S, A SP = Π S Ax, D Π P, A SS = Π S Ax, D Π S. Note that A SP = A P S, while both A SS and A P P are self-adjoint. Furthermore A SP, and hence A P S, is of order 1, in that A P S : H s+1 R n H s R n for s 1..1 To see this, we write A P S w = 1 div A Π S w = 1 [div, A] Π S w,. where we use div Π S. Each component of the composition then takes the form 1 b Π S w, where b denotes a first-order derivative of µ or λ, which being Lipschitz acts as a multiplier on H s for 1 s 1. This establishes?? for such s. If instead one writes A P S w = Π P A div 1 w and commutes one derivative from the right past A, one obtains?? for the range s. Our main conjugation result is the following. Lemma.1 There exists an operator K of order 1, in the sense that such that K : H s R n H s+1 R n for 1 s 1,.3 Π P I + KAx, Dw = A P P I + Kw + R P w,.4 Π S I + KAx, Dw = A SS I + Kw + R S w,.5 where R P and R S are bounded maps on L R n. Additionally, K P P = K SS =, KSP = K P S,.6 where K P P = Π P K Π P, K SP = Π S K Π P, etc. 3
4 To establish??-??, we expand A and K in terms of their P and S components. Terms such as Π P K SP and A P P K SP vanish. By?? and??, the terms K P S A SP and K SP A P S are bounded on L R n, hence are absorbed into R. Using??, the conditions??-?? thus reduce to the following conditions K P S A SS A P P K P S + A P S = R,.7 K SP A P P A SS K SP + A SP = R,.8 meaning precisely that the left hand sides are bounded maps on L R n. Condition?? is equivalent to?? by taking adjoints, so we will establish??. Consider first the case of smooth µ and λ. Since K P S is of order 1, and A P P and A SS of order, the equality?? reduces to saying that the top i.e. first order terms in the symbol expansion for the left hand side cancel out. By??, this implies that σ 1 K P S = σ 1 A P S µx + λx ξ..9 We will also reduce?? to a principal symbol calculation for the case of µ and λ of regularity C 1,1. Care needs to be taken, however, since by?? the symbol of A P S, hence K P S, is only Lipschitz in x, whereas A P P and A SS are of second order. That we can reduce matters to a symbol calculation depends critically on the fact that K P S has a particular structure, which in turn depends on the special form of A P S. We consider the operator A P S more closely. Since div Π S =, by?? we can write A P S w i = = = n j,k=1 n j,k=1 1 i j k µ k Π S w j + j Π S w k 1 i j k µ j Π S w k n i k µ Π S w k k=1 n j,k=1 1 i j j k µ Π S w k In going from the first to the second line we used that k kπ S w k =, together with symmetry of the expression in j, k. The second term in the third row is a bounded operator on L ; hence we may write Motivated by??, we now define A P S w = µ Π S w + R. K P S w = ΦD 1 µ + λ 1 µ Π S w, and set K SP = K P S so that?? holds. Here, Φξ is a smooth cutoff to ξ M, { 1, ξ M, Φξ =, ξ M, with M depending on the C 1,1 norms of µ and λ. Condition?? holds since λ + µ 1 µ is a multiplier on H s for 1 s 1. Furthermore, by taking M large, we will have Kw H s 1 w Hs, 1 s 1..1 Finally, we note that since K P S w is the gradient of a function, it follows that K P S = Π P K P S = K P S Π S. 4
5 where To evaluate the composition A P P K P S, we observe that we may write A P P = µ + λ Π P + n R k k,.11 k=1 R k : H s R n H s R n, 1 s 1. This follows from the fact that the commutator [Π P k, µ] is bounded on H s R n over this range, which is a consequence of the Calderón commutator theorem; see Calderón Since k K P S is bounded on L, it follows that A P P K P S w = µ + λφd µ + λ 1 µ Π S w + Rw = ΦD µ + λ µ + λ 1 µ Π S w + Rw. To evaluate the composition K P S A SS, we first observe that we can write A SS w = Π S µ Π S w + A 1 w, where A 1 : L H 1. To see this, we note that div Π S, so the λ terms in A do not contribute to A SS. For the same reason we can write Π S k µ i Π S w k = ΠS k µ i Π S w k k which maps L H 1. Finally, we write Π S k µ k Π S w i = ΠS µ Π S w i Π S k µ k Π S w i, k and observe that the last term maps L H 1. By??, the composition K P S A 1 is bounded on L, hence we are reduced to showing that ΦD 1 µ + λ 1 µ divµ Π S w = ΦD µ + λ 1 µ µ Π S w + Rw,.1 where we write Π S = div, and understand that acts on each component of a vector, i.e. as a scalar operator. The relation?? is derived as a consequence of the identities k V W = V W + div W V W V.13 V div W = div V div W div W div V + W div V.14 where in the term div W V the div pairs with. Setting V = µ + λ 1 µ, and W = µ Π S w, then?? is equivalent to ΦD 1 V divw V W = Rw. By??-?? we can write the left hand side as ΦD 1 W divv div V div W + divw divv divw V. k The last two terms, ΦD 1 div W div V W V 5
6 are L bounded operators applied to W, hence of the form Rw, since V is Lipschitz. Since div W = µ Π S w, the same holds for ΦD 1 divv div W. The remaining term is where we use the special structure, since it involves two derivatives of V. The key observation is that divv = div µ + λ 1 µ L R n, since µ = div µ, and hence the expression in fact involves only second derivatives of µ and λ. Consequently, w ΦD 1 W divv : L R n H 1 R n. 3 The Parametrix and Decoupling of Modes We will realize the parametrix construction for?? using parametrices for the purely polarized equations, following the procedure introduced in Smith Theorem 3.1 There exist strongly continuous one parameter families of operators c P t : H s R n H s R n, s P t : H s 1 R n H s R n, t c P t : H s R n H s 1 R n, t s P t : H s 1 R n H s 1 R n, for s and t [ t, t ], such that the solution to the Cauchy problem t A P P vp = G P, v P t= = f P, t v P t= = g P, where Π S f P = Π S g P = Π S G P =, and is given by f P H s R n, g P H s 1 R n, G P L 1 [ t, t ], H s 1 R n, v P t, = c P tf P + s P tg P + s P t s G P s, ds. The same result holds if P is exchanged with S, and both results hold for all t >. Proof. Following pages of Smith 1998 we produce families of operators ct and st, respectively bounded and regularizing of degree 1 on H s R n for all s, with c = I, t c =, s =, t s = I, such that for s 3, t µ + λ ct = T t : H s R n H s 1 R n, t µ + λ st = T 1 t : H s 1 R n H s 1 R n. 3.1 We set c P t = Π P ctπ P, s P t = Π P stπ P, and observe that, for s, t A P P cp t = T P, t : H s R n H s 1 R n, t A P P sp t = T P,1 t : H s 1 R n H s 1 R n. 3. To see this, we use the decomposition?? to see that A P P c P t µ + λ c P t : H s R n H s 1 R n, s. As in the discussion following??, the commutator [Π P, µ + λ ] : H s R n H s 1 R n, s, 6
7 and?? follows from??. The exact evolution operators c P and s P are obtained from c P and s P following pages of Smith 1998, using a convergent iteration. For example, where T t is the solution to c P t = c P t + s P t s T s ds, T t + T P,1 t s T s ds = T P, t. We observe here that, since c P t, s P t and T P, t are orthogonal to Π S on both the right and left, it follows that so is c P t. The bounds on t c P t, etc., follow similarly, using the fact that t ct and t st respectively map H s R n H s 1 R n and H s 1 R n H s 1 R n for all s. Suppose that w is the solution to the linear elastic system??. We pose By?? we can then write v = I + Kw, v P = Π P v, v S = Π S v. t v P = Π P I + KAx, Dw = A P P I + Kw + Rw = A P P v P + RI + K 1 v, where at the last stage we use?? to invert I + K on the range of Sobolev spaces we will consider for w. Similar considerations apply to v S, hence we have the system t A P P vp = R P v, v P t= = f P, t v P t= = g P, t A SS vs = R S v, v S t= = f S, 3.3 t v S t= = g S, where R P and R S are bounded on L R n. Here, we have set f P = Π P I + Kf and g P = Π P I + Kg, and similarly for f S, g S. The system?? is equivalent to the integral system, v P t = c P t f P + s P t g P + v S t = c S t f S + s S t g S + s P t sr P vs, ds, s S t sr S vs, ds. We now restrict attention to data f H s R n, g H s 1 R n, with s 1. Then c P t f P + s P t g P C [ t, t ], H s R n, and the same holds with P replaced by S. For s 1, then s P t sr P is a bounded map on H s R n, and similarly with P replaced by S, so the Volterra type system?? can be solved by iteration to yield a solution v = v P + v S C [ t, t ], H s R n. The solution takes the form v = c P t f P + s P t g P + c S t f S + s S t g S + rt f, g, 3.4 7
8 where rt : L R n H 1 R n H 1 R n and t rt : L R n H 1 R n L R n, with both rt and t rt strongly continuous in t. We next write f = I + Kf, g = I + Kg, and w = I + K 1 v. By??, Consequently, I + K 1 I : H s R n H s+1 R n for 1 s 1. I + K 1 c P ti + K c P t : L R n H 1 R n. By similar results for the other terms, we have the following corollary where rt is of a different form than above. Theorem 3. Let wt, x solve the Cauchy problem??, with data f L R n, g H 1 R n. f P, g P = Π P f, g, and f S, g S = Π S f, g. Then w may be written as wt, = c P tf P + s P tg P + c S tf S + s S tg S + rtf, g, where rt is bounded from L R n H 1 R n to H 1 R n, and t rt is bounded from L R n H 1 R n to L R n. Both rt and t rt are strongly continuous in t. In particular, if the initial data for w is a pure P mode, in that f S = g S =, then Π S w = Π S rtf, g C [ t, t ], H 1 R n, provided that f P L R n and g P H 1 R n. This can be viewed as saying that the P S coupling is regularizing of order 1, and the same holds for pure S modes. This decoupling can also be expressed in terms of mapping properties of the components for the solution operator to??. An immediate result of Theorem?? is the following. Corollary 3.3 Let the solution to?? be written in the form wt, = ctf + stg, and decompose c = c P P + c P S + c SP + c SS, s = s P P + s P S + s SP + s SS, where c P P = Π P c Π P, etc. Then with c P, c S, s P, s S as in Theorem??, we have Let c P P t c P t : L R n H 1 R n, c P S t : L R n H 1 R n, s P P t s P t : H 1 R n H 1 R n, s P S t : H 1 R n H 1 R n. All results hold with P and S interchanged, and are strongly continuous in t. Furthermore, in each case the time derivative of the operator is a bounded map from the indicated domain into L R n. For example, t c P P t t c P t : L R n L R n, with strongly continuous dependence on t. The corollary follows from Theorem?? since all operators listed arise as compositions of rt with the projections Π P and Π S. We remark that if the stronger regularity assumption is made, µ, λ C k,1 with k, then the range of Sobolev indexes s on which rt is regularizing of order 1 increases, and Corollary?? thus applies to a 8
9 range of H s R n. Alternatively, by adding successive terms to K, the system?? for v = I + Kw can be decoupled to higher order, that is, with R P and R S of negative order; see Taylor 1975 for the case of smooth µ, λ where R P and R S are smoothing operators. This requires adding lower order terms to A P P and A SS. The order of gain in regularity in Corollary?? cannot be increased, however, since it is limited by the term of order 1 in K. Acknowledgments The research of Brytik, de Hoop, and Uhlmann was supported in part under NSF CMG grant EAR ; Uhlmann was also partly supported by a Chancellor Professorship at UC Berkeley and a Senior Clay Award. The research of Smith was supported under NSF grant DMS References Calderón, A Commutators of singular integral operators. Proc. Nat. Acad. Sci. U.S.A. 53: Dencker, N On the propagation of polarization sets for systems of real principal type. J. Functional Analysis 46: Langston, C.A Structure under Mount Rainier, Washington, inferred from teleseismic body waves. J. Geophys. Res. 84: Lin, C.L., Nakamura, G., Uhlmann, G., Wang, J.N. 1. Quantitative unique continuation for the Lamé system with less regular coefficients. Preprint, arxiv:15.338v1. Smith, H.F A parametrix construction for wave equations with C 1,1 coefficients. Ann. Inst. Fourier, Grenoble 4: Stolk, C.. Microlocal analysis of a seismic linearized inverse problem. Wave Motion 3:67 9. Stolk, C., de Hoop, M.V.. Microlocal analysis of seismic inverse scattering in anisotropic elastic media. Comm. Pure Appl. Math. 55: Stolk, C., de Hoop, M.V. 6. Seismic inverse scattering in the downward continuation approach. Wave Motion 43: Taylor, M Reflection of singularities of solutions to systems of differential equations. Comm. Pure Appl. Math. 8: Wang, J.N On the relation between singularities of coefficients and singularities of reflected waves in the Lamé system. Inverse Problems 14: Vinnik, L Detection of waves converted from P to SV in the mantle. Phys. Earth Planet. Inter. 15:
Liouville-type theorem for the Lamé system with singular coefficients
Liouville-type theorem for the Lamé system with singular coefficients Blair Davey Ching-Lung Lin Jenn-Nan Wang Abstract In this paper, we study a Liouville-type theorem for the Lamé system with rough coefficients
More informationPOINTWISE BOUNDS ON QUASIMODES OF SEMICLASSICAL SCHRÖDINGER OPERATORS IN DIMENSION TWO
POINTWISE BOUNDS ON QUASIMODES OF SEMICLASSICAL SCHRÖDINGER OPERATORS IN DIMENSION TWO HART F. SMITH AND MACIEJ ZWORSKI Abstract. We prove optimal pointwise bounds on quasimodes of semiclassical Schrödinger
More informationRANDOM PROPERTIES BENOIT PAUSADER
RANDOM PROPERTIES BENOIT PAUSADER. Quasilinear problems In general, one consider the following trichotomy for nonlinear PDEs: A semilinear problem is a problem where the highest-order terms appears linearly
More informationWave equation on manifolds and finite speed of propagation
Wave equation on manifolds and finite speed of propagation Ethan Y. Jaffe Let M be a Riemannian manifold (without boundary), and let be the (negative of) the Laplace-Beltrami operator. In this note, we
More informationRecovery of anisotropic metrics from travel times
Purdue University The Lens Rigidity and the Boundary Rigidity Problems Let M be a bounded domain with boundary. Let g be a Riemannian metric on M. Define the scattering relation σ and the length (travel
More informationA REMARK ON AN EQUATION OF WAVE MAPS TYPE WITH VARIABLE COEFFICIENTS
A REMARK ON AN EQUATION OF WAVE MAPS TYPE WITH VARIABLE COEFFICIENTS DAN-ANDREI GEBA Abstract. We obtain a sharp local well-posedness result for an equation of wave maps type with variable coefficients.
More informationOn the local existence for an active scalar equation in critical regularity setting
On the local existence for an active scalar equation in critical regularity setting Walter Rusin Department of Mathematics, Oklahoma State University, Stillwater, OK 7478 Fei Wang Department of Mathematics,
More informationMICROLOCAL ANALYSIS METHODS
MICROLOCAL ANALYSIS METHODS PLAMEN STEFANOV One of the fundamental ideas of classical analysis is a thorough study of functions near a point, i.e., locally. Microlocal analysis, loosely speaking, is analysis
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationCurriculum Vitae for Hart Smith - Updated April, 2016
Curriculum Vitae for Hart Smith - Updated April, 2016 Education Princeton University, 1984 1988. Ph.D. in Mathematics, January, 1989. Thesis supervisor: E. M. Stein. Thesis title: The subelliptic oblique
More informationTopics in Harmonic Analysis Lecture 6: Pseudodifferential calculus and almost orthogonality
Topics in Harmonic Analysis Lecture 6: Pseudodifferential calculus and almost orthogonality Po-Lam Yung The Chinese University of Hong Kong Introduction While multiplier operators are very useful in studying
More informationThe oblique derivative problem for general elliptic systems in Lipschitz domains
M. MITREA The oblique derivative problem for general elliptic systems in Lipschitz domains Let M be a smooth, oriented, connected, compact, boundaryless manifold of real dimension m, and let T M and T
More informationNew Identities for Weak KAM Theory
New Identities for Weak KAM Theory Lawrence C. Evans Department of Mathematics University of California, Berkeley Abstract This paper records for the Hamiltonian H = p + W (x) some old and new identities
More information12x + 18y = 30? ax + by = m
Math 2201, Further Linear Algebra: a practical summary. February, 2009 There are just a few themes that were covered in the course. I. Algebra of integers and polynomials. II. Structure theory of one endomorphism.
More informationMARY ANN HORN be decoupled into three wave equations. Thus, we would hope that results, analogous to those available for the wave equation, would hold
Journal of Mathematical Systems, Estimation, and Control Vol. 8, No. 2, 1998, pp. 1{11 c 1998 Birkhauser-Boston Sharp Trace Regularity for the Solutions of the Equations of Dynamic Elasticity Mary Ann
More informationIntroduction to Seismic Imaging
Introduction to Seismic Imaging Alison Malcolm Department of Earth, Atmospheric and Planetary Sciences MIT August 20, 2010 Outline Introduction Why we image the Earth How data are collected Imaging vs
More informationPart 2 Introduction to Microlocal Analysis
Part 2 Introduction to Microlocal Analysis Birsen Yazıcı & Venky Krishnan Rensselaer Polytechnic Institute Electrical, Computer and Systems Engineering August 2 nd, 2010 Outline PART II Pseudodifferential
More informationLecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems
University of Warwick, EC9A0 Maths for Economists Peter J. Hammond 1 of 45 Lecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems Peter J. Hammond latest revision 2017 September
More informationIncoming and disappearaing solutions of Maxwell s equations. Université Bordeaux I
1 / 27 Incoming and disappearaing solutions of Maxwell s equations Vesselin Petkov (joint work with F. Colombini and J. Rauch) Université Bordeaux I MSRI, September 9th, 2010 Introduction 2 / 27 0. Introduction
More informationCOMP 558 lecture 18 Nov. 15, 2010
Least squares We have seen several least squares problems thus far, and we will see more in the upcoming lectures. For this reason it is good to have a more general picture of these problems and how to
More informationOn stable inversion of the attenuated Radon transform with half data Jan Boman. We shall consider weighted Radon transforms of the form
On stable inversion of the attenuated Radon transform with half data Jan Boman We shall consider weighted Radon transforms of the form R ρ f(l) = f(x)ρ(x, L)ds, L where ρ is a given smooth, positive weight
More informationMath 302 Outcome Statements Winter 2013
Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a
More informationUpon successful completion of MATH 220, the student will be able to:
MATH 220 Matrices Upon successful completion of MATH 220, the student will be able to: 1. Identify a system of linear equations (or linear system) and describe its solution set 2. Write down the coefficient
More informationWe use the overhead arrow to denote a column vector, i.e., a number with a direction. For example, in three-space, we write
1 MATH FACTS 11 Vectors 111 Definition We use the overhead arrow to denote a column vector, ie, a number with a direction For example, in three-space, we write The elements of a vector have a graphical
More informationMATH 240 Spring, Chapter 1: Linear Equations and Matrices
MATH 240 Spring, 2006 Chapter Summaries for Kolman / Hill, Elementary Linear Algebra, 8th Ed. Sections 1.1 1.6, 2.1 2.2, 3.2 3.8, 4.3 4.5, 5.1 5.3, 5.5, 6.1 6.5, 7.1 7.2, 7.4 DEFINITIONS Chapter 1: Linear
More informationDiscreteness of Transmission Eigenvalues via Upper Triangular Compact Operators
Discreteness of Transmission Eigenvalues via Upper Triangular Compact Operators John Sylvester Department of Mathematics University of Washington Seattle, Washington 98195 U.S.A. June 3, 2011 This research
More informationMobile Robotics 1. A Compact Course on Linear Algebra. Giorgio Grisetti
Mobile Robotics 1 A Compact Course on Linear Algebra Giorgio Grisetti SA-1 Vectors Arrays of numbers They represent a point in a n dimensional space 2 Vectors: Scalar Product Scalar-Vector Product Changes
More informationFFTs in Graphics and Vision. The Laplace Operator
FFTs in Graphics and Vision The Laplace Operator 1 Outline Math Stuff Symmetric/Hermitian Matrices Lagrange Multipliers Diagonalizing Symmetric Matrices The Laplacian Operator 2 Linear Operators Definition:
More informationA SHARP STABILITY ESTIMATE IN TENSOR TOMOGRAPHY
A SHARP STABILITY ESTIMATE IN TENSOR TOMOGRAPHY PLAMEN STEFANOV 1. Introduction Let (M, g) be a compact Riemannian manifold with boundary. The geodesic ray transform I of symmetric 2-tensor fields f is
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationMATH 583A REVIEW SESSION #1
MATH 583A REVIEW SESSION #1 BOJAN DURICKOVIC 1. Vector Spaces Very quick review of the basic linear algebra concepts (see any linear algebra textbook): (finite dimensional) vector space (or linear space),
More informationImage amplitudes in reverse time migration/inversion
Image amplitudes in reverse time migration/inversion Chris Stolk 1, Tim Op t Root 2, Maarten de Hoop 3 1 University of Amsterdam 2 University of Twente 3 Purdue University TRIP Seminar, October 6, 2010
More informationPart 2 Introduction to Microlocal Analysis
Part 2 Introduction to Microlocal Analysis Birsen Yazıcı& Venky Krishnan Rensselaer Polytechnic Institute Electrical, Computer and Systems Engineering March 15 th, 2010 Outline PART II Pseudodifferential(ψDOs)
More informationCan There Be a General Theory of Fourier Integral Operators?
Can There Be a General Theory of Fourier Integral Operators? Allan Greenleaf University of Rochester Conference on Inverse Problems in Honor of Gunther Uhlmann UC, Irvine June 21, 2012 How I started working
More informationShort note on compact operators - Monday 24 th March, Sylvester Eriksson-Bique
Short note on compact operators - Monday 24 th March, 2014 Sylvester Eriksson-Bique 1 Introduction In this note I will give a short outline about the structure theory of compact operators. I restrict attention
More informationNOTES FOR CARDIFF LECTURES ON MICROLOCAL ANALYSIS
NOTES FOR CARDIFF LECTURES ON MICROLOCAL ANALYSIS JARED WUNSCH Note that these lectures overlap with Alex s to a degree, to ensure a smooth handoff between lecturers! Our notation is mostly, but not completely,
More informationComplex geometrical optics solutions for Lipschitz conductivities
Rev. Mat. Iberoamericana 19 (2003), 57 72 Complex geometrical optics solutions for Lipschitz conductivities Lassi Päivärinta, Alexander Panchenko and Gunther Uhlmann Abstract We prove the existence of
More informationProceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005
Proceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005 SOME INVERSE SCATTERING PROBLEMS FOR TWO-DIMENSIONAL SCHRÖDINGER
More informationGEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS
Methods in Geochemistry and Geophysics, 36 GEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS Michael S. ZHDANOV University of Utah Salt Lake City UTAH, U.S.A. 2OO2 ELSEVIER Amsterdam - Boston - London
More informationMATH 167: APPLIED LINEAR ALGEBRA Chapter 3
MATH 167: APPLIED LINEAR ALGEBRA Chapter 3 Jesús De Loera, UC Davis February 18, 2012 Orthogonal Vectors and Subspaces (3.1). In real life vector spaces come with additional METRIC properties!! We have
More informationRecent progress on the explicit inversion of geodesic X-ray transforms
Recent progress on the explicit inversion of geodesic X-ray transforms François Monard Department of Mathematics, University of Washington. Geometric Analysis and PDE seminar University of Cambridge, May
More informationLECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI
LECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI Georgian Technical University Tbilisi, GEORGIA 0-0 1. Formulation of the corresponding
More informationWEYL S THEOREM FOR PAIRS OF COMMUTING HYPONORMAL OPERATORS
WEYL S THEOREM FOR PAIRS OF COMMUTING HYPONORMAL OPERATORS SAMEER CHAVAN AND RAÚL CURTO Abstract. Let T be a pair of commuting hyponormal operators satisfying the so-called quasitriangular property dim
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 informationMathematical foundations - linear algebra
Mathematical foundations - linear algebra Andrea Passerini passerini@disi.unitn.it Machine Learning Vector space Definition (over reals) A set X is called a vector space over IR if addition and scalar
More informationGLOBAL WELL-POSEDNESS FOR NONLINEAR NONLOCAL CAUCHY PROBLEMS ARISING IN ELASTICITY
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 55, pp. 1 7. ISSN: 1072-6691. UL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu GLOBAL WELL-POSEDNESS FO NONLINEA NONLOCAL
More informationExtended isochron rays in prestack depth (map) migration
Extended isochron rays in prestack depth (map) migration A.A. Duchkov and M.V. de Hoop Purdue University, 150 N.University st., West Lafayette, IN, 47907 e-mail: aduchkov@purdue.edu (December 15, 2008)
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationR. Courant and D. Hilbert METHODS OF MATHEMATICAL PHYSICS Volume II Partial Differential Equations by R. Courant
R. Courant and D. Hilbert METHODS OF MATHEMATICAL PHYSICS Volume II Partial Differential Equations by R. Courant CONTENTS I. Introductory Remarks S1. General Information about the Variety of Solutions.
More informationON STRICHARTZ ESTIMATES FOR SCHRÖDINGER OPERATORS IN COMPACT MANIFOLDS WITH BOUNDARY. 1. Introduction
ON STRICHARTZ ESTIMATES FOR SCHRÖDINGER OPERATORS IN COMPACT MANIFOLDS WITH BOUNDARY MATTHEW D. BLAIR, HART F. SMITH, AND CHRISTOPHER D. SOGGE 1. Introduction Let (M, g) be a Riemannian manifold of dimension
More informationShear Stresses and Displacement for Strike-slip Dislocation in an Orthotropic Elastic Half-space with Rigid Surface
International Journal of Applied Science-Research and Review (IJAS) www.ijas.org.uk Original Article Shear Stresses and Displacement for Strike-slip Dislocation in an Orthotropic Elastic Half-space with
More information4.2. ORTHOGONALITY 161
4.2. ORTHOGONALITY 161 Definition 4.2.9 An affine space (E, E ) is a Euclidean affine space iff its underlying vector space E is a Euclidean vector space. Given any two points a, b E, we define the distance
More informationLecture 2: Some basic principles of the b-calculus
Lecture 2: Some basic principles of the b-calculus Daniel Grieser (Carl von Ossietzky Universität Oldenburg) September 20, 2012 Summer School Singular Analysis Daniel Grieser (Oldenburg) Lecture 2: Some
More informationCUTOFF RESOLVENT ESTIMATES AND THE SEMILINEAR SCHRÖDINGER EQUATION
CUTOFF RESOLVENT ESTIMATES AND THE SEMILINEAR SCHRÖDINGER EQUATION HANS CHRISTIANSON Abstract. This paper shows how abstract resolvent estimates imply local smoothing for solutions to the Schrödinger equation.
More informationMath 350 Fall 2011 Notes about inner product spaces. In this notes we state and prove some important properties of inner product spaces.
Math 350 Fall 2011 Notes about inner product spaces In this notes we state and prove some important properties of inner product spaces. First, recall the dot product on R n : if x, y R n, say x = (x 1,...,
More informationMath Linear Algebra II. 1. Inner Products and Norms
Math 342 - Linear Algebra II Notes 1. Inner Products and Norms One knows from a basic introduction to vectors in R n Math 254 at OSU) that the length of a vector x = x 1 x 2... x n ) T R n, denoted x,
More informationSPECTRAL THEORY EVAN JENKINS
SPECTRAL THEORY EVAN JENKINS Abstract. These are notes from two lectures given in MATH 27200, Basic Functional Analysis, at the University of Chicago in March 2010. The proof of the spectral theorem for
More informationAsymptotic wave-splitting in anisotropic linear acoustics
arxiv:0908.1507v1 [math-ph] 11 Aug 2009 Asymptotic wave-splitting in anisotropic linear acoustics Abstract B. L. G. Jonsson and M. Norgren Electromagnetic Engineering, Royal Institute of Technology, Stockholm,
More informationFractional order operators on bounded domains
on bounded domains Gerd Grubb Copenhagen University Geometry seminar, Stanford University September 23, 2015 1. Fractional-order pseudodifferential operators A prominent example of a fractional-order pseudodifferential
More informationClassical differential geometry of two-dimensional surfaces
Classical differential geometry of two-dimensional surfaces 1 Basic definitions This section gives an overview of the basic notions of differential geometry for twodimensional surfaces. It follows mainly
More informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
More informationExistence of minimizers for the pure displacement problem in nonlinear elasticity
Existence of minimizers for the pure displacement problem in nonlinear elasticity Cristinel Mardare Université Pierre et Marie Curie - Paris 6, Laboratoire Jacques-Louis Lions, Paris, F-75005 France Abstract
More informationThe X-ray transform for a non-abelian connection in two dimensions
The X-ray transform for a non-abelian connection in two dimensions David Finch Department of Mathematics Oregon State University Corvallis, OR, 97331, USA Gunther Uhlmann Department of Mathematics University
More informationEQUIVARIANT AND FRACTIONAL INDEX OF PROJECTIVE ELLIPTIC OPERATORS. V. Mathai, R.B. Melrose & I.M. Singer. Abstract
j. differential geometry x78 (2008) 465-473 EQUIVARIANT AND FRACTIONAL INDEX OF PROJECTIVE ELLIPTIC OPERATORS V. Mathai, R.B. Melrose & I.M. Singer Abstract In this note the fractional analytic index,
More informationSeismic inverse scattering by reverse time migration
Seismic inverse scattering by reverse time migration Chris Stolk 1, Tim Op t Root 2, Maarten de Hoop 3 1 University of Amsterdam 2 University of Twente 3 Purdue University MSRI Inverse Problems Seminar,
More informationFIXED POINT ITERATIONS
FIXED POINT ITERATIONS MARKUS GRASMAIR 1. Fixed Point Iteration for Non-linear Equations Our goal is the solution of an equation (1) F (x) = 0, where F : R n R n is a continuous vector valued mapping in
More informationMath 108b: Notes on the Spectral Theorem
Math 108b: Notes on the Spectral Theorem From section 6.3, we know that every linear operator T on a finite dimensional inner product space V has an adjoint. (T is defined as the unique linear operator
More informationParaproducts and the bilinear Calderón-Zygmund theory
Paraproducts and the bilinear Calderón-Zygmund theory Diego Maldonado Department of Mathematics Kansas State University Manhattan, KS 66506 12th New Mexico Analysis Seminar April 23-25, 2009 Outline of
More informationNONLOCAL DIFFUSION EQUATIONS
NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,
More informationVelocity averaging a general framework
Outline Velocity averaging a general framework Martin Lazar BCAM ERC-NUMERIWAVES Seminar May 15, 2013 Joint work with D. Mitrović, University of Montenegro, Montenegro Outline Outline 1 2 L p, p >= 2 setting
More information14 Singular Value Decomposition
14 Singular Value Decomposition For any high-dimensional data analysis, one s first thought should often be: can I use an SVD? The singular value decomposition is an invaluable analysis tool for dealing
More informationASYMPTOTIC BEHAVIOR OF GENERALIZED EIGENFUNCTIONS IN N-BODY SCATTERING
ASYMPTOTIC BEHAVIOR OF GENERALIZED EIGENFUNCTIONS IN N-BODY SCATTERING ANDRAS VASY Abstract. In this paper an asymptotic expansion is proved for locally (at infinity) outgoing functions on asymptotically
More informationMathematical foundations - linear algebra
Mathematical foundations - linear algebra Andrea Passerini passerini@disi.unitn.it Machine Learning Vector space Definition (over reals) A set X is called a vector space over IR if addition and scalar
More informationA LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE. 1.
A LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE THOMAS CHEN AND NATAŠA PAVLOVIĆ Abstract. We prove a Beale-Kato-Majda criterion
More informationTHE KADISON-SINGER PROBLEM AND THE UNCERTAINTY PRINCIPLE
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 0002-9939(XX)0000-0 THE KADISON-SINGER PROBLEM AND THE UNCERTAINTY PRINCIPLE PETER G. CASAZZA AND ERIC WEBER Abstract.
More informationWork of Lars Hörmander. Michael Taylor
Work of Lars Hörmander Michael Taylor Lars Hörmander was one of the giants of twentieth century mathematics. He made a big splash with his Ph.D. thesis, much of which was published in a 1955 Acta Mathematica
More informationSharp estimates for a class of hyperbolic pseudo-differential equations
Results in Math., 41 (2002), 361-368. Sharp estimates for a class of hyperbolic pseudo-differential equations Michael Ruzhansky Abstract In this paper we consider the Cauchy problem for a class of hyperbolic
More informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
More informationReview problems for MA 54, Fall 2004.
Review problems for MA 54, Fall 2004. Below are the review problems for the final. They are mostly homework problems, or very similar. If you are comfortable doing these problems, you should be fine on
More informationChapter 3. Determinants and Eigenvalues
Chapter 3. Determinants and Eigenvalues 3.1. Determinants With each square matrix we can associate a real number called the determinant of the matrix. Determinants have important applications to the theory
More informationON A CLASS OF NONSMOOTH COMPOSITE FUNCTIONS
MATHEMATICS OF OPERATIONS RESEARCH Vol. 28, No. 4, November 2003, pp. 677 692 Printed in U.S.A. ON A CLASS OF NONSMOOTH COMPOSITE FUNCTIONS ALEXANDER SHAPIRO We discuss in this paper a class of nonsmooth
More informationMicrolocal Analysis : a short introduction
Microlocal Analysis : a short introduction Plamen Stefanov Purdue University Mini Course, Fields Institute, 2012 Plamen Stefanov (Purdue University ) Microlocal Analysis : a short introduction 1 / 25 Introduction
More informationA new class of pseudodifferential operators with mixed homogenities
A new class of pseudodifferential operators with mixed homogenities Po-Lam Yung University of Oxford Jan 20, 2014 Introduction Given a smooth distribution of hyperplanes on R N (or more generally on a
More informationNovember 18, 2013 ANALYTIC FUNCTIONAL CALCULUS
November 8, 203 ANALYTIC FUNCTIONAL CALCULUS RODICA D. COSTIN Contents. The spectral projection theorem. Functional calculus 2.. The spectral projection theorem for self-adjoint matrices 2.2. The spectral
More informationOn a class of pseudodifferential operators with mixed homogeneities
On a class of pseudodifferential operators with mixed homogeneities Po-Lam Yung University of Oxford July 25, 2014 Introduction Joint work with E. Stein (and an outgrowth of work of Nagel-Ricci-Stein-Wainger,
More informationRepresentation theory and quantum mechanics tutorial Spin and the hydrogen atom
Representation theory and quantum mechanics tutorial Spin and the hydrogen atom Justin Campbell August 3, 2017 1 Representations of SU 2 and SO 3 (R) 1.1 The following observation is long overdue. Proposition
More informationRandom walk on groups
Random walk on groups Bob Hough September 22, 2015 Bob Hough Random walk on groups September 22, 2015 1 / 12 Set-up G a topological (finite) group, with Borel sigma algebra B P(G) the set of Borel probability
More informationTechnische Universität Dresden
Technische Universität Dresden Als Typoskript gedruckt Herausgeber: Der Rektor arxiv:1502.05496v1 [math.fa] 19 Feb 2015 On a class of block operator matrices in system theory. Sascha Trostorff Institut
More informationANSWERS (5 points) Let A be a 2 2 matrix such that A =. Compute A. 2
MATH 7- Final Exam Sample Problems Spring 7 ANSWERS ) ) ). 5 points) Let A be a matrix such that A =. Compute A. ) A = A ) = ) = ). 5 points) State ) the definition of norm, ) the Cauchy-Schwartz inequality
More informationarxiv: v3 [math.ap] 1 Sep 2017
arxiv:1603.0685v3 [math.ap] 1 Sep 017 UNIQUE CONTINUATION FOR THE SCHRÖDINGER EQUATION WITH GRADIENT TERM YOUNGWOO KOH AND IHYEOK SEO Abstract. We obtain a unique continuation result for the differential
More informationSharp Sobolev Strichartz estimates for the free Schrödinger propagator
Sharp Sobolev Strichartz estimates for the free Schrödinger propagator Neal Bez, Chris Jeavons and Nikolaos Pattakos Abstract. We consider gaussian extremisability of sharp linear Sobolev Strichartz estimates
More informationExample Linear Algebra Competency Test
Example Linear Algebra Competency Test The 4 questions below are a combination of True or False, multiple choice, fill in the blank, and computations involving matrices and vectors. In the latter case,
More informationExtension and Representation of Divergence-free Vector Fields on Bounded Domains. Tosio Kato, Marius Mitrea, Gustavo Ponce, and Michael Taylor
Extension and Representation of Divergence-free Vector Fields on Bounded Domains Tosio Kato, Marius Mitrea, Gustavo Ponce, and Michael Taylor 1. Introduction Let Ω R n be a bounded, connected domain, with
More informationAn inverse scattering problem in random media
An inverse scattering problem in random media Pedro Caro Joint work with: Tapio Helin & Matti Lassas Computational and Analytic Problems in Spectral Theory June 8, 2016 Outline Introduction and motivation
More informationResolvent differences and their spectral estimates
Resolvent differences and their spectral estimates Gerd Grubb Copenhagen University M.Z. Solomyak 80 years, Imperial College September 2011 Plan 1. Spectral upper estimates 1. Spectral upper estimates
More informationProperties of the Scattering Transform on the Real Line
Journal of Mathematical Analysis and Applications 58, 3 43 (001 doi:10.1006/jmaa.000.7375, available online at http://www.idealibrary.com on Properties of the Scattering Transform on the Real Line Michael
More informationBASIC ALGORITHMS IN LINEAR ALGEBRA. Matrices and Applications of Gaussian Elimination. A 2 x. A T m x. A 1 x A T 1. A m x
BASIC ALGORITHMS IN LINEAR ALGEBRA STEVEN DALE CUTKOSKY Matrices and Applications of Gaussian Elimination Systems of Equations Suppose that A is an n n matrix with coefficents in a field F, and x = (x,,
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 informationThe inverse conductivity problem with power densities in dimension n 2
The inverse conductivity problem with power densities in dimension n 2 François Monard Guillaume Bal Dept. of Applied Physics and Applied Mathematics, Columbia University. June 19th, 2012 UC Irvine Conference
More information