Acoustic Wave Equation
|
|
- Coral Williams
- 5 years ago
- Views:
Transcription
1 Acoustic Wave Equation Sjoerd de Ridder (most of the slides) & Biondo Biondi January 16 th 2011
2 Table of Topics Basic Acoustic Equations Wave Equation Finite Differences Finite Difference Solution Pseudospectral Solution Stability and Accuracy Green s function Perturbation Representation Born Approximation
3 Basic linearized acoustic equations in lossless, isotropic, non flowing media Linearized - Linear for small perturbation on a static state. Lossless - Material parameters are independent of time. Isotropic - Material response independent of direction. Non flowing - No material derivative Equation of motion (three equations for three components) Acoustic stress-strain relationship (a pressure-rate strain-rate relation) ρ t υ i + i p = f i (1) ρ t p + i υ i = q (2)
4 Fields p = p(x, t) υ i = υ i (x, t) pressure i component of velocity Sources q = q(x, t) f i = f i volume injection rate i component of external force Medium Parameters κ = κ(x) ρ = ρ(x) compressibility density
5 Wave Equation Solve equations (1) and (2) for pressure ρ i ρ 1 i p ρκ t 2 p = ρ i ρ 1 f i ρ t q, (3) or i 2 p ρκ t 2 p = ρ i ρ 1 f i ρ t q + ρ 1 i ρ i p. (4) Thus in a constant density and sourceless medium i 2 p c 2 t 2 p = 0, (5) with wave velocity c = c(x) = κρ, κ = κ(x), ρ = ρ 0.
6 Finite Differences Derivation of finite difference stencils for F (s) s Expand F (s + s) in Taylor series F (s + s) = F (s) + 1 1! sf (s) s + i=2 1 i! i sf (s) { s} i (6) Express F (s) s as a function of... s F (s) = 1 s {F (s + s) F (s)} i=2 1 i! i sf (s) { s} i 1 (7) this is a forward finite difference stencil.
7 Expand F (s + s) and F (s s) in Taylor series F (s + s) = F (s) + 1 1! sf (s) s + i=2 1 i! i sf (s) { s} i (8) F (s s) = F (s) 1 1! sf (s) s + i=2 Substract equations (9) from (8), express s F (s) = 1 2 s {F (s + s) F (s s)} this is a centered finite difference stencil. i=1 1 i! i sf (s) { s} i (9) F (s) s as a function of... 1 (1 + 2i)! 1+2i s F (s) { s} 2i (10)
8 or last, s F (s) in a backward finite difference stencil from equation (9) as s F (s) = 1 s {F (s) F (s s)} i=2 1 i! i sf (s) { s} i 1 (11)
9 Derivation of finite difference stencils for 2 s F (s) Add equation (9) and (8), express 2 s F (s) as a function of... s 2 F (s) = 1 {F (s s) 2F (s) + F (s + s)} + ( s) 2 1 (2 + 2i)! 2+2i s F (s) { s} 2i (12) i=1 This is a centered finite difference. Forward and backward finite difference stencils for 2 s F (s) can be obtained from combinations of Taylor series for F (s + s) and F (s + 2 s), or F (s s) and F (s 2 s) respectively.
10 Finite Difference Solution of WE Wave equation, FD 2 nd -order in space h = y = h 2 P(x, t) 1 c 2 (x) 2 t P(x, t) = 1 ( h) P x (t) 1 c 2 (x) 2 t P x (t) (13) Laplacian, FD 4 th -order in space 2 P(x, t) 1 c 2 (x) 2 t P(x, t) = 1 12( h) c 2 (x) 2 t P x (t) (14) P x (t)
11 Laplacian, FD 4 th -order in space, isotropic 2 P(x, t) = 1 α 12( h) P x (t) + 1 β 12( h) P x (t) (15) β = 1 α 2, for example: α = 1 β = 0, α = 1/2 β = 1/4 or α = 2/3 β = 1/6.
12 Wave equation, FD 2 nd -order time stepping 2 t P(x, t) c 2 (x) 2 P(x, t) = 1 t P t (x) c 2 (x) 2 P t (x) Solve for P t+1 (x) (16) P t+1 (x) = 2P t (x) P t 1 (x) + t 2 c 2 (x) 2 P t (x) (17) Wave equation, FD 4 th -order time stepping Include the 4 th -order derivative from equation (12), by substituting the wave equation (Dablain, 1986), as 4 t P(x, t) = 2 t 2 t P(x, t) = 2 t c 2 (x) 2 P(x, t) (18)
13 Pseudospectral (Fourier) methods Laplacian computed using FFTs: c 2 (x) 2 P t (x) c 2 (x) FFT { 1 } k 2 FFT [P t (x)] Wave equation, FD 2 nd -order time stepping and pseudospectral Laplacian: P t+1 (x) = 2P t (x) P t 1 (x) + t 2 c 2 (x) FFT { 1 } k 2 FFT [P t (x)]
14 Stability and accuracy of explicit methods Courant number : Cour = maximum velocity. c max t min( x, y, z) where c max is the Courant-Friedrichs-Lewy (CFL) condition : Cour 1 it is a necessary, but not sufficient condition for a stable explicit extrapolator. Numerical dispersion causes c P c, where c P = ω/ k is the effective phase velocity of numerically propagated waves
15 Stability and accuracy analysis of pseudospectral methods Substitute [ ( a generic )] plane wave solution: exp i kx + ωt Dispersion relation: ω = «2 sin 1 ± c t k 2 t Phase velocity: c P = ω 2 sin 1 k = «± c t k 2 t k For stability it must be c t k 2 1: 1D: Maximum k equal to Nyquist wavenumber k Nyq = π/ x stability requires Cour 2/π D: k max = 2k Nyq stability requires Cour 2/π D: k max = 3k Nyq stability requires Cour 2/ 3π 0.367
16 Stability and accuracy of 2 nd -order in time and space Substitute [ ( a generic )] plane wave solution: exp i kx + ωt Dispersion relation: ω =» 2 sin 1 c t x q sin 2 kx x ( 2 )+sin 2 kz z ( 2 ) t Phase velocity (worst case at k x = 0 or k z = 0): c P = ω k x = 2 sin 1 [ c t kx x sin( x 2 )] tk x For stability the argument of sin 1 must be between -1 and 1: 1D: Cour 1 2D: Worst case at k x = k z = k Nyq : Cour 2/
17 Observations Stability Stability constraint becomes more stringent with higher dimensions FD more stable than pseudospectral because errors in the spatial derivatives slows down high frequencies. Dispersion Pseudospectral High frequencies (wavenumbers) arrive before low frequencies (wavenumbers). Dispersion gets worse as the Courant number increases. FD High frequencies (wavenumbers) tend to arrive after low frequencies (wavenumbers). Dispersion gets better as the Courant number increases.
18 Frequency dispersion with finite-differences schemes
19 Frequency dispersion with finite-differences schemes
20 Frequency dispersion with finite-differences schemes
21 Frequency dispersion with finite-differences schemes
22 Frequency dispersion with pseudspectral Laplacian
23 Wavelength dispersion with finite-differences schemes
24 Wavelength dispersion with finite-differences schemes
25 Wavelength dispersion with finite-differences schemes
26 Wavelength dispersion with finite-differences schemes
27 Wavelength dispersion with pseudspectral Laplacian
28 Green s function Introduce Green s function for a constant density and sourceless medium equation (5) by a point source term acting at t = 0 and x = x s 2 i G c 2 2 t G = δ(x x s )δ(t), (19) where G = G(x, x s, t) is the Green s function. The solution for pressure to another forcing function for example s = sx, t can be represented as p(x, t) = G(x, x, t t )s(x, t )dx dt (20)
29 Perturbation Representation Represent the medium velocity as a background velocity and a perturbation c 2 (x) = c 2 b (x) [1 + α(x)] (21) Substitution into equation (19) gives i 2 G(x, x s, t) c 2 b (x) 2 t G(x, x s, t) = (22) δ(x x s )δ(t) + α(x)c 2 b (x) 2 t G(x, x s, t),
30 Introducing G b (x, x s, t) as a solution to 2 i G b (x, x s, t) c 2 b (x) 2 t G b (x, x s, t) = δ(x x s )δ(t), (23) we see that is we represent the full solution as a sum of the background solution plus a perturbed solution as G(x, x s, t) = G b (x, x s, t) + G p (x, x s, t). (24)
31 Equation (22) can be thus written as 2 i G p (x, x s, t) c 2 b (x) 2 t G p (x, x s, t) = α(x)c 2 b (x) 2 t G(x, x s, t). (25) Note the forcing function dependent on medium parameter α. Thus using a representation as (20) for G p (x, x s, t) we find for G(x, x s, t) G(x, x s, t) = G b (x, x s, t) (26) G b (x, x, t t )α(x )c 2 b (x ) 2 t G(x, x s, t )dx dt
32 Born Approximation The Born approximation is made in the perturbation representation by substituting the total field under the integral for the background field. G(x, x s, t) = G b (x, x s, t) (27) G b (x, x, t t )α(x )c 2 b (x ) 2 t G b (x, x s, t )dx dt This is an explicit representation for G(x, x s, t). The perturbation can represent a (single additional) scattered wavefield as G s (x, x s, t) = d(x, x s, t) = G b (x, x, t t )α(x )c 2 b (x ) t 2 G b (x, x s, t )dx dt. (28)
The Finite Difference Method
The Finite Difference Method Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Heiner Igel Computational Seismology 1 / 32 Outline 1 Introduction Motivation
More informationWaveform inversion and time-reversal imaging in attenuative TI media
Waveform inversion and time-reversal imaging in attenuative TI media Tong Bai 1, Tieyuan Zhu 2, Ilya Tsvankin 1, Xinming Wu 3 1. Colorado School of Mines 2. Penn State University 3. University of Texas
More informationWave-equation tomography for anisotropic parameters
Wave-equation tomography for anisotropic parameters Yunyue (Elita) Li and Biondo Biondi ABSTRACT Anisotropic models are recognized as more realistic representations of the subsurface where complex geological
More informationIntroduction to Seismology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 12.510 Introduction to Seismology Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 12.510 Introduction
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 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 informationChapter 3. Finite Difference Methods for Hyperbolic Equations Introduction Linear convection 1-D wave equation
Chapter 3. Finite Difference Methods for Hyperbolic Equations 3.1. Introduction Most hyperbolic problems involve the transport of fluid properties. In the equations of motion, the term describing the transport
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 13
REVIEW Lecture 12: Spring 2015 Lecture 13 Grid-Refinement and Error estimation Estimation of the order of convergence and of the discretization error Richardson s extrapolation and Iterative improvements
More informationA 27-point scheme for a 3D frequency-domain scalar wave equation based on an average-derivative method
Geophysical Prospecting 04 6 58 77 doi: 0./365-478.090 A 7-point scheme for a 3D frequency-domain scalar wave equation based on an average-derivative method Jing-Bo Chen Key Laboratory of Petroleum Resources
More informationLowrank finite-differences and lowrank Fourier finite-differences for seismic wave extrapolation in the acoustic approximation a
Lowrank finite-differences and lowrank Fourier finite-differences for seismic wave extrapolation in the acoustic approximation a a Published in Geophysical Journal International, 193, 960-969, (2013) Xiaolei
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 informationOn the Earthquake-Source Numerical Implementation in the Seismic Wave Equation
Journal of Earthquake Engineering ISSN: 1363-2469 (Print) 1559-808X (Online) Journal homepage: http://www.tandfonline.com/loi/ueqe20 On the Earthquake-Source Numerical Implementation in the Seismic Wave
More informationSeismic wavepropagation concepts applied to the interpretation of marine controlled-source electromagnetics
Seismic wavepropagation concepts applied to the interpretation of marine controlled-source electromagnetics Rune Mittet, EMGS SUMMARY are very useful for the interpretation of seismic data. Moreover, these
More informationElastic wave-equation migration for laterally varying isotropic and HTI media. Richard A. Bale and Gary F. Margrave
Elastic wave-equation migration for laterally varying isotropic and HTI media Richard A. Bale and Gary F. Margrave a Outline Introduction Theory Elastic wavefield extrapolation Extension to laterally heterogeneous
More informationWave Phenomena Physics 15c. Lecture 11 Dispersion
Wave Phenomena Physics 15c Lecture 11 Dispersion What We Did Last Time Defined Fourier transform f (t) = F(ω)e iωt dω F(ω) = 1 2π f(t) and F(w) represent a function in time and frequency domains Analyzed
More informationTopic 5-1: Introduction to Phonons Kittel pages: 91, 92
Topic 5-1: Introduction to Phonons Kittel pages: 91, 92 Summary: In this video we introduce the concept that atoms are not rigid, fixed points within the lattice. Instead we treat them as quantum harmonic
More informationMIT (Spring 2014)
18.311 MIT (Spring 014) Rodolfo R. Rosales May 6, 014. Problem Set # 08. Due: Last day of lectures. IMPORTANT: Turn in the regular and the special problems stapled in two SEPARATE packages. Print your
More informationWebster s horn model on Bernoulli flow
Webster s horn model on Bernoulli flow Aalto University, Dept. Mathematics and Systems Analysis January 5th, 2018 Incompressible, steady Bernoulli principle Consider a straight tube Ω R 3 havin circular
More informationLecture 17: Initial value problems
Lecture 17: Initial value problems Let s start with initial value problems, and consider numerical solution to the simplest PDE we can think of u/ t + c u/ x = 0 (with u a scalar) for which the solution
More informationSound Waves Sound Waves:
3//18 Sound Waves Sound Waves: 1 3//18 Sound Waves Linear Waves compression rarefaction Inference of Sound Wave equation: Sound Waves We look at small disturbances in a compressible medium (note: compressible
More informationStolt residual migration for converted waves
Stanford Exploration Project, Report 11, September 18, 1, pages 1 6 Stolt residual migration for converted waves Daniel Rosales, Paul Sava, and Biondo Biondi 1 ABSTRACT P S velocity analysis is a new frontier
More informationMagnetohydrodynamic Waves
Magnetohydrodynamic Waves Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 17, 2016 These slides are largely based off of 4.5 and 4.8 of The Physics of
More informationOne-step extrapolation method for reverse time migration
GEOPHYSICS VOL. 74 NO. 4 JULY-AUGUST 2009; P. A29 A33 5 FIGS. 10.1190/1.3123476 One-step extrapolation method for reverse time migration Yu Zhang 1 and Guanquan Zhang 2 ABSTRACT We have proposed a new
More informationThe Schroedinger equation
The Schroedinger equation After Planck, Einstein, Bohr, de Broglie, and many others (but before Born), the time was ripe for a complete theory that could be applied to any problem involving nano-scale
More informationReview of Fundamental Equations Supplementary notes on Section 1.2 and 1.3
Review of Fundamental Equations Supplementary notes on Section. and.3 Introduction of the velocity potential: irrotational motion: ω = u = identity in the vector analysis: ϕ u = ϕ Basic conservation principles:
More informationLongitudinal Waves. waves in which the particle or oscillator motion is in the same direction as the wave propagation
Longitudinal Waves waves in which the particle or oscillator motion is in the same direction as the wave propagation Longitudinal waves propagate as sound waves in all phases of matter, plasmas, gases,
More informationStanford Exploration Project, Report 105, September 5, 2000, pages
Stanford Exploration Project, Report 5, September 5, 2, pages 57 66 56 Stanford Exploration Project, Report 5, September 5, 2, pages 57 66 Short Note Simple factorization of implicit finite-difference
More informationFundamentals of Fluid Dynamics: Waves in Fluids
Fundamentals of Fluid Dynamics: Waves in Fluids Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/ tzielins/ Institute
More informationTraveltime sensitivity kernels: Banana-doughnuts or just plain bananas? a
Traveltime sensitivity kernels: Banana-doughnuts or just plain bananas? a a Published in SEP report, 103, 61-68 (2000) James Rickett 1 INTRODUCTION Estimating an accurate velocity function is one of the
More informationTHE PHYSICS OF WAVES CHAPTER 1. Problem 1.1 Show that Ψ(x, t) = (x vt) 2. is a traveling wave.
CHAPTER 1 THE PHYSICS OF WAVES Problem 1.1 Show that Ψ(x, t) = (x vt) is a traveling wave. Show thatψ(x, t) is a wave by substitutioninto Equation 1.1. Proceed as in Example 1.1. On line version uses Ψ(x,
More informationFourier finite-difference wave propagation a
Fourier finite-difference wave propagation a a Published in Geophysics, 76, T123-T129, (2011) Xiaolei Song and Sergey Fomel ABSTRACT We introduce a novel technique for seismic wave extrapolation in time.
More informationA frequency-space 2-D scalar wave extrapolator using extended 25-point finite-difference operator
GEOPHYSICS, VOL. 63, NO. 1 (JANUARY-FEBRUARY 1998); P. 289 296, 7 FIGS. A frequency-space 2-D scalar wave extrapolator using extended 25-point finite-difference operator Changsoo Shin and Heejeung Sohn
More informationExplicit extrapolators and common azimuth migration
Explicit extrapolators and common azimuth migration Eric Dussaud October 9th, 004 Agenda Downward continuation via explicit extrapolation. Design of explicit extrapolators. Design of 1-D convolution operators.
More informationA second look at waves
A second loo at waves ravelling waves A first loo at Amplitude Modulation (AM) Stationary and reflected waves Lossy waves: dispersion & evanescence I thin this is the MOS IMPORAN of my eight lectures,
More informationAsymptotic Behavior of Waves in a Nonuniform Medium
Available at http://pvamuedu/aam Appl Appl Math ISSN: 1932-9466 Vol 12, Issue 1 June 217, pp 217 229 Applications Applied Mathematics: An International Journal AAM Asymptotic Behavior of Waves in a Nonuniform
More information1 Introduction to MATLAB
L3 - December 015 Solving PDEs numerically (Reports due Thursday Dec 3rd, carolinemuller13@gmail.com) In this project, we will see various methods for solving Partial Differential Equations (PDEs) using
More informationPositive Definite Kernels: Opportunities and Challenges
Positive Definite Kernels: Opportunities and Challenges Michael McCourt Department of Mathematical and Statistical Sciences University of Colorado, Denver CUNY Mathematics Seminar CUNY Graduate College
More informationPseudo-acoustic wavefield propagation for anisotropic media
Anisotropic wavefield propagation Pseudo-acoustic wavefield propagation for anisotropic media Ben D. Wards, Gary F. Margrave, and Michael P. Lamoureux ABSTRACT Reverse-time migration (RTM) is a powerful
More informationLinear Advection Equation
Linear Advection Equation Numerical Methods Laboratory Modelling the Climate System ARC Centre of Excellence for Climate System Science 2nd Annual Winter School 18 June 2013 1 Analytical Solution The linear
More information3D acoustic wave modeling with a time-space domain dispersion-relation-based Finite-difference scheme
P-8 3D acoustic wave odeling with a tie-space doain dispersion-relation-based Finite-difference schee Yang Liu * and rinal K. Sen State Key Laboratory of Petroleu Resource and Prospecting (China University
More informationPEAT SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity
PEAT8002 - SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity Nick Rawlinson Research School of Earth Sciences Australian National University Anisotropy Introduction Most of the theoretical
More informationq t = F q x. (1) is a flux of q due to diffusion. Although very complex parameterizations for F q
! Revised Tuesday, December 8, 015! 1 Chapter 7: Diffusion Copyright 015, David A. Randall 7.1! Introduction Diffusion is a macroscopic statistical description of microscopic advection. Here microscopic
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 informationS=z y. p(x, t) v. = ρ x (x,t)
LTUR 2: One-imensional Traveling Waves Main Points - xponential and sine-wave solutions to the one-dimensional wave equation. - The distributed compressibility and mass in acoustic plane waves are analogous
More informationWAVES CP4 REVISION LECTURE ON. The wave equation. Traveling waves. Standing waves. Dispersion. Phase and group velocities.
CP4 REVISION LECTURE ON WAVES The wave equation. Traveling waves. Standing waves. Dispersion. Phase and group velocities. Boundary effects. Reflection and transmission of waves. !"#$%&''(%)*%+,-.%/%+,01%
More informationAmplitude preserving AMO from true amplitude DMO and inverse DMO
Stanford Exploration Project, Report 84, May 9, 001, pages 1 167 Amplitude preserving AMO from true amplitude DMO and inverse DMO Nizar Chemingui and Biondo Biondi 1 ABSTRACT Starting from the definition
More informationFundamentals of Acoustics
Fundamentals of Acoustics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy of Sciences
More informationSolution of Additional Exercises for Chapter 4
1 1. (1) Try V (x) = 1 (x 1 + x ). Solution of Additional Exercises for Chapter 4 V (x) = x 1 ( x 1 + x ) x = x 1 x + x 1 x In the neighborhood of the origin, the term (x 1 + x ) dominates. Hence, the
More informationLINEAR DISPERSIVE WAVES
LINEAR DISPERSIVE WAVES Nathaniel Egwu December 16, 2009 Outline 1 INTRODUCTION Dispersion Relations Definition of Dispersive Waves 2 SOLUTION AND ASYMPTOTIC ANALYSIS The Beam Equation General Solution
More informationPeriodic functions: simple harmonic oscillator
Periodic functions: simple harmonic oscillator Recall the simple harmonic oscillator (e.g. mass-spring system) d 2 y dt 2 + ω2 0y = 0 Solution can be written in various ways: y(t) = Ae iω 0t y(t) = A cos
More informationSound Propagation through Media. Nachiketa Tiwari Indian Institute of Technology Kanpur
Sound Propagation through Media Nachiketa Tiwari Indian Institute of Technology Kanpur LECTURE-17 HARMONIC PLANE WAVES Introduction In this lecture, we discuss propagation of 1-D planar waves as they travel
More informationThe Pseudospectral Method
The Pseudospectral Method Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Heiner Igel Computational Seismology 1 / 34 Outline 1 Infinite Order Finite Differences
More informationStability and dispersion analysis of high order FDTD methods for Maxwell s equations in dispersive media
Contemporary Mathematics Volume 586 013 http://dx.doi.org/.90/conm/586/11666 Stability and dispersion analysis of high order FDTD methods for Maxwell s equations in dispersive media V. A. Bokil and N.
More informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
More informationTodd Arbogast. Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and Sciences (ICES)
CONSERVATIVE CHARACTERISTIC METHODS FOR LINEAR TRANSPORT PROBLEMS Todd Arbogast Department of Mathematics and, (ICES) The University of Texas at Austin Chieh-Sen (Jason) Huang Department of Applied Mathematics
More informationIntegration of Rational Functions by Partial Fractions
Title Integration of Rational Functions by MATH 1700 MATH 1700 1 / 11 Readings Readings Readings: Section 7.4 MATH 1700 2 / 11 Rational functions A rational function is one of the form where P and Q are
More informationONE-WAY AND ONE-RETURN APPROXIMATIONS (DE WOLF APPROXIMATION) FOR FAST ELASTIC WAVE MODELING IN COMPLEX MEDIA
ADVANCES IN GEOPHYSICS, VOL. 48, CHAPTER 5 ONE-WAY AND ONE-RETURN APPROXIMATIONS (DE WOLF APPROXIMATION) FOR FAST ELASTIC WAVE MODELING IN COMPLEX MEDIA RU-SHAN WU, XIAO-BI XIE AND XIAN-YUN WU Modeling
More informationModel Order Reduction and Stability Enforcement of Finite-Difference Time-Domain Equations Beyond the CFL Limit. Xihao Li
Model Order Reduction and Stability Enforcement of Finite-Difference Time-Domain Equations Beyond the CFL Limit by Xihao Li A thesis submitted in conformity with the requirements for the degree of Master
More informationWave Propagation in Heterogeneous Media: Born and Rytov Approximations. Chris Sherman
Wave Propagation in Heterogeneous Media: Born and Rytov Approximations Chris Sherman Stochastic Scalar Wave Equation Wave Equation: &! " %& 1 # t V x ' ) () u(x,t) = 0 (*) Velocity Perturbation: = V o
More informationwhere d is the vibration direction of the displacement and c is the wave velocity. For a fixed time t,
3 Plane waves 3.1 Plane waves in unbounded solid Consider a plane wave propagating in the direction with the unit vector p. The displacement of the plane wave is assumed to have the form ( u i (x, t) =
More informationPlane Wave: Introduction
Plane Wave: Introduction According to Mawell s equations a timevarying electric field produces a time-varying magnetic field and conversely a time-varying magnetic field produces an electric field ( i.e.
More informationBASIC WAVE CONCEPTS. Reading: Main 9.0, 9.1, 9.3 GEM 9.1.1, Giancoli?
1 BASIC WAVE CONCEPTS Reading: Main 9.0, 9.1, 9.3 GEM 9.1.1, 9.1.2 Giancoli? REVIEW SINGLE OSCILLATOR: The oscillation functions you re used to describe how one quantity (position, charge, electric field,
More informationWaves, the Wave Equation, and Phase Velocity. We ll start with optics. The one-dimensional wave equation. What is a wave? Optional optics texts: f(x)
We ll start with optics Optional optics texts: Waves, the Wave Equation, and Phase Velocity What is a wave? f(x) f(x-) f(x-) f(x-3) Eugene Hecht, Optics, 4th ed. J.F. James, A Student's Guide to Fourier
More informationEfficient gradient computation of a misfit function for FWI using the adjoint method
Efficient gradient computation of a misfit function for FWI using the adjoint method B. Galuzzi, department of Earth Sciences A. Desio, University of Milan E. Stucchi, department of Earth Sciences A. Desio,
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 informationSeismic Imaging. William W. Symes. C. I. M. E. Martina Franca September
Seismic Imaging William W. Symes C. I. M. E. Martina Franca September 2002 www.trip.caam.rice.edu 1 0 offset (km) -4-3 -2-1 1 2 time (s) 3 4 5 How do you turn lots of this... (field seismogram from the
More informationSUMMARY. fractured reservoirs in the context of exploration seismology.
Modeling scattering of cross-well seismic waves using Radiative Transfer Theory Josimar A. Da Silva Jr, Oleg V. Poliannikov and Michael Fehler, Earth Resources Laboratory / M.I.T SUMMARY We model P and
More informationProf. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit
Prof. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit Central concepts: Phase velocity: velocity with which surfaces of constant phase move Group velocity: velocity with which slow
More informationCONSTRUCTION AND ANALYSIS OF WEIGHTED SEQUENTIAL SPLITTING FDTD METHODS FOR THE 3D MAXWELL S EQUATIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 5 Number 6 Pages 747 784 c 08 Institute for Scientific Computing and Information CONSTRUCTION AND ANALYSIS OF WEIGHTED SEQUENTIAL SPLITTING
More informationIntegration of Rational Functions by Partial Fractions
Title Integration of Rational Functions by Partial Fractions MATH 1700 December 6, 2016 MATH 1700 Partial Fractions December 6, 2016 1 / 11 Readings Readings Readings: Section 7.4 MATH 1700 Partial Fractions
More informationIntroduction to Waves in Structures. Mike Brennan UNESP, Ilha Solteira São Paulo Brazil
Introduction to Waves in Structures Mike Brennan UNESP, Ilha Solteira São Paulo Brazil Waves in Structures Characteristics of wave motion Structural waves String Rod Beam Phase speed, group velocity Low
More informationFundamental Solution
Fundamental Solution onsider the following generic equation: Lu(X) = f(x). (1) Here X = (r, t) is the space-time coordinate (if either space or time coordinate is absent, then X t, or X r, respectively);
More informationExam in TMA4195 Mathematical Modeling Solutions
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Exam in TMA495 Mathematical Modeling 6..07 Solutions Problem a Here x, y are two populations varying with time
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 informationCourse on Inverse Problems
California Institute of Technology Division of Geological and Planetary Sciences March 26 - May 25, 2007 Course on Inverse Problems Albert Tarantola Institut de Physique du Globe de Paris Lesson XVI CONCLUSION
More informationSeismic inverse scattering by reverse time migration
Seismic inverse scattering by reverse time migration Chris Stolk 1, Tim Op t Root 2, Maarten de Hoop 3 1 University of Amsterdam 2 University of Twente 3 Purdue University MSRI Inverse Problems Seminar,
More informationCourse on Inverse Problems Albert Tarantola
California Institute of Technology Division of Geological and Planetary Sciences March 26 - May 25, 27 Course on Inverse Problems Albert Tarantola Institut de Physique du Globe de Paris CONCLUSION OF THE
More informationLinear second-order differential equations with constant coefficients and nonzero right-hand side
Linear second-order differential equations with constant coefficients and nonzero right-hand side We return to the damped, driven simple harmonic oscillator d 2 y dy + 2b dt2 dt + ω2 0y = F sin ωt We note
More informationComputational Seismology: An Introduction
Computational Seismology: An Introduction Aim of lecture: Understand why we need numerical methods to understand our world Learn about various numerical methods (finite differences, pseudospectal methods,
More informationWe briefly discuss two examples for solving wave propagation type problems with finite differences, the acoustic and the seismic problem.
Excerpt from GEOL557 Numerical Modeling of Earth Systems by Becker and Kaus 2016 1 Wave propagation Figure 1: Finite difference discretization of the 2D acoustic problem. We briefly discuss two examples
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 information3.3. Phase and Amplitude Errors of 1-D Advection Equation
3.3. Phase and Amplitude Errors of 1-D Advection Equation Reading: Duran section 2.4.2. Tannehill et al section 4.1.2. The following example F.D. solutions of a 1D advection equation show errors in both
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 informationNumerical Hydraulics
ETH Zurich, Fall 2017 Numerical Hydraulics Assignment 2 Numerical solution of shallow water wave propagation (www.surfertoday.com) 1 Introduction 1.1 Equations Understanding the propagation of shallow
More informationDynamic Meteorology (lecture 13, 2016)
Dynamic Meteorology (lecture 13, 2016) Topics Chapter 3, lecture notes: High frequency waves (no rota;on) Boussinesq approxima4on Normal mode analysis of the stability of hydrosta4c balance Buoyancy waves
More information6.1. Linearized Wave Equations in a Uniform Isotropic MHD Plasma. = 0 into Ohm s law yields E 0
Chapter 6. Linear Waves in the MHD Plasma 85 Chapter 6. Linear Waves in the MHD Plasma Topics or concepts to learn in Chapter 6:. Linearize the MHD equations. The eigen-mode solutions of the MHD waves
More informationThermoacoustic Devices
Thermoacoustic Devices Peter in t panhuis Sjoerd Rienstra Han Slot 9th May 2007 Introduction Thermoacoustics All effects in acoustics in which heat conduction and entropy variations play a role. (Rott,
More informationNon-Continuum Energy Transfer: Phonons
Non-Continuum Energy Transfer: Phonons D. B. Go Slide 1 The Crystal Lattice The crystal lattice is the organization of atoms and/or molecules in a solid simple cubic body-centered cubic hexagonal a NaCl
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationCourant and all that. Consistency, Convergence Stability Numerical Dispersion Computational grids and numerical anisotropy
Consistency, Convergence Stability Numerical Dispersion Computational grids and numerical anisotropy c dx ε 1 he goal of this lecture is to understand how to find suitable parameters (i.e., grid spacing
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 informationCalculation of Reflection and Transmission Coefficients in scuff-transmission
Calculation of Reflection and Transmission Coefficients in scuff-transmission Homer Reid May 9, 2015 Contents 1 The Setup 2 2 Scattering coefficients from surface currents 4 2.1 Computation of b(q).........................
More informationBasics on Numerical Methods for Hyperbolic Equations
Basics on Numerical Methods for Hyperbolic Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8,
More informationTitle: Exact wavefield extrapolation for elastic reverse-time migration
Title: Exact wavefield extrapolation for elastic reverse-time migration ummary: A fundamental step of any wave equation migration algorithm is represented by the numerical projection of the recorded data
More information(TRAVELLING) 1D WAVES. 1. Transversal & Longitudinal Waves
(TRAVELLING) 1D WAVES 1. Transversal & Longitudinal Waves Objectives After studying this chapter you should be able to: Derive 1D wave equation for transversal and longitudinal Relate propagation speed
More informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
More informationPhysics 505 Homework No. 12 Solutions S12-1
Physics 55 Homework No. 1 s S1-1 1. 1D ionization. This problem is from the January, 7, prelims. Consider a nonrelativistic mass m particle with coordinate x in one dimension that is subject to an attractive
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationElastodynamic single-sided homogeneous Green's function representation: Theory and examples
Delft University of Technology Elastodynamic single-sided homogeneous Green's function representation: Theory and examples Wapenaar, Kees; Reinicke Urruticoechea, Christian DOI 10.3997/2214-4609.201701220
More informationThe Evolution of Large-Amplitude Internal Gravity Wavepackets
The Evolution of Large-Amplitude Internal Gravity Wavepackets Sutherland, Bruce R. and Brown, Geoffrey L. University of Alberta Environmental and Industrial Fluid Dynamics Laboratory Edmonton, Alberta,
More information