Exponentially Convergent Sparse Discretizations and Application to Near Surface Geophysics
|
|
- Estella Barnett
- 6 years ago
- Views:
Transcription
1 Exponentially Convergent Sparse Discretizations and Application to Near Surface Geophysics Murthy N. Guddati North Carolina State University November 9, 017
2 Outline Part 1: Impedance Preserving Discretization Part : Absorbing Boundary Conditions (1-sided DtN map) Joint work with Tassoulas, Druksin, Lim, Zahid, Savadatti, Thirunavukkarasu Part 3: Complex-length FEM for finite domains (-sided DtN map) Joint work with Druskin, Vaziri Astaneh Part 4: Inversion for Near-surface Geophysics Joint work with Vaziri Astaneh
3 Part1. Impedance Preserving Discretization
4 Model Problem 3D wave equation in free space u u u 1 u 0 x y z c t Fourier transform in t, y, z, with u Ue ik yik zit y z U x 0 ku k k y kz is complex valued c Exact solution: ikx ikx U Ae Be, where k is the horizontal wavenumber i( kxk yk zt ) ikx y z. U e u e is a plane/evanescent wave 4
5 Finite Element Solution on a Uniform Grid in x FE discretization of: U ku x 0 Element contribution matrix with uniform element size of h: k elem A B kh h B A 6 3 A B 1 kh 1 h 3 1 kh 1 h 6 Assembly results in the difference equation: BU AU BU 0 j1 j j1 5
6 Changing Mesh Size: Reflections A simple analysis using two uniform meshes with different element sizes (h, H), but the same material What happens when a right propagating wave hits the interface? Exact solution just passes through Finite element solution reflections due to impedance mismatch R Z Z H H Z Z h h Z Z H h : discrete impedance of left domain : discrete impedance of right domain 6
7 Computing Discrete Impedance (Half-space Stiffness) Basic idea: discrete half-space + finite element = discrete half-space A B U ZU A Z B U 0 A Z B B A Z h U B A Z 0 h 0 h 0 h 1 0 h U1 0 A B 1 kh 1 h 3 1 kh 1 h 6 Zh A B ik 1 kh 1 Error term Z h. depends on element size, resulting impedance mismatch when the element size changes, resulting in reflections 7
8 Optimal Integration for Minimizing Reflection Error Minimize the error in impedance by using generalized integration rules A B h h 4 k h k h Minimize the error term by choosing kh Zh A B ik Error term The error in impedance is completely eliminated! No more reflections Formally valid for more general nd order equations (anisotropic, viscoelasticity etc., electromagnetics etc. G, 006, CMAME) Linear elements + midpoint integration = Impedance Preserving Discretization 8
9 Part. Absorbing Boundary Conditions Perfectly Matched Discrete Layers
10 Perfectly Matched Discrete Layers Impedance Preserving Discretization of PML Perfectly Matched Layers (PML) (Berenger, 1994; Chew et.al. 1995) Step 1: Bend the domain into complex space Reduced reflection into the interior Interior (real x) PML Region (imaginary or complex x) RPML ikl e Step : discretize PMDL domain (in complex space) Impedance is no longer preserved; perfect matching is destroyed Requires a large number of carefully chosen PML layers Impedance preserving discretization comes to the rescue! Impedance is preserved/matched, irrespective of element length, small, large, real, complex Perfectly Matching Discrete Layers (PMDL) Discretize with 3-5 complex-length linear finite elements No discretization error, but truncation causes reflections. The reflection coefficient is derived as R PMDL jnlayer 1 iklj / j1 1 iklj / 10
11 PMDL vs PML: Effectiveness of Midpoint Integration PMDL with 3 layers PML with 3 layers 11
12 PMDL: Some More Old Results Impedance preservation property is valid for any equation that is linear and second order in space (G, CMAME, 006) Elastic and other complicated wave equations (G, Lim & Zahid, 007) PMDL with 5 layers PML with 5 layers Evanescent waves can be treated effectively Padded PMDL contains large real lengths with midpoint integration (Zahid & G, CMAME, 006) 1
13 Salient Features of PMDL Exponential convergence R jnlayer k k j j1 k k j Near optimal discretization Optimal: need staggered grids (with Druskin et al., 003) Links PML to rational ABCs Lindman, Engquist-Majda, Higdon and variants (e.g. CRBC) We started this from E-M/Higdon ABCs (G, Tassoulas, 000) Extensions to corners is straightforward Additional advantage: Provides solutions to some difficult cases Backpropagating waves: anisotropy PML for discrete/periodic media 13
14 PMDL for Backpropagating Waves Opposing signs of phase and group velocities Backpropagating waves grow in the PML region PML cannot work! (Bécache, Fauqueux and Joly, 003) Reduced Reflection into the interior PML result: radiation in anisotropic media Interior PML Region A counter-intuitive idea: make the reflections in PML region decay faster than the growth of the incident wave Works only with PMDL: needs impedance preserving discretization! Result from PMDL after the fix Savadatti & G (01), J Comp. Phys. 14
15 Anisotropic elasticity Tilted Elliptic Case Arbitrary parameters Ideal Slowness Stable parameters Savadatti & G (01), J Comp. Phys. 15
16 Anisotropic elasticity Non-elliptic Case Traditional mesh Two different coordinated materials Savadatti & G (01), J Comp. Phys. 16
17 PMDL for Periodic Media (after discretization) Periodic media has internal reflections and transmissions Constructive interference leads to long-range propagation PML s complex stretching spoils this balance and internal reflections and transmissions get mixed up! Basic Ideas (Discrete/Periodic PMDL): Periodic media = Discrete vector wave equation (vector size = ndof in a cell) Discrete vector equation = impedance preserving discretization of more complicated wave equation Apply PMDL on the complicated wave equation results in impedance matching for periodic media Open problem: stability for complex problems PML for Lattice Waves: 7% reflections w/0 PML layers Discrete PMDL: less than 1% error w/ 4 PMDL layers G & Thirunavukkarasu, JCP (009), Waves
18 Part 3. Two-Sided DtN Map Complex-length Finite Element Method
19 Facilitating the Approximation of -Sided DtN Map u Consider the equation: z A B Exact -sided DtN map: Kexact = B A u 0, z (0, L) By definition, exact DtN Map is impedance preserving: A B Z exact Consider impedance preserving discretization of the interval: A B Kexact =, A B Z B A exact Error in A and B would be similar since: A B Z A B exact Approximating two-sided map reduces to approximating one-sided map Better derivation based on Crank-Nicolson discretization of the propagator 19
20 exp(+il) exp(-il) 1D Helmholtz Equation z 0 u z u 0 z 1 st Order Form v u / z L u 0 i iu z v / i i 0 v / i Eigenvalues i z Downward waves: 0
21 Propagator Approximation approximate P n L 1 u 0 L... u 1 1 u1 u0 1 0 i u1 u0 L v v i 0 v v st Order Form v u / z, v v / i Crank- Nicolson L n ul u 0 i u z v i 0 v Crank- Nicolson Valid for both Downward AND upward waves! 1
22 Padé Approximant Imaginary L 1 L... j exp( ) Padé j d L d P d j j d 0 0 L n ( j 0,, n) Real Complex-Length Finite Element Method n j0 ( n j)! ( ) j x 0 Lj L / x j! ( n j)! j
23 Complex-Length FEM: Exponential Convergence f 1 L 1 Laplace Equation Helmholtz Equation (can t beat Nyquist limit)
24 Complex-Length FEM: Some Observations Exponentially convergent Piecewise linear interpolation sparse computation Edges do not move ( L j =L) can be combined with other types of meshes for other subdomains Mesh is not bent outside (Re(L j )>0) Order of elements do not matter! more on this later. With refinement, and proper ordering, mesh converges to a smooth curve on the complex plane n j0 j exp( ) Padé j d L d P d j j d 0 0 ( n j)! ( ) j x 0 Lj L / x j! ( n j)! j
25 Energy Conservation and Eigenvalue Problems Do complex lengths lead to energy absorption, like PML? No, due to conjugate pair of lengths decay grows back! -sided DtN Map is Hermitian Do complex lengths lead to complex eigenvalues of K with respect to M? No. Eigenvalues are real and positive! Eigenvectors are complex (K and M are complex symmetric) u z u 0 ( ) K M 0
26 Element Ordering u c exp( i x) 6
27 Part 4. Near Surface Geophysical Site Characterization using Guided Wave Inversion
28 Guided Wave Dispersion Low Frequency High Frequency Dispersion Curve 8
29 Spectral Analysis FFT Phase Velocity (m/s) = Frequency (1/s) Wavenumber (1/m) Experimental Dispersion Curve 9
30 Medium Characterization Inverse Identification Experimental Dispersion Curve Iteratively Minimize N experimantal predicted i i E c c i1 Initial Guess Optimization Scheme Gradient Based, e.g. Newton-like Methods Global Search, e.g. Genetic Algorithm Forward Problem: Predicted Dispersion Curve 30
31 Forward Modeling State of the Art 1 u u T u u 1 Drrr Drzr Dr u Drz Dzz Dz u r r r z z r z r 1 T u T u 1 Dr Dz Du Iu 0 r r z r Discretize z direction Hankel Transform r direction Quadratic Eigenvalue Problem det K 0 k c p k 31
32 Reducing the Problem Size: CFEM+PMDL Complex-Length FEM (Finite Layers) Perfectly Matched discrete Layers (Halfspace) 3
33 Forward Modeling: CFEM vs. FEM Error in dispersion curve 33
34 Inversion: Experimental Dispersion Curve Experimental Dispersion Curve Inverse Identification 40 Geophones 36 Geophones 1 Geophones 5 th Mode Experimental Dispersion Curve 4 th Mode 3 rd Mode nd Mode 1 st (fundamental) Mode 34
35 Challenge No analytical derivative Rough misfit function Minimize N experimantal predicted i i E c c i1 Existing approach: Finite Difference Method (FDM) Expensive: Multiple computations of dispersion curve Slow convergence: Oscillatory gradient 35
36 Proposed Derivative for Experimental Curve Analytical Derivative k L( K m j ) i m ( K k ) j L i R R Alternative Approach for Eigenvector R K 1 P k k k i k experimental 6 ( 10 ) 36
37 Proposed Derivative j 1 1 ( P K )( K / m j )( K P) 1 1 ( P K )( K / ki )( K P) eff. ki m Analytical vs. FDM Gradient 37
38 Inversion Results: Synthetic Examples Analytical Gradient FDM Gradient 38
39 Inversion Results: 14-Layer Soil Profile Analytical Gradient FDM Gradient Iterations (Existing) Iterations (Proposed) CPU Time (Proposed) CPU Time (Existing) s s Experimental data from: J Xia et al., J. Environ. Eng. Geophys., 5.3, 1-13 (000) 39
40 Conclusions Discretization that perfectly preserves the impedance is possible Linear FEM with midpoint integration preserves impedance Related to Crank-Nicolson discretization of the propagator Preserves the evanescence in PML region Absorbing Boundary Conds.: Perfectly Matched Discrete Layers (PMDL) Exponential convergence Link to other ABCs we can get the best of both worlds Facilitates stable ABCs for some backward propagating waves Formally extensible to discrete periodic media Open questions: Parameters of discretization for stability and accuracy 40
41 Conclusions (Contd.) Two-sided DtN Map Exponential convergence on the edges is possible with linear interpolation: Complex-length Finite Element Method (CFEM) Impedance preserving discretization is the key! Currently based on Padé approximant; could be further optimized Open questions: further theoretical understanding; extensions to variable coefficients and higher dimensions? Guided Wave Inversion Forward modeling: a good application of CFEM Approximate differentiation of the effective dispersion curve facilitates faster convergence and efficient gradient computation Future work: Bayesian and hybrid inversion 41
42 Thank you! Impedance Preserving Discretization G (006), Arbitrarily wide angle wave equations for complex media, CMAME Perfectly Matched Discrete Layers G, Tassoulas (000), Continued fraction absorbing boundary conditions for the wave equation, J Comp. Acoustics. Asvadurov, Druskin, G, Knizhnerman (003), On optimal finite-difference approximation of PML, SINUM. G, Lim (006), Continued fraction absorbing boundary conditions for convex polygonal domains, IJNME. Thirunavukkarasu, G (011), Absorbing boundary conditions for time-harmonic wave propagation in discretized domains, CMAME. Savadatti, G (01), Accurate absorbing boundary conditions for anisotropic elastic media. parts 1/, JCP. Complex-length FEM G, Druskin, Vaziri Astaneh (016), Exponential Convergence through Linear Finite Element Discretization of Stratified Subdomains, JCP. Vaziri Astaneh, G (016), Efficient Computation of Dispersion Curves for Multilayered Waveguides and Half-Spaces, CMAME. Guided Wave Inversion Vaziri Astaneh, G (016), Improved Algorithms for Inversion of Surface Waves Using Multistation Analysis, Geoph. J. Intl. This work is supported by National Science Foundation DMS , CMMI
CONTINUED-FRACTION ABSORBING BOUNDARY CONDITIONS FOR THE WAVE EQUATION
Journal of Computational Acoustics, Vol. 8, No. 1 (2) 139 156 c IMACS CONTINUED-FRACTION ABSORBING BOUNDARY CONDITIONS FOR THE WAVE EQUATION MURTHY N. GUDDATI Department of Civil Engineering, North Carolina
More informationElectromagnetic wave propagation. ELEC 041-Modeling and design of electromagnetic systems
Electromagnetic wave propagation ELEC 041-Modeling and design of electromagnetic systems EM wave propagation In general, open problems with a computation domain extending (in theory) to infinity not bounded
More informationAn implicit scheme is currently being used with the CFABC. A full explicit scheme is not possible since the mass matrix is singular (the absorbing
Abstract LIM, KENG WIT. Absorbing Boundary Conditions for Corner Regions. (Under the direction of Murthy N. Guddati.) This thesis contains the work that extends the continued fraction absorbing boundary
More informationABSTRACT. THIRUNAVUKKARASU, SENGANAL. Absorbing Boundary Conditions for Molecular Dynamics. (Under the direction of Professor Murthy N. Guddati).
ABSTRACT THIRUNAVUKKARASU, SENGANAL. Absorbing Boundary Conditions for Molecular Dynamics. (Under the direction of Professor Murthy N. Guddati). With the goal of minimizing the domain size for Molecular
More informationA PML absorbing boundary condition for 2D viscoacoustic wave equation in time domain: modeling and imaging
A PML absorbing boundary condition for 2D viscoacoustic wave equation in time domain: modeling and imaging Ali Fathalian and Kris Innanen Department of Geoscience, University of Calgary Summary The constant-q
More informationAn Introduction to Radiation Boundary Conditions for Time-Domain Scattering Problems
An Introduction to Radiation Boundary Conditions for Time-Domain Scattering Problems Tom Hagstrom Southern Methodist University Page 1 of 42 Contributors to aspects of this work: T. Warburton, Rice, D.
More informationA new second order absorbing boundary layer formulation for anisotropic-elastic wavefeld simulation
A new second order absorbing boundary layer formulation for anisotropic-elastic wavefeld simulation Junxiao Li, Kris Innanen and Bing Wang University of Calgary China University of Petroleum-Beijing Summary
More informationUnphysical negative values of the anelastic SH plane wave energybased transmission coefficient
Shahin Moradi and Edward S. Krebes Anelastic energy-based transmission coefficient Unphysical negative values of the anelastic SH plane wave energybased transmission coefficient ABSTRACT Computing reflection
More informationThe Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator
The Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator Martin J. Gander and Achim Schädle Mathematics Section, University of Geneva, CH-, Geneva, Switzerland, Martin.gander@unige.ch
More informationA PRECONDITIONER FOR THE HELMHOLTZ EQUATION WITH PERFECTLY MATCHED LAYER
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 A PRECONDITIONER FOR THE HELMHOLTZ EQUATION WITH PERFECTLY
More informationComplete Radiation Boundary Conditions, Optimal Grids, and Hermite Methods
Complete Radiation Boundary Conditions, Optimal Grids, and Hermite Methods Tom Hagstrom Southern Methodist University Contributors to aspects of this work: CRBC: T. Warburton, Virginia Tech M. de Castro
More informationWe G Model Reduction Approaches for Solution of Wave Equations for Multiple Frequencies
We G15 5 Moel Reuction Approaches for Solution of Wave Equations for Multiple Frequencies M.Y. Zaslavsky (Schlumberger-Doll Research Center), R.F. Remis* (Delft University) & V.L. Druskin (Schlumberger-Doll
More informationSolving Time-Harmonic Scattering Problems by the Ultra Weak Variational Formulation
Introduction Solving Time-Harmonic Scattering Problems by the Ultra Weak Variational Formulation Plane waves as basis functions Peter Monk 1 Tomi Huttunen 2 1 Department of Mathematical Sciences University
More informationSome Problems in the Numerical Analysis of Elastic Waves
Some Problems in the Numerical Analysis of Elastic Waves Tom Hagstrom Southern Methodist University Collaborators: Daniel Appelö, University of Colorado Daniel Baffet, Basel Jeff Banks, RPI Jacobo Bielak,
More informationAn acoustic wave equation for orthorhombic anisotropy
Stanford Exploration Project, Report 98, August 1, 1998, pages 6?? An acoustic wave equation for orthorhombic anisotropy Tariq Alkhalifah 1 keywords: Anisotropy, finite difference, modeling ABSTRACT Using
More informationMultiscale methods for time-harmonic acoustic and elastic wave propagation
Multiscale methods for time-harmonic acoustic and elastic wave propagation Dietmar Gallistl (joint work with D. Brown and D. Peterseim) Institut für Angewandte und Numerische Mathematik Karlsruher Institut
More informationBeam Propagation Method Solution to the Seminar Tasks
Beam Propagation Method Solution to the Seminar Tasks Matthias Zilk The task was to implement a 1D beam propagation method (BPM) that solves the equation z v(xz) = i 2 [ 2k x 2 + (x) k 2 ik2 v(x, z) =
More informationPerfectly Matched Layer Finite Element Simulation of Parasitic Acoustic Wave Radiation in Microacoustic Devices
Perfectly Matched Layer Finite Element Simulation of Parasitic Acoustic Wave Radiation in Microacoustic Devices Markus Mayer, Sabine Zaglmayr, Karl Wagner and Joachim Schöberl EPCOS AG, Surface Acoustic
More informationFINITE-DIFFERENCE FREQUENCY-DOMAIN ANALYSIS OF NOVEL PHOTONIC
FINITE-DIFFERENCE FREQUENCY-DOMAIN ANALYSIS OF NOVEL PHOTONIC WAVEGUIDES Chin-ping Yu (1) and Hung-chun Chang (2) (1) Graduate Institute of Electro-Optical Engineering, National Taiwan University, Taipei,
More informationInfinite boundary element absorbing boundary for wave propagation simulations
GEOPHYSICS, VOL. 65, NO. 2 (MARCH-APRIL 2000); P. 596 602, 10 FIGS. Infinite boundary element absorbing boundary for wave propagation simulations Li-Yun Fu and Ru-Shan Wu ABSTRACT In the boundary element
More informationAdditive Schwarz method for scattering problems using the PML method at interfaces
Additive Schwarz method for scattering problems using the PML method at interfaces Achim Schädle 1 and Lin Zschiedrich 2 1 Zuse Institute, Takustr. 7, 14195 Berlin, Germany schaedle@zib.de 2 Zuse Institute,
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationNumerical analysis of ultrasonic guided waves propagation in highly. attenuative viscoelastic material. Li Hong, Wang Qingfeng
Proceedings of the 8th International Conference on Sensing Technology, Sep. -4, 4, Liverpool, UK Numerical analysis of ultrasonic guided waves propagation in highly attenuative viscoelastic material Li
More informationIntroduction to PML in time domain
Introduction to PML in time domain Alexander Thomann Introduction to PML in time domain - Alexander Thomann p.1 Overview 1 Introduction 2 PML in one dimension Classical absorbing layers One-dimensional
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationAn optimized perfectly matched layer for the Schrödinger equation
An optimized perfectly matched layer for the Schrödinger equation Anna Nissen Gunilla Kreiss June 6, 009 Abstract A perfectly matched layer (PML) for the Schrödinger equation using a modal ansatz is presented.
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationFull Wave Analysis of RF Signal Attenuation in a Lossy Rough Surface Cave Using a High Order Time Domain Vector Finite Element Method
Progress In Electromagnetics Research Symposium 2006, Cambridge, USA, March 26-29 425 Full Wave Analysis of RF Signal Attenuation in a Lossy Rough Surface Cave Using a High Order Time Domain Vector Finite
More informationFinite Element Analysis of Acoustic Scattering
Frank Ihlenburg Finite Element Analysis of Acoustic Scattering With 88 Illustrations Springer Contents Preface vii 1 The Governing Equations of Time-Harmonic Wave Propagation, 1 1.1 Acoustic Waves 1 1.1.1
More informationNear-optimal perfectly matched layers for indefinite Helmholtz problems. Vladimir Druskin, Stefan Güttel and Leonid Knizhnerman.
Near-optimal perfectly matched layers for indefinite Helmholtz problems Vladimir Druskin, Stefan Güttel and Leonid Knizhnerman November 203 MIMS EPrint: 203.53 Manchester Institute for Mathematical Sciences
More informationA Quasi-Optimal Non-Overlapping Domain Decomposition Algorithm for the Helmholtz Equation
A Quasi-Optimal Non-Overlapping Domain Decomposition Algorithm for the Helmholtz Equation Y. Boubendir, X. Antoine, C. Geuzaine June 28, 2011 Abstract This paper presents a new non-overlapping domain decomposition
More informationA Quasi-Optimal Non-Overlapping Domain Decomposition Algorithm for the Helmholtz Equation
A Quasi-Optimal Non-Overlapping Domain Decomposition Algorithm for the Helmholtz Equation Y. Boubendir, X. Antoine, C. Geuzaine May 8, 2012 Abstract This paper presents a new non-overlapping domain decomposition
More informationComparison of a Finite Difference and a Mixed Finite Element Formulation of the Uniaxial Perfectly Matched Layer
Comparison of a Finite Difference and a Mixed Finite Element Formulation of the Uniaxial Perfectly Matched Layer V. A. Bokil a and M. W. Buksas b Center for Research in Scientific Computation a North Carolina
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 informationPseudospectral and High-Order Time-Domain Forward Solvers
Pseudospectral and High-Order Time-Domain Forward Solvers Qing H. Liu G. Zhao, T. Xiao, Y. Zeng Department of Electrical and Computer Engineering Duke University DARPA/ARO MURI Review, August 15, 2003
More informationLecture 16: Relaxation methods
Lecture 16: Relaxation methods Clever technique which begins with a first guess of the trajectory across the entire interval Break the interval into M small steps: x 1 =0, x 2,..x M =L Form a grid of points,
More informationAppendix B: The Transfer Matrix Method
Y D Chong (218) PH441: Quantum Mechanics III Appendix B: The Transfer Matrix Method The transfer matrix method is a numerical method for solving the 1D Schrödinger equation, and other similar equations
More informationD034 On the Effects of Anisotropy in Marine CSEM
D034 On the Effects of Anisotropy in Marine CSEM L.O. Løseth* (Statoil ASA, formerly Dept. of Physics, NTNU), B. Ursin (Dept. of Appl. Geophys., NTNU) & L. Amundsen (Statoil ASA) SUMMARY In marine CSEM
More informationNonstandard Finite Difference Time Domain Algorithm for Berenger s Perfectly Matched Layer
ACES JOURNAL, VOL. 6, NO., FEBRUARY 011 153 Nonstandard Finite Difference Time Domain Algorithm for Berenger s Perfectly Matched Layer Naoki Okada and James B. Cole Graduate School of Systems and Information
More informationA Fast Fourier transform based direct solver for the Helmholtz problem. AANMPDE 2017 Palaiochora, Crete, Oct 2, 2017
A Fast Fourier transform based direct solver for the Helmholtz problem Jari Toivanen Monika Wolfmayr AANMPDE 2017 Palaiochora, Crete, Oct 2, 2017 Content 1 Model problem 2 Discretization 3 Fast solver
More informationSource Free Surface x
Finite-dierence time-domain model for elastic waves in the ground Christoph T. Schroeder and Waymond R. Scott, Jr. School of Electrical and Computer Engineering Georgia Institute of Technology Atlanta,
More informationarxiv: v1 [math.na] 19 Nov 2018
A Bivariate Spline Solution to the Exterior Helmholtz Equation and its Applications Shelvean Kapita Ming-Jun Lai arxiv:1811.07833v1 [math.na] 19 Nov 2018 November 20, 2018 Abstract We explain how to use
More informationA new method for the solution of scattering problems
A new method for the solution of scattering problems Thorsten Hohage, Frank Schmidt and Lin Zschiedrich Konrad-Zuse-Zentrum Berlin, hohage@zibde * after February 22: University of Göttingen Abstract We
More informationModelling in photonic crystal structures
Modelling in photonic crystal structures Kersten Schmidt MATHEON Nachwuchsgruppe Multiscale Modelling and Scientific Computing with PDEs in collaboration with Dirk Klindworth (MATHEON, TU Berlin) Holger
More informationSurface Waves and Free Oscillations. Surface Waves and Free Oscillations
Surface waves in in an an elastic half spaces: Rayleigh waves -Potentials - Free surface boundary conditions - Solutions propagating along the surface, decaying with depth - Lamb s problem Surface waves
More informationContents. 1 Basic Equations 1. Acknowledgment. 1.1 The Maxwell Equations Constitutive Relations 11
Preface Foreword Acknowledgment xvi xviii xix 1 Basic Equations 1 1.1 The Maxwell Equations 1 1.1.1 Boundary Conditions at Interfaces 4 1.1.2 Energy Conservation and Poynting s Theorem 9 1.2 Constitutive
More informationEffects of tool eccentricity on wave dispersion properties in borehole acoustic logging while drilling
Effects of tool eccentricity on wave dispersion properties in borehole acoustic logging while drilling Yibing Zheng and M. Nafi Toksöz Earth Resources Laboratory Dept. of Earth, Atmospheric, and Planetary
More informationEfficient outdoor sound propagation modelling in time-domain
Efficient outdoor sound propagation modelling in time-domain Timothy Van Renterghem Acoustics Group Department of Information Technology Ghent University/iMinds CEAS, X-Noise : Atm. and ground effects
More informationSpatio-Temporal Characterization of Bio-acoustic Scatterers in Complex Media
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Spatio-Temporal Characterization of Bio-acoustic Scatterers in Complex Media Karim G. Sabra, School of Mechanical Engineering,
More informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
More informationOPAC102. The Acoustic Wave Equation
OPAC102 The Acoustic Wave Equation Acoustic waves in fluid Acoustic waves constitute one kind of pressure fluctuation that can exist in a compressible fluid. The restoring forces responsible for propagating
More informationModel reduction of wave propagation
Model reduction of wave propagation via phase-preconditioned rational Krylov subspaces Delft University of Technology and Schlumberger V. Druskin, R. Remis, M. Zaslavsky, Jörn Zimmerling November 8, 2017
More informationUsing windscreens to improve the efficiency of noise barriers in wind: finite-difference time-domain simulations
Using windscreens to improve the efficiency of noise barriers in wind: finite-difference time-domain simulations PACS REFERENCE: 43.8.Js, 43.8.Gq, 43.5.Gf Van Renterghem Timothy; Botteldooren Dick University
More informationApplication of Nodal Discontinuous Glaerkin Methods in Acoustic Wave Modeling
1 Application of Nodal Discontinuous Glaerkin Methods in Acoustic Wave Modeling Xin Wang ABSTRACT This work will explore the discontinuous Galerkin finite element method (DG-FEM) for solving acoustic wave
More informationCHAPTER 9 ELECTROMAGNETIC WAVES
CHAPTER 9 ELECTROMAGNETIC WAVES Outlines 1. Waves in one dimension 2. Electromagnetic Waves in Vacuum 3. Electromagnetic waves in Matter 4. Absorption and Dispersion 5. Guided Waves 2 Skip 9.1.1 and 9.1.2
More informationShifted Laplace and related preconditioning for the Helmholtz equation
Shifted Laplace and related preconditioning for the Helmholtz equation Ivan Graham and Euan Spence (Bath, UK) Collaborations with: Paul Childs (Schlumberger Gould Research), Martin Gander (Geneva) Douglas
More informationdirection normal to the interface. As such, the PML medium has been touted as the best absorbing boundary condition for numerical solution of partial
Perfectly Matched Layers in the Discretized Space: An Analysis and Optimization y W. C. Chew and J. M. Jin Center for Computational Electromagnetics Electromagnetics Laboratory Department of Electrical
More informationOn domain decomposition preconditioners for finite element approximations of the Helmholtz equation using absorption
On domain decomposition preconditioners for finite element approximations of the Helmholtz equation using absorption Ivan Graham and Euan Spence (Bath, UK) Collaborations with: Paul Childs (Emerson Roxar,
More informationSpectral Processing. Misha Kazhdan
Spectral Processing Misha Kazhdan [Taubin, 1995] A Signal Processing Approach to Fair Surface Design [Desbrun, et al., 1999] Implicit Fairing of Arbitrary Meshes [Vallet and Levy, 2008] Spectral Geometry
More information2072 J. Acoust. Soc. Am. 102 (4), October /97/102(4)/2072/11/$ Acoustical Society of America 2072
The perfectly matched layer for acoustic waves in absorptive media Qing-Huo Liu and Jianping Tao Klipsch School of Electrical and Computer Engineering, New Mexico State University, Las Cruces, New Mexico
More informationRiesz bases of Floquet modes in semi-infinite periodic waveguides and implications
Riesz bases of Floquet modes in semi-infinite periodic waveguides and implications Thorsten Hohage joint work with Sofiane Soussi Institut für Numerische und Angewandte Mathematik Georg-August-Universität
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More information1.1 A Scattering Experiment
1 Transfer Matrix In this chapter we introduce and discuss a mathematical method for the analysis of the wave propagation in one-dimensional systems. The method uses the transfer matrix and is commonly
More informationLecture 2 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell
Lecture Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell 1. Dispersion Introduction - An electromagnetic wave with an arbitrary wave-shape
More informationThe Pseudospectral Method
The Pseudospectral Method Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Heiner Igel Computational Seismology 1 / 43 Outline 1 Introduction Motivation History
More informationLinearized AVO in viscoelastic media Shahpoor Moradi,Kristopher A. Innanen, University of Calgary, Department of Geoscience, Calgary, Canada
Shahpoor Moradi,Kristopher. Innanen, University of Calgary, Department of Geoscience, Calgary, Canada SUMMRY Study of linearized reflectivity is very important for amplitude versus offset VO) analysis.
More informationUSAGE OF NUMERICAL METHODS FOR ELECTROMAGNETIC SHIELDS OPTIMIZATION
October 4-6, 2007 - Chiinu, Rep.Moldova USAGE OF NUMERICAL METHODS FOR ELECTROMAGNETIC SHIELDS OPTIMIZATION Ionu- P. NICA, Valeriu Gh. DAVID, /tefan URSACHE Gh. Asachi Technical University Iai, Faculty
More informationA. Safaeinili and D. E. Chimenti Center for Nondestructive Evaluation Iowa State University Ames, Ia
FLOQUET ANALYSIS OF LAMB WAVES PROPAGATING IN PERIODICALLY-LAYERED COMPOSITES A. Safaeinili and D. E. Chimenti Center for Nondestructive Evaluation Iowa State University Ames, Ia. 50011 INTRODUCTION In
More informationComputer simulation of multiscale problems
Progress in the SSF project CutFEM, Geometry, and Optimal design Computer simulation of multiscale problems Axel Målqvist and Daniel Elfverson University of Gothenburg and Uppsala University Umeå 2015-05-20
More informationHigh-order non-reflecting boundary conditions for dispersive wave problems in stratified media
High-order non-reflecting boundary conditions for dispersive wave problems in stratified media V. van ~oolen', B. eta' & D. ~ivoli~ '~e~artment of Mathematics, Naval Postgraduate School, USA Department
More informationFinite Difference and Finite Element Methods
Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element
More informationFrozen light in photonic crystals with degenerate band edge
Frozen light in photonic crystals with degenerate band edge Alex Figotin and Ilya Vitebskiy Department of Mathematics, University of California, Irvine, California 92697, USA Received 9 October 2006; published
More informationComparison between least-squares reverse time migration and full-waveform inversion
Comparison between least-squares reverse time migration and full-waveform inversion Lei Yang, Daniel O. Trad and Wenyong Pan Summary The inverse problem in exploration geophysics usually consists of two
More informationWe have the pleasure to inform you that, after the reviewing process, your paper entitled:
ΔΙΕΥΘΥΝΣΗ askarlat@geo.auth.gr ΣΥΓΓΡΑΦΕΑ ΕΠΙΚΟΙΝΩΝΙΑΣ Patras, 1-03-2010 Dear colleague, We have the pleasure to inform you that, after the reviewing process, your paper entitled: IMPLEMENTATION OF A NON-SPLITTING
More informationA perfectly matched layer absorbing boundary condition for the second-order seismic wave equation
Geophys. J. Int. (3) 15, 1 153 A perfectly matched layer absorbing boundary condition for the second-order seismic wave equation Dimitri Komatitsch and Jeroen Tromp Seismological Laboratory, California
More informationApplication of the perfectly matched layer (PML) absorbing boundary condition to elastic wave propagation
Application of the perfectly matched layer (PML) absorbing boundary condition to elastic wave propagation Frank D. Hastings, a) John B. Schneider, and Shira L. Broschat School of Electrical Engineering
More informationDedicated to my father Mr. Ambrose Duru ( )
Dedicated to my father Mr. Ambrose Duru (1939 2001) List of papers This thesis is based on the following papers, which are referred to in the text by their Roman numerals. I II III IV V K. Duru and G.
More informationHigh Frequency Scattering by Convex Polygons Stephen Langdon
Bath, October 28th 2005 1 High Frequency Scattering by Convex Polygons Stephen Langdon University of Reading, UK Joint work with: Simon Chandler-Wilde Steve Arden Funded by: Leverhulme Trust University
More informationConstrained Minimization and Multigrid
Constrained Minimization and Multigrid C. Gräser (FU Berlin), R. Kornhuber (FU Berlin), and O. Sander (FU Berlin) Workshop on PDE Constrained Optimization Hamburg, March 27-29, 2008 Matheon Outline Successive
More informationElastic wavefield separation for VTI media
CWP-598 Elastic wavefield separation for VTI media Jia Yan and Paul Sava Center for Wave Phenomena, Colorado School of Mines ABSTRACT The separation of wave modes from isotropic elastic wavefields is typically
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 informationCHAPTER 4 ELECTROMAGNETIC WAVES IN CYLINDRICAL SYSTEMS
CHAPTER 4 ELECTROMAGNETIC WAVES IN CYLINDRICAL SYSTEMS The vector Helmholtz equations satisfied by the phasor) electric and magnetic fields are where. In low-loss media and for a high frequency, i.e.,
More informationA far-field based T-matrix method for three dimensional acoustic scattering
ANZIAM J. 50 (CTAC2008) pp.c121 C136, 2008 C121 A far-field based T-matrix method for three dimensional acoustic scattering M. Ganesh 1 S. C. Hawkins 2 (Received 14 August 2008; revised 4 October 2008)
More informationTLM method and acoustics
INTERNATIONAL JOURNAL OF NUMERICAL MODELLING: ELECTRONIC NETWORKS, DEVICES AND FIELDS Int. J. Numer. Model. 2001; 14:171}183 (DOI: 10.1002/jnm.405) TLM method and acoustics Jorge A. PortmH * and Juan A.
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 informationMain Menu SUMMARY INTRODUCTION
- a new method of solution Lasse Amundsen, Børge Arntsen, Arne Reitan, Eirik Ø Dischler, and Bjørn Ursin StatoilHydro Reseach Center, and NTNU SUMMARY We present a wave propagation method rigorous in one-way
More informationElectromagnetic fields and waves
Electromagnetic fields and waves Maxwell s rainbow Outline Maxwell s equations Plane waves Pulses and group velocity Polarization of light Transmission and reflection at an interface Macroscopic Maxwell
More informationThe exploding-reflector concept for ground-penetrating-radar modeling
ANNALS OF GEOPHYSCS, VOL. 45, N. 3/4, June/August 2002 The exploding-reflector concept for ground-penetrating-radar modeling José M. Carcione ( 1 ), Laura Piñero Feliciangeli ( 2 ) and Michela Zamparo
More informationStochastic Collocation Methods for Polynomial Chaos: Analysis and Applications
Stochastic Collocation Methods for Polynomial Chaos: Analysis and Applications Dongbin Xiu Department of Mathematics, Purdue University Support: AFOSR FA955-8-1-353 (Computational Math) SF CAREER DMS-64535
More informationAn eigenvalue method using multiple frequency data for inverse scattering problems
An eigenvalue method using multiple frequency data for inverse scattering problems Jiguang Sun Abstract Dirichlet and transmission eigenvalues have important applications in qualitative methods in inverse
More informationincident wave reflected wave region 1 region 2 z=0 plane transmitted wave
A 3-D Perfectly Matched Medium from Modied Maxwell's Equations with Stretched Coordinates y Weng Cho Chew William H. Weedon Electromagnetics Laboratory Department of Electrical Computer Engineering University
More informationIntegration, differentiation, and root finding. Phys 420/580 Lecture 7
Integration, differentiation, and root finding Phys 420/580 Lecture 7 Numerical integration Compute an approximation to the definite integral I = b Find area under the curve in the interval Trapezoid Rule:
More informationImplicit 3-D depth migration by wavefield extrapolation with helical boundary conditions
Stanford Exploration Project, Report 97, July 8, 1998, pages 1 13 Implicit 3-D depth migration by wavefield extrapolation with helical boundary conditions James Rickett, Jon Claerbout, and Sergey Fomel
More informationPerfectly matched layer for second-order time-domain elastic wave equation: formulation and stability
arxiv:32.3722v [physics.comp-ph] 3 Dec 23 Perfectly matched layer for second-order time-domain elastic wave equation: formulation and stability Hisham Assi, Richard S. C. Cobbold Institute of Biomaterials
More informationFastMMLib: a generic Fast Multipole Method library. Eric DARRIGRAND. IRMAR Université de Rennes 1
FastMMLib: a generic Fast Multipole Method library Eric DARRIGRAND joint work with Yvon LAFRANCHE IRMAR Université de Rennes 1 Outline Introduction to FMM Historically SVD principle A first example of
More informationScattering of electromagnetic waves by thin high contrast dielectrics II: asymptotics of the electric field and a method for inversion.
Scattering of electromagnetic waves by thin high contrast dielectrics II: asymptotics of the electric field and a method for inversion. David M. Ambrose Jay Gopalakrishnan Shari Moskow Scott Rome June
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 informationPeriodic Structures in FDTD
EE 5303 Electromagnetic Analsis Using Finite Difference Time Domain Lecture #19 Periodic Structures in FDTD Lecture 19 These notes ma contain coprighted material obtained under fair use rules. Distribution
More informationA Fast Iterative Solver for Scattering by Elastic Objects in Layered Media
A Fast Iterative Solver for Scattering by Elastic Objects in Layered Media K. Ito and J. Toivanen Center for Research in Scientific Computation, North Carolina State University, Raleigh, North Carolina
More informationMulti-Domain Approaches for the Solution of High-Frequency Time-Harmonic Propagation Problems
Académie universitaire Wallonie Europe Université de Liège Faculté des Sciences Appliquées Collège de doctorat en Électricité, électronique et informatique Multi-Domain Approaches for the Solution of High-Frequency
More information