Towards a Dimensionally Adaptive Inversion of the Magnetotelluric Problem
|
|
- Sheena Henry
- 5 years ago
- Views:
Transcription
1 1 Towards a Dimensionally Adaptive Inversion of the Magnetotelluric Problem J. Alvarez-Aramberri 1 D. Pardo 2 H. Barucq 3 1 University of the Basque Country (UPV/EHU), Bilbao, Spain, Inria team-project Magique-3D, France 2 University of the Basque Country (UPV/EHU), BCAM, and Ikerbasque, Bilbao, Spain. 3 EPC Magique-3D, Inria, LMA, University of Pau, France Jul 18, 2014 The Third BCAM Workshop on Computational Mathematics, BCAM.
2 The Magnetotelluric Method 2
3 The Magnetotelluric Method 3 The solar wind produces a radiation pressure that causes a compression on the day-side and a tail on the night-side onto the magnetosphere and due to this interaction, hydromagnetic waves are created. When those waves reach the ionosphere, they induce an EM field that works as power source in magnetotellurics. Range (distance): From few meters to hundreds of kilometers. Range (frequency): Between Hz. Applications: Hydrocarbon (oil and gas) exploration, geothermal exploration, earthquake precursor prediction, mining exploration, as well as hydrocarbon and groundwater monitoring.
4 The Magnetotelluric Method 4 Figure: Typical MT problem.
5 The Magnetotelluric Method 5 Figure: Sketch of the MT method.
6 Motivation 6 Figure: 2D MT problem. Blue rectangle: Natural source. Red crosses: Receivers.
7 Motivation 7 Figure: Primary Problem: 1D layered media with known analytical solution E 1D.
8 Motivation 8 Figure: Secondary Problem: 2D media. Blue rectangle: New source. Red crosses: Receivers.
9 Outline 9 1 Mathematical Modelization: Direct Problem 2 Inverse Problem: Dimensionally Adaptive Algorithm 3 Numerical Results
10 Outline 10 1 Mathematical Modelization: Direct Problem 2 Inverse Problem: Dimensionally Adaptive Algorithm 3 Numerical Results
11 Direct Problem 11 Maxwell s equations in frequency domain E = iωµh H = (ρ 1 + iωɛ)e + J imp (ɛe) = ρ f (µh) = 0, When E = E(x, z) and H = H(x, z) 2 uncoupled modes: Transverse Electric (TE): (E y, H x, H z ) Transverse Magnetic (TM): (H y, E x, E z ) ( ) µ 1 E k 2 E = jωj imp.
12 Direct Problem 12 We focus on the TE mode and we solve the equation for E y with a hp-finite Element Method (FEM): where k 2 = ω 2 ε jωρ 1. µ 1 E y k 2 E y = jωj imp y, (1) Multiply equation (1) by the complex conjugate of a test function F V (Ω) = H 1 Γ D (Ω) = {F L 2 (Ω) : F ΓD = 0, F L 2 (Ω)} is the space of admissible test functions and integrating by parts (omitting the subscript y from now on): Find E E D + V (Ω), such that: ( F ) T µ 1 E F k 2 E = jω Ω Ω Ω F J imp F V (Ω),
13 Direct Problem 13 With the state scalar-valued solution function E, we compute the Quantity of Interest L i (E), a linear and continuous functional in E associated to the i-th receiver and defined as: L i (E) = 1 Ω R i Ω R i E dω, where Ω R i is the domain occupied by the i-th receiver. From Maxwell s equations: H x = 1 E jωµ z. The impedance Z i = Z i yx and the apparent resistivity ρ i ap are: Z i = Li (E) L i (H x ) ρ i ap = Zi ωµ.
14 Outline 14 1 Mathematical Modelization: Direct Problem 2 Inverse Problem: Dimensionally Adaptive Algorithm 3 Numerical Results
15 Definition of the Inverse Problem 15 We define the following Misfit Function: C(σ) = 1 2M M g i (ρ) dobs i 2 i g(ρ) = {g i (ρ)} M i=1 is the vector with the i-th component being equal to the quantity of interest corresponding to the i-th receiver. In particular, g i (ρ) = Z i (ρ) or g i (ρ) = ρ i ap(σ) d = {dobs i }M i=1 is the vector of EM measurements obtained at different receivers.
16 Definition of the Inverse Problem 16 Piecewise constant distribution of conductivity: P ρ(x) = ρ i χ i (x), where χ i (x) are Heaviside shape functions. i Model Parameters: ρ = {ρ 1,..., ρ P } Problem to Solve: ρ = arg min l ρ u C(ρ), where the vectors l and u represent lower and upper bounds on the resistivity.
17 Solving the Inverse Problem 17 Solve Direct MT Problem Experimental data at receivers ρ Numerical solution at receivers Comparison Update the value of ρ ρ : Solution to the IP
18 Solving the Inverse Problem 18 Solve Direct MT Problem Experimental data at receivers ρ Numerical solution at receivers Comparison Update the value of ρ ρ : Solution to the IP
19 Dimensionally Adaptive Algorithm 19 Consider a reference 1D resistivity model ρ 1D, with known (analytical) solution. Let Ω 2D be the domain where the 2D inhomogeneities are located. Then, the conductivity distribution can be represented as the sum ρ = ρ 1D + ρ 2D, being ρ 2D zero outside Ω 2D. Divide the solution into the primary E 1D and secondary E 2D fields : E = E 1D + E 2D.
20 Dimensionally Adaptive Algorithm 20 Figure: 2D MT problem. Blue rectangle: Natural source. Red crosses: Receivers.
21 Dimensionally Adaptive Algorithm 21 Figure: Primary Problem: 1D layered media with known analytical solution E 1D.
22 Dimensionally Adaptive Algorithm 22 Figure: Secondary Problem: Solution E 2D computed with hp-fem. Blue rectangle: New source. Red crosses: Receivers.
23 Dimensionally Adaptive Algorithm 23 Compute E 1D with an analytic solution. Compute E 2D solving the equation for the secondary field with the hp-fem. ) (µ 1 E 2D k 2 E 2D = (ρ 2D ) 1 E 1D, Its variational formulation is given by (E 1D already known): Find E 2D E D + V, such that: ( F ) T µ 1 E 2D F k 2 E 2D = Ω Ω Ω F (ρ 2D ) 1 E 1D, F V. We obtain E = E 1D + E 2D.
24 Dimensionally Adaptive Algorithm 24 Define a new Misfit Function: C(ρ) = 1 2M M g i (ρ) dobs i 2 + β(ρ ρ 1D) i where ρ 1D is the solution of the 1D Inverse Problem. Step1: Solve the Inverse Problem for the 1D layered media. Step2: Solve the Inverse Problem for the 2D problem. Use the secondary field formulation for the direct problem. Use the 1D solution of the Inverse Problem as a regularization parameter.
25 Advantages of the New Approach 25 MT Direct Problem is computationally cheaper. We only need to accurately solve the secondary field variations, which in general are more localized than the total fields. Therefore, it is generally possible to use a coarser grid. Smaller computational domain: With this formulation, we avoid modeling the original source, which in MT implies to deal with big computational domains. Only the 46% of elements needed. Nonlinear Inverse Problems are usually ill-posed. It allows the design of more robust algorithms.
26 Outline 26 1 Mathematical Modelization: Direct Problem 2 Inverse Problem: Dimensionally Adaptive Algorithm 3 Numerical Results
27 Model Problem: Magnetotelluric Problem 27 Figure: TE mode driven by the impressed electric current J imp = ωŷ ρ 1 ρ 2 ρ 3 Model Model Model Model Table: Different models for the formation of the subsurface (Ohm-m).
28 Validation of the Solution 28 Apparent resistivity (Ohm m) Model 1 Model 2 Model 3 Model Frequency (Hz.) 10 0 Relative error in percent Model 1 Model 2 Model 3 Model Frequency (Hz.) 10 0 Figure: Apparent resistivity against frequency for the numerical solution (left) and the relative error for different subsurface formations for a range of frequencies (right). ρ 1 ρ 2 ρ 3 Model Model Model Model
29 Formation with an Inhomogeneity (Target) Figure: TE mode driven by J imp = ωŷ. Target s width: 10km. ρ 1 ρ 2 ρ 3 ρ 4 Layered Layered Target Table: Different models of the formation of the subsurface. (Ohm-m.) 29
30 Formation with an Inhomogeneity (Target) Figure: TE mode driven by J imp = ωŷ. Target s width: 10km. ρ 1 ρ 2 ρ 3 ρ 4 Layered Layered Target Table: Different models of the formation of the subsurface. (Ohm-m.) 30
31 Solution with the Secondary Field Formulation 31 Apparent resistivity (Ohm m) Layered 1 Layered 2 Target Position (km.) Figure: Apparent resistivity at different distances for three different subsurface formations. In all cases ω = 10 3 Hz.
32 Conclusions / Future (Current) Work 32 The secondary field formulation provides accurate solutions to the MT problem. These solutions are obtained with lower computational cost (only the half of the elements needed). Implement the adaptive dimensional algorithm for the inversion using the secondary field formulation for the Direct Problem.
Impedance Transmission Conditions for the Electric Potential across a Highly Conductive Casing
Impedance Transmission Conditions for the Electric Potential across a Highly Conductive Casing Hélène Barucq, Aralar Erdozain, David Pardo, Victor Péron To cite this version: Hélène Barucq, Aralar Erdozain,
More informationMagnetotelluric (MT) Method
Magnetotelluric (MT) Method Dr. Hendra Grandis Graduate Program in Applied Geophysics Faculty of Mining and Petroleum Engineering ITB Geophysical Methods Techniques applying physical laws (or theory) to
More informationFast Inversion of Logging-While-Drilling (LWD) Resistivity Measurements
Fast Inversion of Logging-While-Drilling (LWD) Resistivity Measurements David Pardo 1 Carlos Torres-Verdín 2 1 University of the Basque Country (UPV/EHU) and Ikerbasque, Bilbao, Spain. 2 The University
More informationA Brief Introduction to Magnetotellurics and Controlled Source Electromagnetic Methods
A Brief Introduction to Magnetotellurics and Controlled Source Electromagnetic Methods Frank Morrison U.C. Berkeley With the help of: David Alumbaugh Erika Gasperikova Mike Hoversten Andrea Zirilli A few
More informationAntennas and Propagation. Chapter 2: Basic Electromagnetic Analysis
Antennas and Propagation : Basic Electromagnetic Analysis Outline Vector Potentials, Wave Equation Far-field Radiation Duality/Reciprocity Transmission Lines Antennas and Propagation Slide 2 Antenna Theory
More informationAn Experimental Magnetotelluric Short Sounding During High Sun Spot Activity Cycle in 2013
Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 An Experimental Magnetotelluric Short Sounding During High Sun Spot Activity Cycle in 2013 Lendriadi AGUNG and Imam B.
More informationA Nested Dissection Parallel Direct Solver. for Simulations of 3D DC/AC Resistivity. Measurements. Maciej Paszyński (1,2)
A Nested Dissection Parallel Direct Solver for Simulations of 3D DC/AC Resistivity Measurements Maciej Paszyński (1,2) David Pardo (2), Carlos Torres-Verdín (2) (1) Department of Computer Science, AGH
More informationMT Prospecting. Map Resistivity. Determine Formations. Determine Structure. Targeted Drilling
MT Prospecting Map Resistivity Determine Formations Determine Structure Targeted Drilling Cross-sectional interpretation before and after an MT survey of a mineral exploration prospect containing volcanic
More information5 Electromagnetic Waves
5 Electromagnetic Waves 5.1 General Form for Electromagnetic Waves. In free space, Maxwell s equations are: E ρ ɛ 0 (5.1.1) E + B 0 (5.1.) B 0 (5.1.3) B µ 0 ɛ 0 E µ 0 J (5.1.4) In section 4.3 we derived
More informationMultiphysics, Inversion, and Petroleum
CIC biogune Research Line II: Multiphysics, Inversion, and Petroleum Research Professor at BCAM Team: D. Pardo, I. Garay, I. Andonegui, J. Álvarez Collaborators: P. de la Hoz, M. Paszynski, L.E. García-Castillo,
More informationUniform Plane Waves Page 1. Uniform Plane Waves. 1 The Helmholtz Wave Equation
Uniform Plane Waves Page 1 Uniform Plane Waves 1 The Helmholtz Wave Equation Let s rewrite Maxwell s equations in terms of E and H exclusively. Let s assume the medium is lossless (σ = 0). Let s also assume
More informationANTENNAS. Vector and Scalar Potentials. Maxwell's Equations. E = jωb. H = J + jωd. D = ρ (M3) B = 0 (M4) D = εe
ANTENNAS Vector and Scalar Potentials Maxwell's Equations E = jωb H = J + jωd D = ρ B = (M) (M) (M3) (M4) D = εe B= µh For a linear, homogeneous, isotropic medium µ and ε are contant. Since B =, there
More informationELE3310: Basic ElectroMagnetic Theory
A summary for the final examination EE Department The Chinese University of Hong Kong November 2008 Outline Mathematics 1 Mathematics Vectors and products Differential operators Integrals 2 Integral expressions
More informationECE 107: Electromagnetism
ECE 107: Electromagnetism Set 7: Dynamic fields Instructor: Prof. Vitaliy Lomakin Department of Electrical and Computer Engineering University of California, San Diego, CA 92093 1 Maxwell s equations Maxwell
More informationThe sun small resistivity enables it to absorb large amount of energy from weak magnetic fields.
The sun small resistivity enables it to absorb large amount of energy from weak magnetic fields. The solar cycle is the mechanism by witch the galactic center transfer energy to the sun by magnetic fields.
More informationPhysics 506 Winter 2004
Physics 506 Winter 004 G. Raithel January 6, 004 Disclaimer: The purpose of these notes is to provide you with a general list of topics that were covered in class. The notes are not a substitute for reading
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION Supplementary online material for Bai et al., (2). EHS3D MT data collection Broadband magnetotelluric (MT) data were recorded on profiles P, P2 and P4 in the frequency band -.5
More informationFDTD for 1D wave equation. Equation: 2 H Notations: o o. discretization. ( t) ( x) i i i i i
FDTD for 1D wave equation Equation: 2 H = t 2 c2 2 H x 2 Notations: o t = nδδ, x = iδx o n H nδδ, iδx = H i o n E nδδ, iδx = E i discretization H 2H + H H 2H + H n+ 1 n n 1 n n n i i i 2 i+ 1 i i 1 = c
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 informationScalar electromagnetic integral equations
Scalar electromagnetic integral equations Uday K Khankhoje Abstract This brief note derives the two dimensional scalar electromagnetic integral equation starting from Maxwell s equations, and shows how
More informationEM waves: energy, resonators. Scalar wave equation Maxwell equations to the EM wave equation A simple linear resonator Energy in EM waves 3D waves
EM waves: energy, resonators Scalar wave equation Maxwell equations to the EM wave equation A simple linear resonator Energy in EM waves 3D waves Simple scalar wave equation 2 nd order PDE 2 z 2 ψ (z,t)
More informationACM 106a: Lecture 1 Agenda
1 ACM 106a: Lecture 1 Agenda Introduction to numerical linear algebra Common problems First examples Inexact computation What is this course about? 2 Typical numerical linear algebra problems Systems of
More information3D MAGNETOTELLURIC CHARACTERIZATION OF THE COSO GEOTHERMAL FIELD
PROCEEDINGS, Thirtieth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 31-February 2, 2005 SGP-TR-176 3D MAGNETOTELLURIC CHARACTERIZATION OF THE COSO GEOTHERMAL
More information1 Fundamentals of laser energy absorption
1 Fundamentals of laser energy absorption 1.1 Classical electromagnetic-theory concepts 1.1.1 Electric and magnetic properties of materials Electric and magnetic fields can exert forces directly on atoms
More informationAn inverse problem for eddy current equations
An inverse problem for eddy current equations Alberto Valli Department of Mathematics, University of Trento, Italy In collaboration with: Ana Alonso Rodríguez, University of Trento, Italy Jessika Camaño,
More informationWaves. Daniel S. Weile. ELEG 648 Waves. Department of Electrical and Computer Engineering University of Delaware. Plane Waves Reflection of Waves
Waves Daniel S. Weile Department of Electrical and Computer Engineering University of Delaware ELEG 648 Waves Outline Outline Introduction Let s start by introducing simple solutions to Maxwell s equations
More informationECE 6340 Intermediate EM Waves. Fall Prof. David R. Jackson Dept. of ECE. Notes 6
ECE 6340 Intermediate EM Waves Fall 016 Prof. David R. Jackson Dept. of ECE Notes 6 1 Power Dissipated by Current Work given to a collection of electric charges movg an electric field: ( qe ) ( ρ S E )
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 informationVerification and Validation Calculations Using the STAR Geophysical Postprocessor Suite
SAIC-03/1040 Final Report Verification and Validation Calculations Using the STAR Geophysical Postprocessor Suite Prepared by: J. W. Pritchett December 2003 Verification and Validation Calculations Using
More informationElectromagnetic Waves For fast-varying phenomena, the displacement current cannot be neglected, and the full set of Maxwell s equations must be used
Electromagnetic Waves For fast-varying phenomena, the displacement current cannot be neglected, and the full set of Maxwell s equations must be used B( t) E = dt D t H = J+ t D =ρ B = 0 D=εE B=µ H () F
More information20 Poynting theorem and monochromatic waves
0 Poynting theorem and monochromatic waves The magnitude of Poynting vector S = E H represents the amount of power transported often called energy flux byelectromagneticfieldse and H over a unit area transverse
More informationEECS 117. Lecture 25: Field Theory of T-Lines and Waveguides. Prof. Niknejad. University of California, Berkeley
EECS 117 Lecture 25: Field Theory of T-Lines and Waveguides Prof. Niknejad University of California, Berkeley University of California, Berkeley EECS 117 Lecture 25 p. 1/2 Waveguides and Transmission Lines
More informationFinite Element Method (FEM)
Finite Element Method (FEM) The finite element method (FEM) is the oldest numerical technique applied to engineering problems. FEM itself is not rigorous, but when combined with integral equation techniques
More informationPlasma waves in the fluid picture I
Plasma waves in the fluid picture I Langmuir oscillations and waves Ion-acoustic waves Debye length Ordinary electromagnetic waves General wave equation General dispersion equation Dielectric response
More informationTIV Contrast Source Inversion of mcsem data
TIV Contrast ource Inversion of mcem data T. Wiik ( ), L. O. Løseth ( ), B. Ursin ( ), K. Hokstad (, ) ( )epartment for Petroleum Technology and Applied Geophysics, Norwegian University of cience and Technology
More informationSEG Houston 2009 International Exposition and Annual Meeting SUMMARY
Simultaneous joint inversion of MT and CSEM data using a multiplicative cost function Aria Abubakar, Maokun Li, Jianguo Liu and Tarek M. Habashy, Schlumberger-Doll Research, Cambridge, MA, USA SUMMARY
More informationElectromagnetic Wave Propagation Lecture 3: Plane waves in isotropic and bianisotropic media
Electromagnetic Wave Propagation Lecture 3: Plane waves in isotropic and bianisotropic media Daniel Sjöberg Department of Electrical and Information Technology September 2016 Outline 1 Plane waves in lossless
More informationAPPENDIX A: Magnetotelluric Data in Relation to San Pedro Mesa Structural. The San Pedro Mesa structural high (discussed in main text of paper) was
Page of DR for GSA Special Paper 9, Chapter, Geophysical constraints APPENDIX A: Magnetotelluric Data in Relation to San Pedro Mesa Structural High The San Pedro Mesa structural high (discussed in main
More informationModel Validation for Geomagnetically Induced Currents Based on Real Data
for Geomagnetically Induced Currents Based on Real Data Maryam Kazerooni Advisor: Prof. Thomas J. Overbye Collaborator: Prof. Hao Zhu Universities of Illinois at Urbana-Champaign May 214 Maryam Kazerooni
More informationModeling of Overhead Power Lines for Broadband PLC Applications.
Modeling of Overhead Power Lines for Broadband PLC Applications. T. A. Papadopoulos, G. K. Papagiannis, D. P. Labridis Power Systems Laboratory Dept. of Electrical & Computer Engineering Aristotle University
More informationAnna Avdeeva Dmitry Avdeev and Marion Jegen. 31 March Introduction to 3D MT inversion code x3di
Anna Avdeeva (aavdeeva@ifm-geomar.de), Dmitry Avdeev and Marion Jegen 31 March 211 Outline 1 Essential Parts of 3D MT Inversion Code 2 Salt Dome Overhang Study with 3 Outline 1 Essential Parts of 3D MT
More informationFinal Ph.D. Progress Report. Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo
Final Ph.D. Progress Report Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo Dissertation Committee: I. Babuska, L. Demkowicz, C. Torres-Verdin, R. Van
More informationElectromagnetic Waves Across Interfaces
Lecture 1: Foundations of Optics Outline 1 Electromagnetic Waves 2 Material Properties 3 Electromagnetic Waves Across Interfaces 4 Fresnel Equations 5 Brewster Angle 6 Total Internal Reflection Christoph
More informationGeophysics Course Introduction to DC Resistivity
NORAD supported project in MRRD covering Capacity Building and Institutional Cooperation in the field of Hydrogeology for Faryab Province Afghanistan Geophysics Course Introduction to DC Resistivity By
More informationQBX and the DPIE for the Maxwell Equations
University of Illinois @ Urbana-Champaign Fall 2017 - CS 598 APK The Goal For this project, the goals were to implement a model for the Python package pytential for tackling the Maxwell Equations for perfect
More informationECE Spring Prof. David R. Jackson ECE Dept. Notes 1
ECE 6341 Spring 16 Prof. David R. Jackson ECE Dept. Notes 1 1 Fields in a Source-Free Region Sources Source-free homogeneous region ( ε, µ ) ( EH, ) Note: For a lossy region, we replace ε ε c ( / ) εc
More informationPolynomial Chaos Approach for Maxwell s Equations in Dispersive Media
Polynomial Chaos Approach for Maxwell s Equations in Dispersive Media Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 15, 2013 Prof. Gibson (OSU) PC-FDTD
More informationSimulation of marine controlled source electromagnetic measurements using a parallel fourier hp-finite element method
Comput Geosci (2011 15:53 67 DOI.07/s596-0-9195-1 ORIGINAL PAPER Simulation of marine controlled source electromagnetic measurements using a parallel fourier hp-finite element method David Pardo Myung
More informationDielectric Slab Waveguide
Chapter Dielectric Slab Waveguide We will start off examining the waveguide properties of a slab of dielectric shown in Fig... d n n x z n Figure.: Cross-sectional view of a slab waveguide. { n, x < d/
More informationMulti-order Vector Finite Element Modeling of 3D Magnetotelluric Data Including Complex Geometry and Anisotropy. Aixa M.
Multi-order Vector Finite Element Modeling of 3D Magnetotelluric Data Including Complex Geometry and Anisotropy Aixa M. Rivera Ríos Thesis submitted for the degree of Doctor of Philosophy in The University
More informationSIMULATION AND INTERPRETATION OF BOREHOLE GEOPHYSICAL MEASUREMENTS USING hp FINTE ELEMENTS
SIMULATION AND INTERPRETATION OF BOREHOLE GEOPHYSICAL MEASUREMENTS USING hp FINTE ELEMENTS hp-fem team: D. Pardo, M. J. Nam, L. Demkowicz, C. Torres-Verdín, V. M. Calo, M. Paszynski, and P. J. Matuszyk
More informationSimulation study of EM radiation f rom Langmuir waves in warm ma gnetized plasmas
Simulation study of EM radiation f rom Langmuir waves in warm ma gnetized plasmas Iver H. Cairns Eun-Hwa Kim Peter A. Robinson School of Physics, University of Sydney, Australia < AIMS > Demonstrate that
More informationViscous capillary fluids in fast rotation
Viscous capillary fluids in fast rotation Centro di Ricerca Matematica Ennio De Giorgi SCUOLA NORMALE SUPERIORE BCAM BASQUE CENTER FOR APPLIED MATHEMATICS BCAM Scientific Seminar Bilbao May 19, 2015 Contents
More informationA resolution comparison of horizontal and vertical magnetic transfer functions
Journal of the Earth and Space Physics, Vol. 40, No. 3, 2014, P. 47-53 A resolution comparison of horizontal and vertical magnetic transfer functions Habibian Dehkordi, B. 1 and Oskooi, B. 2 1 Assitant
More informationElectromagnetically Induced Flows in Water
Electromagnetically Induced Flows in Water Michiel de Reus 8 maart 213 () Electromagnetically Induced Flows 1 / 56 Outline 1 Introduction 2 Maxwell equations Complex Maxwell equations 3 Gaussian sources
More informationModal Interactions in Lossy Dielectric Metamaterial Slabs
Modal Interactions in Lossy Dielectric Metamaterial Slabs A. B. Yakovlev (), G. Lovat (), P. Burghignoli (), and G. W. Hanson () () University of Mississippi () La Sapienza University of Rome () University
More informationChap. 1 Fundamental Concepts
NE 2 Chap. 1 Fundamental Concepts Important Laws in Electromagnetics Coulomb s Law (1785) Gauss s Law (1839) Ampere s Law (1827) Ohm s Law (1827) Kirchhoff s Law (1845) Biot-Savart Law (1820) Faradays
More informationMagnetotelluric Array and Magnetic Magnetotelluric Surveys
Magnetotelluric Array and Magnetic Magnetotelluric Surveys Dr. Wen J Whan JTech Co., Ltd.; wwj@jtech.com.tw Abstract. Surveys using synchronous MT array could further reduce the local noise. Multi variance
More informationBack to basics : Maxwell equations & propagation equations
The step index planar waveguide Back to basics : Maxwell equations & propagation equations Maxwell equations Propagation medium : Notations : linear Real fields : isotropic Real inductions : non conducting
More informationTypical anisotropies introduced by geometry (not everything is spherically symmetric) temperature gradients magnetic fields electrical fields
Lecture 6: Polarimetry 1 Outline 1 Polarized Light in the Universe 2 Fundamentals of Polarized Light 3 Descriptions of Polarized Light Polarized Light in the Universe Polarization indicates anisotropy
More informationDownloaded 08/29/13 to Redistribution subject to SEG license or copyright; see Terms of Use at
New approach to 3D inversion of MCSEM and MMT data using multinary model transform Alexander V. Gribenko and Michael S. Zhdanov, University of Utah and TechnoImaging SUMMARY Marine controlled-source electromagnetic
More informationSIMULATION OF RESISTIVITY LOGGING-WHILE-DRILLING (LWD) MEASUREMENTS USING A SELF-ADAPTIVE GOAL-ORIENTED HP FINITE ELEMENT METHOD.
SIMULATION OF RESISTIVITY LOGGING-WHILE-DRILLING (LWD) MEASUREMENTS USING A SELF-ADAPTIVE GOAL-ORIENTED HP FINITE ELEMENT METHOD. D. PARDO,, L. DEMKOWICZ, C. TORRES-VERDÍN, AND M. PASZYNSKI, Abstract.
More informationGeneral review: - a) Dot Product
General review: - a) Dot Product If θ is the angle between the vectors a and b, then a b = a b cos θ NOTE: Two vectors a and b are orthogonal, if and only if a b = 0. Properties of the Dot Product If a,
More informationHelmholtz Wave Equation TE, TM, and TEM Modes Rect Rectangular Waveguide TE, TM, and TEM Modes Cyl Cylindrical Waveguide.
Waveguides S. R. Zinka zinka@vit.ac.in School of Electronics Engineering Vellore Institute of Technology April 26, 2013 Outline 1 Helmholtz Wave Equation 2 TE, TM, and TEM Modes Rect 3 Rectangular Waveguide
More informationApproximation SEM-DG pour les problèmes d ondes elasto-acoustiques
Approximation SEM-DG pour les problèmes d ondes elasto-acoustiques Helene Barucq, Henri Calandra, Aurélien Citrain, Julien Diaz, Christian Gout To cite this version: Helene Barucq, Henri Calandra, Aurélien
More informationOverview in Images. S. Lin et al, Nature, vol. 394, p , (1998) T.Thio et al., Optics Letters 26, (2001).
Overview in Images 5 nm K.S. Min et al. PhD Thesis K.V. Vahala et al, Phys. Rev. Lett, 85, p.74 (000) J. D. Joannopoulos, et al, Nature, vol.386, p.143-9 (1997) T.Thio et al., Optics Letters 6, 197-1974
More informationFourier Series Expansion in a Non-Orthogonal System of Coordinates for the Simulation of 3D AC Borehole Resistivity Measurements
Fourier Series Expansion in a Non-Orthogonal System of Coordinates for the Simulation of 3D AC Borehole Resistivity Measurements D. Pardo[a], C. Torres-Verdín[a], M. J. Nam[a], M. Paszynski[b], and V.
More informationReduced MHD. Nick Murphy. Harvard-Smithsonian Center for Astrophysics. Astronomy 253: Plasma Astrophysics. February 19, 2014
Reduced MHD Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 19, 2014 These lecture notes are largely based on Lectures in Magnetohydrodynamics by Dalton
More informationFEM: Domain Decomposition and Homogenization for Maxwell s Equations of Large Scale Problems
FEM: and Homogenization for Maxwell s Equations of Large Scale Problems Karl Hollaus Vienna University of Technology, Austria Institute for Analysis and Scientific Computing February 13, 2012 Outline 1
More informationELE 3310 Tutorial 10. Maxwell s Equations & Plane Waves
ELE 3310 Tutorial 10 Mawell s Equations & Plane Waves Mawell s Equations Differential Form Integral Form Faraday s law Ampere s law Gauss s law No isolated magnetic charge E H D B B D J + ρ 0 C C E r dl
More information22 Phasor form of Maxwell s equations and damped waves in conducting media
22 Phasor form of Maxwell s equations and damped waves in conducting media When the fields and the sources in Maxwell s equations are all monochromatic functions of time expressed in terms of their phasors,
More informationPROBLEM 1 (15 points) In a Cartesian coordinate system, assume the magnetic flux density
PROBLEM 1 (15 points) In a Cartesian coordinate system, assume the magnetic flux density varies as ( ) where is a constant, is the unit vector in x direction. a) Sketch the magnetic flux density and the
More informationModelling Electric Fields In Ireland And UK For Space Weather Applications
Modelling Electric Fields In Ireland nd UK For Space Weather pplications Joan Campanyà 1, Peter Gallagher 1, Seán Blake 1, Mark Gibbs 2, David Jackson 2, Ciarán Beggan 3, Gemma S. Richardson 3, Colin Hogg
More informationElectromagnetic Wave Propagation Lecture 2: Uniform plane waves
Electromagnetic Wave Propagation Lecture 2: Uniform plane waves Daniel Sjöberg Department of Electrical and Information Technology March 25, 2010 Outline 1 Plane waves in lossless media General time dependence
More informationIII. Spherical Waves and Radiation
III. Spherical Waves and Radiation Antennas radiate spherical waves into free space Receiving antennas, reciprocity, path gain and path loss Noise as a limit to reception Ray model for antennas above a
More informationSummary of Beam Optics
Summary of Beam Optics Gaussian beams, waves with limited spatial extension perpendicular to propagation direction, Gaussian beam is solution of paraxial Helmholtz equation, Gaussian beam has parabolic
More informationc 2006 Society for Industrial and Applied Mathematics
SIAM J. APPL. MATH. Vol. 66, No. 6, pp. 2085 206 c 2006 Society for Industrial and Applied Mathematics TWO-DIMENSIONAL HIGH-ACCURACY SIMULATION OF RESISTIVITY LOGGING-WHILE-DRILLING (LWD) MEASUREMENTS
More informationElectromagnetic Waves Retarded potentials 2. Energy and the Poynting vector 3. Wave equations for E and B 4. Plane EM waves in free space
Electromagnetic Waves 1 1. Retarded potentials 2. Energy and the Poynting vector 3. Wave equations for E and B 4. Plane EM waves in free space 1 Retarded Potentials For volume charge & current = 1 4πε
More information6.1 Quadratic Expressions, Rectangles, and Squares. 1. What does the word quadratic refer to? 2. What is the general quadratic expression?
Advanced Algebra Chapter 6 - Note Taking Guidelines Complete each Now try problem in your notes and work the problem 6.1 Quadratic Expressions, Rectangles, and Squares 1. What does the word quadratic refer
More informationConsistency analysis of a 1D Finite Volume scheme for barotropic Euler models
Consistency analysis of a 1D Finite Volume scheme for barotropic Euler models F Berthelin 1,3, T Goudon 1,3, and S Mineaud,3 1 Inria, Sophia Antipolis Méditerranée Research Centre, Proect COFFEE Inria,
More informationMagnetotelluric Survey in an Extremely Noisy Environment at the Pohang Low-Enthalpy Geothermal Area, Korea
Proceedings World Geothermal Congress 2005 Antalya, Turkey, 24-29 April 2005 Magnetotelluric Survey in an Extremely Noisy Environment at the Pohang Low-Enthalpy Geothermal Area, Korea Toshihiro Uchida
More informationTheory of Electromagnetic Nondestructive Evaluation
EE 6XX Theory of Electromagnetic NDE: Theoretical Methods for Electromagnetic Nondestructive Evaluation 1915 Scholl Road CNDE Ames IA 50011 Graduate Tutorial Notes 2004 Theory of Electromagnetic Nondestructive
More informationJoel A. Shapiro January 21, 2010
Joel A. shapiro@physics.rutgers.edu January 21, 20 rmation Instructor: Joel Serin 325 5-5500 X 3886, shapiro@physics Book: Jackson: Classical Electrodynamics (3rd Ed.) Web home page: www.physics.rutgers.edu/grad/504
More informationOn Round-off Error for Adaptive Finite Element Methods
Available online at www.sciencedirect.com Procedia Computer Science 9 (2012 ) 1474 1483 International Conference on Computational Science, ICCS 2012 On Round-off Error for Adaptive Finite Element Methods
More informationHandling Nonlinearity by the Polarization Method and the Newton-Raphson Technique
Handling Nonlinearity by the Polarization Method and the Newton-Raphson Technique M. Kuczmann Laboratory of Electromagnetic Fields, Department of Telecommunications, Széchenyi István University, Egyetem
More informationA Multiphysics Framework Using hp-finite Elements for Electromagnetics Applications
Univ. Carlos III of Madrid A Multiphysics Framework Using hp-finite Elements for Electromagnetics Applications D. Pardo, L. E. García-Castillo, M. J. Nam, C. Torres-Verdín Team: D. Pardo, M. J. Nam, V.
More information1-2 Vector and scalar potentials
1-2 Vector and scalar potentials z a f a r q a q J y x f Maxwell's equations μ (1) ε (2) ρ (3) (4) - Derivation of E and H field by source current J Approach 1 : take (1) and plug into (2) μ ε μ μ (5)
More informationTheoretical Foundation of 3D Alfvén Resonances: Time Dependent Solutions
Theoretical Foundation of 3D Alfvén Resonances: Time Dependent Solutions Tom Elsden 1 Andrew Wright 1 1 Dept Maths & Stats, University of St Andrews DAMTP Seminar - 8th May 2017 Outline Introduction Coordinates
More informationDispersive numerical schemes for Schrödinger equations
Dispersive numerical schemes for Schrödinger equations Enrique Zuazua joint work with L. Ignat zuazua@bcamath.org Basque Center for Applied Mathematics (BCAM), Bilbao, Basque Country, Spain IMA Workshop:
More informationReflection/Refraction
Reflection/Refraction Page Reflection/Refraction Boundary Conditions Interfaces between different media imposed special boundary conditions on Maxwell s equations. It is important to understand what restrictions
More informationA Three-Fluid Approach to Model Coupling of Solar Wind-Magnetosphere-Ionosphere- Thermosphere
A Three-Fluid Approach to Model Coupling of Solar Wind-Magnetosphere-Ionosphere- Thermosphere P. Song Center for Atmospheric Research University of Massachusetts Lowell V. M. Vasyliūnas Max-Planck-Institut
More informationAdvanced Engineering Electromagnetics, ECE750 LECTURE 11 THE FDTD METHOD PART III
Advanced Engineering Electromagnetics, ECE750 LECTURE 11 THE FDTD METHOD PART III 1 11. Yee s discrete algorithm Maxwell s equations are discretized using central FDs. We set the magnetic loss equal to
More informationECE 222b Applied Electromagnetics Notes Set 4b
ECE b Applied Electromagnetics Notes Set 4b Instructor: Prof. Vitali Lomain Department of Electrical and Computer Engineering Universit of California, San Diego 1 Uniform Waveguide (1) Wave propagation
More informationG.Escamez 1,2, F.Sirois 3, M.Tousignant 3, V.Lahtinen 4, A.Stenvall 4, B.Ramdane 2, G.Meunier 2, A.Badel 2,C-E Bruzek 1, P.
30/06/2016 1 3-D NUMERICAL MODELLING OF AC LOSSES IN MULTI-FILAMENTARY MGB 2 WIRES G.Escamez 1,2, F.Sirois 3, M.Tousignant 3, V.Lahtinen 4, A.Stenvall 4, B.Ramdane 2, G.Meunier 2, A.Badel 2,C-E Bruzek
More informationBE 3600 BIOMEDICAL INSTRUMENTATION (LAB) - WEEK 2
BE 3600 BIOMEDICAL INSTRUMENTATION (LAB) - WEEK 2 Principles of Biosensors OBJECTIVE: Learn the various modalities for construction of sensors utilizing electrical and optical measurement techniques. EXPERIMENT
More informationFast inversion of 3D Borehole Resistivity Measurements using Model Reduction Techniques based on 1D Semi-Analytical Solutions.
Fast inversion of 3D Borehole Resistivity Measurements using Model Reduction Techniques based on 1D Semi-Analytical Solutions. Aralar Erdozain Supervised by: Hélène Barucq, David Pardo and Victor Péron
More informationElectromagnetic optics!
1 EM theory Electromagnetic optics! EM waves Monochromatic light 2 Electromagnetic optics! Electromagnetic theory of light Electromagnetic waves in dielectric media Monochromatic light References: Fundamentals
More informationMagnetospheric Electric Fields at Mercury
Magnetospheric Electric Fields at Mercury Lars G. Blomberg Space and Plasma Physics School of Electrical Engineering Royal Institute of Technology (KTH) Stockholm MESSENGER BepiColombo Workshop, Boulder,
More informationClass 15 : Electromagnetic Waves
Class 15 : Electromagnetic Waves Wave equations Why do electromagnetic waves arise? What are their properties? How do they transport energy from place to place? Recap (1) In a region of space containing
More informationGUIDED MICROWAVES AND OPTICAL WAVES
Chapter 1 GUIDED MICROWAVES AND OPTICAL WAVES 1.1 Introduction In communication engineering, the carrier frequency has been steadily increasing for the obvious reason that a carrier wave with a higher
More information