Inverse source estimation problems in magnetostatics
|
|
- Bertha Clarke
- 5 years ago
- Views:
Transcription
1 Inverse source estimation problems in magnetostatics ERNSI Workshop, Lyon, 2017 Juliette Leblond INRIA Sophia Antipolis (France), Team APICS Joint work with: Laurent Baratchart, Sylvain Chevillard, Doug Hardin, Eduardo Lima, Jean-Paul Marmorat, Dmitry Ponomarev
2 Context APICS team, Inria Sophia Antipolis, France Earth, Atmospheric & Planetary Sciences Dep., MIT Cambridge, Massachusetts, USA Center for Constructive Approximation Vanderbilt University, Nashville, Tennessee, USA Planetary sciences, paleomagnetism, remanent magnetization of ancient rocks history and future of Earth magnetic field Magnetization not measurable measures of generated magnetic field inverse problems, non destructive inspection
3 How do rocks acquire magnetization? Igneous rocks (from Earth volcanoes, lava, magma; basalts) Thermoremanent magnetization, ferromagnetic particles (small magnets) follow the magnetic field:
4 How do rocks acquire magnetization? Igneous rocks (from Earth volcanoes, lava, magma; basalts) Thermoremanent magnetization, ferromagnetic particles (small magnets) follow the magnetic field:
5 How do rocks acquire magnetization? Igneous rocks (from Earth volcanoes, lava, magma; basalts) Thermoremanent magnetization, ferromagnetic particles (small magnets) follow the magnetic field:
6 How do rocks acquire magnetization? Igneous rocks (from Earth volcanoes, lava, magma; basalts) Thermoremanent magnetization, ferromagnetic particles (small magnets) follow the magnetic field:
7 How do rocks acquire magnetization? Igneous rocks (from Earth volcanoes, lava, magma; basalts) Thermoremanent magnetization, ferromagnetic particles (small magnets) follow the magnetic field:
8 How do rocks acquire magnetization? Igneous rocks (from Earth volcanoes, lava, magma; basalts) Thermoremanent magnetization, ferromagnetic particles (small magnets) follow the magnetic field: Can be subsequently altered, under high pressure or temperature
9 SQUID microscope Scanning magnetic microscopes: for weakly magnetized small samples (MIT, EAPS) pedestal + sensor sapphire window map of the vertical component of the (tiny) magnetic field on a rectangular region Q above the sample S
10 General scheme Thin sample planar (rectangle) S R 2 M = m = (m1, m 2, m 3 ) t support of unknown magnetization (source term) m B 3 on Q m on S?
11 General scheme Thin sample planar (rectangle) S R 2 M = m = (m1, m 2, m 3 ) t support of unknown magnetization (source term) m B 3 on Q m on S?
12 General scheme Thin sample planar (rectangle) S R 2 M = m = (m1, m 2, m 3 ) t support of unknown magnetization (source term) m B 3 on Q m on S?
13 Hawaiian Basalt example picture of thin sample magnetization recovery measures (nt ) of B3 on Q {h} estimated magnetization (Am2 ) m on S {0} knowing direction of m & after re-magnetization (MIT, EAPS) not always feasible...
14 Problem setting, B 3 = B 3 [m] m generates a magnetic field at x S: B[m](x) = µ 0 U[m](x) where U[m](x) = 1 m(y), x y 4π x y 3 dy S U scalar magnetic potential, µ 0 magnetic constant, 3D gradient
15 Problem setting, B 3 = B 3 [m] m generates a magnetic field at x S: B[m](x) = µ 0 U[m](x) where U[m](x) = 1 m(y), x y 4π x y 3 dy S U scalar magnetic potential, µ 0 magnetic constant, 3D gradient measurements b 3 [m] of vertical component B 3 [m] = µ 0 x3 U[m] performed on square Q R 2 at height h (incomplete data) xi partial derivative w.r.t. x i
16 Problem setting, B 3 = B 3 [m] m generates a magnetic field at x S: B[m](x) = µ 0 U[m](x) where U[m](x) = 1 m(y), x y 4π x y 3 dy S U scalar magnetic potential, µ 0 magnetic constant, 3D gradient measurements b 3 [m] of vertical component B 3 [m] = µ 0 x3 U[m] performed on square Q R 2 at height h (incomplete data) xi partial derivative w.r.t. x i inverse problems: from pointwise values of b 3 [m] on Q recover magnetization m or net moment m on S m 1 m = m 2, m i = m 3 S m i (y) dy average, m i L 2 (S) L 1 (S), i = 1, 2, 3
17 Operators Put X = (x 1, x 2, x 3 ) R 3 grad = = x1 x2 x3 with xi = xi div =, curl = Laplacian = = div (grad) = = 2 x x x3 2
18 Inverse recovery problems Poisson-Laplace equation with source term in divergence form: Maxwell, magnetostatics, time-harmonic, macroscopic U[m] = m in R 3, supp m S {0} B 3 = B 3 [m] = x3 U[m] U[m] and B 3 [m] harmonic functions in {x 3 > 0} b 3 [m] = B 3 [m] Q {h} on Q {h} m or m on S {0}? m = m(y) d y R 3, m L 2 (S, R 3 ), b 3 [m] L 2 (Q) S Integral expression for (x 1, x 2 ) Q (also in terms of Poisson and Riesz transforms) b 3 [m] (x 1, x 2 ) = µ 0 4π m 1 (y) (x 1 y 1 ) + m 2 (y) (x 2 y 2 ) + m 3 (y) x 3 x3 ( S (x1 y 1 ) 2 + (x 2 y 2 ) 2 + x3 2 ) dy 3/2 x 3 =h y = (y 1, y 2 )
19 Inverse recovery problems Data: b 3 [m] on Q net moment m i estimation: uniqueness, unstability regularization (BEP) preliminary step for magnetization recovery (non-uniqueness): mean values m i of m i on S m direction of m
20 Inverse recovery problems Linear map: b 3 : { L 2 (S, R 3 ) L 2 (Q) m B 3 [m] Q {h} = b 3 [m] Magnetization recovery (full inversion): recover m from b 3 [m] ill-posed (non uniqueness) because Ker b 3 {0}, silent sources
21 Inverse recovery problems Linear map: b 3 : { L 2 (S, R 3 ) L 2 (Q) m B 3 [m] Q {h} = b 3 [m] Magnetization recovery (full inversion): recover m from b 3 [m] ill-posed (non uniqueness) because Ker b 3 {0}, silent sources Moment recovery: well-defined i {1, 2, 3} { Ran b3 R b 3 [m] m i
22 Inverse recovery problems Linear map: b 3 : { L 2 (S, R 3 ) L 2 (Q) m B 3 [m] Q {h} = b 3 [m] Magnetization recovery (full inversion): recover m from b 3 [m] ill-posed (non uniqueness) because Ker b 3 {0}, silent sources Moment recovery: well-defined i {1, 2, 3} { Ran b3 R b 3 [m] m i Strategy: linear estimator m of given norm, once and for all find test functions φ i L 2 (Q) such that φ i, b 3 [m] L 2 (Q) m i best constrained approximation problems (BEP) Direct inversion of associated discrete problem: heavy, unstable and asymptotic formulas, for large measurement area Q
23 Strategy for net moment recovery Determine φ i L 2 (Q) such that b 3 [m], φ i Q m i = m, e i S.,. L 2 (Q),.,. L 2 (Q,R 2 ).,. Q,. Q,.,. S.,. S,. S e 1 = (1, 0, 0), e 2 = (0, 1, 0), e 3 = (0, 0, 1) on S b 3 [m], φ i Q = m, b 3 [φ i ] S with adjoint operator b 3 : L2 (Q) L 2 (S, R 3 ) to b 3 linear estimator for moment m i given bound on m S trade-off between prescribed accuracy and norm of φ i in Sobolev space, L 2 norm of gradient of φ stability w.r.t. measurement errors, robustness
24 Tools for analysis Properties of Hilbert Sobolev spaces W 1,2 0 (Ω) L 2 (Ω) Ω R 2 Integral expression of b 3 [m] in terms of Poisson and Riesz transforms Properties of operator b 3 and its adjoint b 3 e i Ker b 3 : silent sources possess vanishing net moments uniqueness of m from b 3 [m] Ran b 3 dense in L 2 (Q)
25 Net moment recovery, strategy no hope to find φ L 2 (Q) such that b 3 [m], φ Q m i = 0 m L 2 (S, R 3 ) because e i Ran b3 and b 3 [m], φ Q m i = m, b3 [φ] e i S however, density result for m L 2 (S, R 3 ) φ n W 1,2 0 (Q) s.t. b 3 [m], φ n Q m i 0 as n inf 1,2 b3 [m], φ Q m i = 0 φ W (Q) 0 but it satisfies 2 φ n Q unstability for m L 2 (S, R 3 ), we can find φ W 1,2 0 (Q) L 2 (Q) s.t. b 3 [m], φ Q m i arbitrarily small at the expense of arbitrarily large 2 φ Q needs for regularization, trade-off (error / constraint)
26 Bounded extremal problem b 3 [m], φ Q m i = m, b 3 [φ] e i S m S b 3 [φ] e i S Best constrained approximation, optimization issue M > 0 unique solution φ W 1,2 0 (Q) to (BEP) b 3 [φ ] e i S = min φ W 1,2 0 (Q), 2 φ Q M b 3 [φ] e i S it satisfies the critical point equation (CPE) and 2 φ Q = M Tykhonov type regularization b 3 b 3 [φ ] λ 2 φ = b 3 [e i ] on Q for unique λ > 0 s.t. 2 φ Q = M also, φ C 1/2 ( Q) Lagrange parameter λ 0 iff constraint M & criterion error 0 resolution algorithm: iterate in λ
27 Ongoing, solving (CPE), preliminary numerics m = (m 1, m 2, m 3 ) on S b 3 [m] on Q synthetic data (matlab) simulated m, b 3 [m] R = 2.55, s = 1.97, h = 0.27 mm, Q / S 1.3, S /h overlapping squares, finite elements P 0 on S Q 1 on Q, for W 1,2 0 (Q) basis functions ψ (CPE) Dirichlet problems for elliptic 2D PDE: φ W 1,2 0 (Q), φ Q = 0 b 3 [φ ], b 3 [ψ] S + λ 2 φ, 2 ψ Q = b 3 [e i ], ψ Q
28 Solving (CPE), preliminary numerics towards a numerical magnetometer, φ = φ ei, i = 1, 2, kb $ 3 [?e1(6)]! e1k L 2 (S;R 3 ) ke1k L 2 (S;R 3 ) kb $ 3 [?e3(6)]! e3k L 2 (S;R 3 ) ke3k L 2 (S;R 3 ) kr?e1 (6)k L 2 (Q;R 2 ) kr?e3 (6)k L 2 (Q;R 2 ) log 10(6) log 10(6) criterion error (cost) on S constraint M on Q b 3 [φ e i ] e i S 2 φ ei Q for i = 1, 3, w.r.t. log 10 (λ)
29 Solving (CPE), preliminary numerics kb $ 3 [?e1(6)]! e1k L 2 (S;R 3 ) 0.45 kb $ 3 [?e3(6)]! e1k L 2 (S;R 3 ) = 10!4 ke1k L 2 (S;R 3 ) = 10!4 ke3k L 2 (S;R 3 ) = 10! = 10! = 10! = 10! = 10! = 10! = 10!8 6 = 10! = 10!10 6 = 10!11 6 = 10! M = kr?e1 (6)k L 2 (Q;R 2 ) = 10! = 10!9 6 = 10! = 10!11 6 = 10! M = kr?e3 (6)k L 2 (Q;R 2 ) L-curves, i = 1 (left), i = 3 (right): criterion error on S w.r.t. constraint M on Q, at values of λ b3 [φe i ] e Q i w.r.t. 2 φ ei S trade-off between error and constraint
30 Solving (CPE), preliminary numerics φ ei (λ 1 ) on Q, for i = 1, 2, 3 λ 1 = φ ei (λ 2 ) on Q, for i = 1, 2, 3 λ 2 = 10 8
31 Solving (CPE), preliminary numerics errors ɛ k on component k = 1, 2, 3 of b 3 [ ] φ ei (λ 1 ) e i on S for i = 1, 2, 3 (rows) λ 1 = 10 11
32 Solving (CPE), preliminary numerics 1st row: actual values of components m i, i = 1, 2, e e e-05 λ = e e e % e e e % e e e % e e e % e e e % e e e % e e e % next rows: estimated values values b 3 [m], φ ei [λ] Q last column: relative errors on moment m
33 Ongoing, next Resolution schemes, numerical analysis (ctn), test noisy & actual data & other resolution schemes: iterative, φ n φ = φ ei n λ ϱ 2 φ n + φ n = ϱ b 3 b3 [ ] φn 1 + φn 1 + ϱ b 3 [e i ] and φ n Q = 0, n 1 in Fourier basis asymptotic estimations product of sin of separated variables Q > S or R > s Local moments determined by b 3 [m] if m of minimal L 2 (S) norm (BEP) for general e Ran b 3 \ Ran b 3 higher order moments Use m for magnetization m estimation
34 References L. Baratchart, D. Hardin, E.A. Lima, E.B. Saff, B.P. Weiss, Characterizing kernels of operators related to thin-plate magnetizations via generalizations of Hodge decompositions, Inverse Problems, 2013 L. Baratchart, S. Chevillard, J. Leblond, Silent and equivalent magnetic distributions on thin plates, Theta Series in Advanced Mathematics, to appear, hal Also, in planetary sciences, paleomagnetism: magnetizations m in L 1 (S, R 3 ) or distributions for 3D samples S, pointwise dipolar sources chondrules, also Moon rocks,... In brain imaging (medical engineering), electroencephalography (EEG), source and conductivity inverse problems, spherical geometry, software FindSources3D ( + MEG)
Estimation de paramètres électrophysiologiques pour
Introduction Estimation de paramètres électrophysiologiques pour l imagerie cérébrale Maureen Clerc Inria Sophia Antipolis, France Paris, 20 jan 2017 Séminaire Jacques-Louis Lions Maureen Clerc (Inria,
More informationMathematical Notes for E&M Gradient, Divergence, and Curl
Mathematical Notes for E&M Gradient, Divergence, and Curl In these notes I explain the differential operators gradient, divergence, and curl (also known as rotor), the relations between them, the integral
More informationUniqueness result for an inverse conductivity recovery problem with application to EEG
Rend. Istit. Mat. Univ. Trieste Volume 48 (2016), 385 406 DOI: 10.13137/2464-8728/13165 Uniqueness result for an inverse conductivity recovery problem with application to EEG Maureen Clerc, Juliette Leblond,
More informationINVERSE POTENTIAL PROBLEMS FOR DIVERGENCE OF MEASURES WITH TOTAL VARIATION REGULARIZATION
INVERSE POTENTIAL PROBLEMS FOR DIVERGENCE OF MEASURES WITH TOTAL VARIATION REGULARIZATION L. BARATCHART, C. VILLALOBOS GUILLEN, D.P. HARDIN, M.C. NORTHINGTON, AND E.B. SAFF arxiv:1809.08334v1 [math.oc]
More informationNitsche-type Mortaring for Maxwell s Equations
Nitsche-type Mortaring for Maxwell s Equations Joachim Schöberl Karl Hollaus, Daniel Feldengut Computational Mathematics in Engineering at the Institute for Analysis and Scientific Computing Vienna University
More informationUltra-high Sensitivity Moment Magnetometry of Geological Samples Using Magnetic Microscopy
Ultra-high Sensitivity Moment Magnetometry of Geological Samples Using Magnetic Microscopy Eduardo A. Lima 1 *, Benjamin P. Weiss 1 1 Department of Earth, Atmospheric, and Planetary Sciences, Massachusetts
More informationPhotometric Stereo: Three recent contributions. Dipartimento di Matematica, La Sapienza
Photometric Stereo: Three recent contributions Dipartimento di Matematica, La Sapienza Jean-Denis DUROU IRIT, Toulouse Jean-Denis DUROU (IRIT, Toulouse) 17 December 2013 1 / 32 Outline 1 Shape-from-X techniques
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 informationCreated by T. Madas VECTOR OPERATORS. Created by T. Madas
VECTOR OPERATORS GRADIENT gradϕ ϕ Question 1 A surface S is given by the Cartesian equation x 2 2 + y = 25. a) Draw a sketch of S, and describe it geometrically. b) Determine an equation of the tangent
More informationA constrained approximation problem arising in parameter identification
Linear Algebra and its Applications 351 35 (00) 487 500 www.elsevier.com/locate/laa A constrained approximation problem arising in parameter identification Birgit Jacob a, Juliette Leblond b, Jean-Paul
More informationNew estimates for the div-curl-grad operators and elliptic problems with L1-data in the half-space
New estimates for the div-curl-grad operators and elliptic problems with L1-data in the half-space Chérif Amrouche, Huy Hoang Nguyen To cite this version: Chérif Amrouche, Huy Hoang Nguyen. New estimates
More informationDetailed Program of the Workshop
Detailed Program of the Workshop Inverse Problems: Theoretical and Numerical Aspects Laboratoire de Mathématiques de Reims, December 17-19 2018 December 17 14h30 Opening of the Workhop 15h00-16h00 Mourad
More informationChapter 1 Mathematical Foundations
Computational Electromagnetics; Chapter 1 1 Chapter 1 Mathematical Foundations 1.1 Maxwell s Equations Electromagnetic phenomena can be described by the electric field E, the electric induction D, the
More informationElectrodynamics and Microwaves 3. Gradient, Curl and Divergence
1 Module 3 Gradient, Divergence and Curl 1. Introduction 2. The operators & 2 3. Gradient 4. Divergence 5. Curl 6. Mathematical expressions for gradient, divergence and curl in different coordinate systems.
More informationLecture 4: Harmonic forms
Lecture 4: Harmonic forms Jonathan Evans 29th September 2010 Jonathan Evans () Lecture 4: Harmonic forms 29th September 2010 1 / 15 Jonathan Evans () Lecture 4: Harmonic forms 29th September 2010 2 / 15
More informationSchur parametrizations of stable all-pass discrete-time systems and balanced realizations. Application to rational L 2 approximation.
Schur parametrizations of stable all-pass discrete-time systems and balanced realizations. Application to rational L 2 approximation. Martine Olivi INRIA, Sophia Antipolis 12 Mai 2004 1 Outline. Introduction
More information송석호 ( 물리학과 )
http://optics.hanyang.ac.kr/~shsong 송석호 ( 물리학과 ) Introduction to Electrodynamics, David J. Griffiths Review: 1. Vector analysis 2. Electrostatics 3. Special techniques 4. Electric fields in mater 5. Magnetostatics
More informationInverse Problems and Optimal Design in Electricity and Magnetism
Inverse Problems and Optimal Design in Electricity and Magnetism P. Neittaanmäki Department of Mathematics, University of Jyväskylä M. Rudnicki Institute of Electrical Engineering, Warsaw and A. Savini
More informationREMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID
REMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID DRAGOŞ IFTIMIE AND JAMES P. KELLIHER Abstract. In [Math. Ann. 336 (2006), 449-489] the authors consider the two dimensional
More informationInequalities of Babuška-Aziz and Friedrichs-Velte for differential forms
Inequalities of Babuška-Aziz and Friedrichs-Velte for differential forms Martin Costabel Abstract. For sufficiently smooth bounded plane domains, the equivalence between the inequalities of Babuška Aziz
More informationINTRODUCTION TO FINITE ELEMENT METHODS ON ELLIPTIC EQUATIONS LONG CHEN
INTROUCTION TO FINITE ELEMENT METHOS ON ELLIPTIC EQUATIONS LONG CHEN CONTENTS 1. Poisson Equation 1 2. Outline of Topics 3 2.1. Finite ifference Method 3 2.2. Finite Element Method 3 2.3. Finite Volume
More information? D. 3 x 2 2 y. D Pi r ^ 2 h, r. 4 y. D 3 x ^ 3 2 y ^ 2, y, y. D 4 x 3 y 2 z ^5, z, 2, y, x. This means take partial z first then partial x
PartialsandVectorCalclulus.nb? D D f, x gives the partial derivative f x. D f, x, n gives the multiple derivative n f x n. D f, x, y, differentiates f successively with respect to x, y,. D f, x, x 2, for
More informationA mathematical model and inversion procedure for Magneto-Acousto-Electric Tomography (MAET)
A mathematical model and inversion procedure for Magneto-Acousto-Electric Tomography (MAET) Leonid Kunyansky University of Arizona, Tucson, AZ Suppported in part by NSF grant DMS-0908243 NSF grant "we
More informationMixed Finite Elements Method
Mixed Finite Elements Method A. Ratnani 34, E. Sonnendrücker 34 3 Max-Planck Institut für Plasmaphysik, Garching, Germany 4 Technische Universität München, Garching, Germany Contents Introduction 2. Notations.....................................
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 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 informationDivergence Theorem and Its Application in Characterizing
Divergence Theorem and Its Application in Characterizing Fluid Flow Let v be the velocity of flow of a fluid element and ρ(x, y, z, t) be the mass density of fluid at a point (x, y, z) at time t. Thus,
More informationChapter 2. Vector Calculus. 2.1 Directional Derivatives and Gradients. [Bourne, pp ] & [Anton, pp ]
Chapter 2 Vector Calculus 2.1 Directional Derivatives and Gradients [Bourne, pp. 97 104] & [Anton, pp. 974 991] Definition 2.1. Let f : Ω R be a continuously differentiable scalar field on a region Ω R
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationDouble Layer Potentials on Polygons and Pseudodifferential Operators on Lie Groupoids
Double Layer Potentials on Polygons and Pseudodifferential Operators on Lie Groupoids joint work with Hengguang Li Yu Qiao School of Mathematics and Information Science Shaanxi Normal University Xi an,
More informationLeast-Squares Finite Element Methods
Pavel В. Bochev Max D. Gunzburger Least-Squares Finite Element Methods Spri ringer Contents Part I Survey of Variational Principles and Associated Finite Element Methods 1 Classical Variational Methods
More informationClassical Fourier Analysis
Loukas Grafakos Classical Fourier Analysis Second Edition 4y Springer 1 IP Spaces and Interpolation 1 1.1 V and Weak IP 1 1.1.1 The Distribution Function 2 1.1.2 Convergence in Measure 5 1.1.3 A First
More informationGradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice
1 Lecture Notes, HCI, 4.1.211 Chapter 2 Gradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice Bastian Goldlücke Computer Vision Group Technical University of Munich 2 Bastian
More informationScale-Invariance of Support Vector Machines based on the Triangular Kernel. Abstract
Scale-Invariance of Support Vector Machines based on the Triangular Kernel François Fleuret Hichem Sahbi IMEDIA Research Group INRIA Domaine de Voluceau 78150 Le Chesnay, France Abstract This paper focuses
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationFast Iterative Solution of Saddle Point Problems
Michele Benzi Department of Mathematics and Computer Science Emory University Atlanta, GA Acknowledgments NSF (Computational Mathematics) Maxim Olshanskii (Mech-Math, Moscow State U.) Zhen Wang (PhD student,
More informationSolving the Laplace Equation by Meshless Collocation Using Harmonic Kernels
Solving the Laplace Equation by Meshless Collocation Using Harmonic Kernels Robert Schaback February 13, 2008 Abstract We present a meshless technique which can be seen as an alternative to the Method
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 informationDivergence Theorem December 2013
Divergence Theorem 17.3 11 December 2013 Fundamental Theorem, Four Ways. b F (x) dx = F (b) F (a) a [a, b] F (x) on boundary of If C path from P to Q, ( φ) ds = φ(q) φ(p) C φ on boundary of C Green s Theorem:
More informationMulti-Parameter Riesz Commutators
Multi-Parameter Riesz Commutators Brett D. Wick Vanderbilt University Department of Mathematics Sixth Prairie Analysis Seminar University of Kansas October 14th, 2006 B. D. Wick (Vanderbilt University)
More informationDefinition and basic properties of heat kernels I, An introduction
Definition and basic properties of heat kernels I, An introduction Zhiqin Lu, Department of Mathematics, UC Irvine, Irvine CA 92697 April 23, 2010 In this lecture, we will answer the following questions:
More informationSpotlight on Laplace s Equation
16 Spotlight on Laplace s Equation Reference: Sections 1.1,1.2, and 1.5. Laplace s equation is the undriven, linear, second-order PDE 2 u = (1) We defined diffusivity on page 587. where 2 is the Laplacian
More information1 Fundamentals. 1.1 Overview. 1.2 Units: Physics 704 Spring 2018
Physics 704 Spring 2018 1 Fundamentals 1.1 Overview The objective of this course is: to determine and fields in various physical systems and the forces and/or torques resulting from them. The domain of
More informationMATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Fall 2011 Professor: Jared Speck
MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Fall 20 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations
More informationMATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Spring 2018 Professor: Jared Speck
MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Spring 208 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations
More informationDivergence Theorem Fundamental Theorem, Four Ways. 3D Fundamental Theorem. Divergence Theorem
Divergence Theorem 17.3 11 December 213 Fundamental Theorem, Four Ways. b F (x) dx = F (b) F (a) a [a, b] F (x) on boundary of If C path from P to Q, ( φ) ds = φ(q) φ(p) C φ on boundary of C Green s Theorem:
More informationClassical Field Theory: Electrostatics-Magnetostatics
Classical Field Theory: Electrostatics-Magnetostatics April 27, 2010 1 1 J.D.Jackson, Classical Electrodynamics, 2nd Edition, Section 1-5 Electrostatics The behavior of an electrostatic field can be described
More informationFinite Element Exterior Calculus and Applications Part V
Finite Element Exterior Calculus and Applications Part V Douglas N. Arnold, University of Minnesota Peking University/BICMR August 15 18, 2015 De Rham complex 0 H 1 (Ω) grad H(curl, Ω) curl H(div, Ω) div
More informationAND NONLINEAR SCIENCE SERIES. Partial Differential. Equations with MATLAB. Matthew P. Coleman. CRC Press J Taylor & Francis Croup
CHAPMAN & HALL/CRC APPLIED MATHEMATICS AND NONLINEAR SCIENCE SERIES An Introduction to Partial Differential Equations with MATLAB Second Edition Matthew P Coleman Fairfield University Connecticut, USA»C)
More informationELECTROMAGNETIC FIELDS AND RELATIVISTIC PARTICLES
ELECTROMAGNETIC FIELDS AND RELATIVISTIC PARTICLES Emil J. Konopinski Professor of Physics Indiana University McGraw-Hill Book Company New York St. Louis San Francisco Auckland Bogota Hamburg Johannesburg
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Electrostatic II Notes: Most of the material presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartolo, Chap... Mathematical Considerations.. The Fourier series and the Fourier
More informationGeneral solution of the Inhomogeneous Div-Curl system and Consequences
General solution o the Inhomogeneous Div-Curl system and Consequences Briceyda Berenice Delgado López Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional March 07 2017 CINVESTAV
More informationModule III: Partial differential equations and optimization
Module III: Partial differential equations and optimization Martin Berggren Department of Information Technology Uppsala University Optimization for differential equations Content Martin Berggren (UU)
More informationON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS
Bendikov, A. and Saloff-Coste, L. Osaka J. Math. 4 (5), 677 7 ON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS ALEXANDER BENDIKOV and LAURENT SALOFF-COSTE (Received March 4, 4)
More informationSCHRIFTENREIHE DER FAKULTÄT FÜR MATHEMATIK. On Maxwell s and Poincaré s Constants. Dirk Pauly SM-UDE
SCHIFTENEIHE DE FAKULTÄT FÜ MATHEMATIK On Maxwell s and Poincaré s Constants by Dirk Pauly SM-UDE-772 2013 On Maxwell s and Poincaré s Constants Dirk Pauly November 11, 2013 Dedicated to Sergey Igorevich
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 informationMaxwell s equations in Carnot groups
Maxwell s equations in Carnot groups B. Franchi (U. Bologna) INDAM Meeting on Geometric Control and sub-riemannian Geometry Cortona, May 21-25, 2012 in honor of Andrey Agrachev s 60th birthday Researches
More informationShort title: Total FETI. Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ Ostrava, Czech Republic
Short title: Total FETI Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ-70833 Ostrava, Czech Republic mail: zdenek.dostal@vsb.cz fax +420 596 919 597 phone
More informationRobust numerical algorithms based on analytic approximation for the solution of inverse problems in annular domains
Robust numerical algorithms based on analytic approximation for the solution of inverse problems in annular domains Mohamed Jaoua, Juliette Leblond, Moncef Mahjoub, Jonathan Partington To cite this version:
More informationCrouzeix-Velte Decompositions and the Stokes Problem
Crouzeix-Velte Decompositions and te Stokes Problem PD Tesis Strauber Györgyi Eötvös Loránd University of Sciences, Insitute of Matematics, Matematical Doctoral Scool Director of te Doctoral Scool: Dr.
More informationSurfaces in spheres: biharmonic and CMC
Surfaces in spheres: biharmonic and CMC E. Loubeau (joint with C. Oniciuc) Université de Bretagne Occidentale, France May 2014 Biharmonic maps Definition Bienergy functional at φ : (M, g) (N, h) E 2 (φ)
More informationsparse and low-rank tensor recovery Cubic-Sketching
Sparse and Low-Ran Tensor Recovery via Cubic-Setching Guang Cheng Department of Statistics Purdue University www.science.purdue.edu/bigdata CCAM@Purdue Math Oct. 27, 2017 Joint wor with Botao Hao and Anru
More informationTravelling bubbles in a moving boundary problem of Hele-Shaw type
Travelling bubbles in a moving boundary problem of Hele-Shaw type G. Prokert, TU Eindhoven Center for Analysis, Scientific Computing, and Applications (CASA) g.prokert@tue.nl Joint work with M. Günther,
More informationFull-wave simulation of tokamak plasma heating
Full-wave simulation of tokamak plasma heating Institut Elie Cartan de Nancy, UMR 7502, Université de Lorraine, CNRS & INRIA (project-team CALVI). B.P. 70239, 54506 Vandœuvre les Nancy Cedex, France. Université
More informationInstitut de Recherche MAthématique de Rennes
LMS Durham Symposium: Computational methods for wave propagation in direct scattering. - July, Durham, UK The hp version of the Weighted Regularization Method for Maxwell Equations Martin COSTABEL & Monique
More informationIndex. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2
Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604
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 informationKirchhoff, Fresnel, Fraunhofer, Born approximation and more
Kirchhoff, Fresnel, Fraunhofer, Born approximation and more Oberseminar, May 2008 Maxwell equations Or: X-ray wave fields X-rays are electromagnetic waves with wave length from 10 nm to 1 pm, i.e., 10
More informationRegularization on Discrete Spaces
Regularization on Discrete Spaces Dengyong Zhou and Bernhard Schölkopf Max Planck Institute for Biological Cybernetics Spemannstr. 38, 72076 Tuebingen, Germany {dengyong.zhou, bernhard.schoelkopf}@tuebingen.mpg.de
More informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 2198-5855 On the divergence constraint in mixed finite element methods for incompressible
More informationA posteriori error estimates in FEEC for the de Rham complex
A posteriori error estimates in FEEC for the de Rham complex Alan Demlow Texas A&M University joint work with Anil Hirani University of Illinois Urbana-Champaign Partially supported by NSF DMS-1016094
More informationIntegrated interpretation of multimodal geophysical data for exploration of geothermal resources Case study: Yamagawa geothermal field in Japan
Integrated interpretation of multimodal geophysical data for exploration of geothermal resources Case study: Yamagawa geothermal field in Japan Masashi Endo*(TechnoImaging), Alex Gribenko (TechnoImaging
More informationMATH 333: Partial Differential Equations
MATH 333: Partial Differential Equations Problem Set 9, Final version Due Date: Tues., Nov. 29, 2011 Relevant sources: Farlow s book: Lessons 9, 37 39 MacCluer s book: Chapter 3 44 Show that the Poisson
More information1 Curvilinear Coordinates
MATHEMATICA PHYSICS PHYS-2106/3 Course Summary Gabor Kunstatter, University of Winnipeg April 2014 1 Curvilinear Coordinates 1. General curvilinear coordinates 3-D: given or conversely u i = u i (x, y,
More informationPDEs in Image Processing, Tutorials
PDEs in Image Processing, Tutorials Markus Grasmair Vienna, Winter Term 2010 2011 Direct Methods Let X be a topological space and R: X R {+ } some functional. following definitions: The mapping R is lower
More informationNecessary conditions for convergence rates of regularizations of optimal control problems
Necessary conditions for convergence rates of regularizations of optimal control problems Daniel Wachsmuth and Gerd Wachsmuth Johann Radon Institute for Computational and Applied Mathematics RICAM), Austrian
More informationThe Data-Dependence graph of Adjoint Codes
The Data-Dependence graph of Adjoint Codes Laurent Hascoët INRIA Sophia-Antipolis November 19th, 2012 (INRIA Sophia-Antipolis) Data-Deps of Adjoints 19/11/2012 1 / 14 Data-Dependence graph of this talk
More informationMA 441 Advanced Engineering Mathematics I Assignments - Spring 2014
MA 441 Advanced Engineering Mathematics I Assignments - Spring 2014 Dr. E. Jacobs The main texts for this course are Calculus by James Stewart and Fundamentals of Differential Equations by Nagle, Saff
More informationProblems of Corner Singularities
Stuttgart 98 Problems of Corner Singularities Monique DAUGE Institut de Recherche MAthématique de Rennes Problems of Corner Singularities 1 Vertex and edge singularities For a polyhedral domain Ω and a
More informationMETHODS FOR SOLVING MATHEMATICAL PHYSICS PROBLEMS
METHODS FOR SOLVING MATHEMATICAL PHYSICS PROBLEMS V.I. Agoshkov, P.B. Dubovski, V.P. Shutyaev CAMBRIDGE INTERNATIONAL SCIENCE PUBLISHING Contents PREFACE 1. MAIN PROBLEMS OF MATHEMATICAL PHYSICS 1 Main
More informationPhysics 6303 Lecture 3 August 27, 2018
Physics 6303 Lecture 3 August 27, 208 LAST TIME: Vector operators, divergence, curl, examples of line integrals and surface integrals, divergence theorem, Stokes theorem, index notation, Kronecker delta,
More informationELECTROMAGNETISM. Second Edition. I. S. Grant W. R. Phillips. John Wiley & Sons. Department of Physics University of Manchester
ELECTROMAGNETISM Second Edition I. S. Grant W. R. Phillips Department of Physics University of Manchester John Wiley & Sons CHICHESTER NEW YORK BRISBANE TORONTO SINGAPORE Flow diagram inside front cover
More informationSparse Recovery Beyond Compressed Sensing
Sparse Recovery Beyond Compressed Sensing Carlos Fernandez-Granda www.cims.nyu.edu/~cfgranda Applied Math Colloquium, MIT 4/30/2018 Acknowledgements Project funded by NSF award DMS-1616340 Separable Nonlinear
More informationMultigrid Methods for Maxwell s Equations
Multigrid Methods for Maxwell s Equations Jintao Cui Institute for Mathematics and Its Applications University of Minnesota Outline Nonconforming Finite Element Methods for a Two Dimensional Curl-Curl
More informationELECTROMAGNETIC SCATTERING FROM PERTURBED SURFACES. Katharine Ott Advisor: Irina Mitrea Department of Mathematics University of Virginia
ELECTROMAGNETIC SCATTERING FROM PERTURBED SURFACES Katharine Ott Advisor: Irina Mitrea Department of Mathematics University of Virginia Abstract This paper is concerned with the study of scattering of
More informationElliptic Regularity. Throughout we assume all vector bundles are smooth bundles with metrics over a Riemannian manifold X n.
Elliptic Regularity Throughout we assume all vector bundles are smooth bundles with metrics over a Riemannian manifold X n. 1 Review of Hodge Theory In this note I outline the proof of the following Fundamental
More informationObservation and Control for Operator Semigroups
Birkhäuser Advanced Texts Basler Lehrbücher Observation and Control for Operator Semigroups Bearbeitet von Marius Tucsnak, George Weiss Approx. 496 p. 2009. Buch. xi, 483 S. Hardcover ISBN 978 3 7643 8993
More informationModeling Magnetic Anomalies with MATLAB. Geodynamics, FA2012 Nate O Flaherty James Holmes
Modeling Magnetic Anomalies with MATLAB Geodynamics, FA2012 Nate O Flaherty James Holmes Introduction History o 1912 Wegner publicly advocates theory of Continental Drift 1943 G.G. Simpson publishes rebuttal,
More informationNumerical Solution of Laplace Equation By Gilberto E. Urroz, October 2004
Numerical Solution of Laplace Equation By Gilberto E. Urroz, October 004 Laplace equation governs a variety of equilibrium physical phenomena such as temperature distribution in solids, electrostatics,
More informationCHAPTER 2. COULOMB S LAW AND ELECTRONIC FIELD INTENSITY. 2.3 Field Due to a Continuous Volume Charge Distribution
CONTENTS CHAPTER 1. VECTOR ANALYSIS 1. Scalars and Vectors 2. Vector Algebra 3. The Cartesian Coordinate System 4. Vector Cartesian Coordinate System 5. The Vector Field 6. The Dot Product 7. The Cross
More informationHardy spaces associated to operators satisfying Davies-Gaffney estimates and bounded holomorphic functional calculus.
Hardy spaces associated to operators satisfying Davies-Gaffney estimates and bounded holomorphic functional calculus. Xuan Thinh Duong (Macquarie University, Australia) Joint work with Ji Li, Zhongshan
More informationRegularity Theory a Fourth Order PDE with Delta Right Hand Side
Regularity Theory a Fourth Order PDE with Delta Right Hand Side Graham Hobbs Applied PDEs Seminar, 29th October 2013 Contents Problem and Weak Formulation Example - The Biharmonic Problem Regularity Theory
More informationOn Acceleration with Noise-Corrupted Gradients. + m k 1 (x). By the definition of Bregman divergence:
A Omitted Proofs from Section 3 Proof of Lemma 3 Let m x) = a i On Acceleration with Noise-Corrupted Gradients fxi ), u x i D ψ u, x 0 ) denote the function under the minimum in the lower bound By Proposition
More information2D X-Ray Tomographic Reconstruction From Few Projections
2D X-Ray Tomographic Reconstruction From Few Projections Application of Compressed Sensing Theory CEA-LID, Thalès, UJF 6 octobre 2009 Introduction Plan 1 Introduction 2 Overview of Compressed Sensing Theory
More informationS12.1 SOLUTIONS TO PROBLEMS 12 (ODD NUMBERS)
OLUTION TO PROBLEM 2 (ODD NUMBER) 2. The electric field is E = φ = 2xi + 2y j and at (2, ) E = 4i + 2j. Thus E = 2 5 and its direction is 2i + j. At ( 3, 2), φ = 6i + 4 j. Thus the direction of most rapid
More informationMULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday.
MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* DOUGLAS N ARNOLD, RICHARD S FALK, and RAGNAR WINTHER Dedicated to Professor Jim Douglas, Jr on the occasion of his seventieth birthday Abstract
More informationTime-Varying Systems; Maxwell s Equations
Time-Varying Systems; Maxwell s Equations 1. Faraday s law in differential form 2. Scalar and vector potentials; the Lorenz condition 3. Ampere s law with displacement current 4. Maxwell s equations 5.
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 informationPartial Differential Equations
Partial Differential Equations Xu Chen Assistant Professor United Technologies Engineering Build, Rm. 382 Department of Mechanical Engineering University of Connecticut xchen@engr.uconn.edu Contents 1
More information21 Laplace s Equation and Harmonic Functions
2 Laplace s Equation and Harmonic Functions 2. Introductory Remarks on the Laplacian operator Given a domain Ω R d, then 2 u = div(grad u) = in Ω () is Laplace s equation defined in Ω. If d = 2, in cartesian
More information