The Elastic Wave Equation. The elastic wave equation

Similar documents
Scattering at an Interface: Oblique Incidence

by Lauren DeDieu Advisor: George Chen

CONSISTENT EARTHQUAKE ACCELERATION AND DISPLACEMENT RECORDS

II. Light is a Ray (Geometrical Optics)

Dynamic Model of the Axially Moving Viscoelastic Belt System with Tensioner Pulley Yanqi Liu1, a, Hongyu Wang2, b, Dongxing Cao3, c, Xiaoling Gai1, d

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4

CH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

Calculation of the Resistance of a Ship Mathematical Formulation. Calculation of the Resistance of a Ship Mathematical Formulation

Different kind of oscillation

Lecture 9: Dynamic Properties

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems

NATIONAL UNIVERSITY OF SINGAPORE PC5202 ADVANCED STATISTICAL MECHANICS. (Semester II: AY ) Time Allowed: 2 Hours

Mechanics Physics 151

TRANSIENTS. Lecture 5 ELEC-E8409 High Voltage Engineering

Cartesian tensors. Order (rank) Scalar. Vector. 3x3 matrix

Chapters 2 Kinematics. Position, Distance, Displacement

, t 1. Transitions - this one was easy, but in general the hardest part is choosing the which variables are state and control variables

2.1 Constitutive Theory

Stochastic Programming handling CVAR in objective and constraint

OPTIMIZATION OF A NONCONVENTIONAL ENGINE EVAPORATOR

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD

Anisotropy and oblique total transmission at a planar negative-index interface

Linear Response Theory: The connection between QFT and experiments

Let s treat the problem of the response of a system to an applied external force. Again,

Mechanics Physics 151

Mechanics Physics 151

( ) () we define the interaction representation by the unitary transformation () = ()

Sound Transmission Throough Lined, Composite Panel Structures: Transversely Isotropic Poro- Elastic Model

DEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL

Solution of a diffusion problem in a non-homogeneous flow and diffusion field by the integral representation method (IRM)

CS286.2 Lecture 14: Quantum de Finetti Theorems II

Normal Random Variable and its discriminant functions

Bandlimited channel. Intersymbol interference (ISI) This non-ideal communication channel is also called dispersive channel

. The geometric multiplicity is dim[ker( λi. number of linearly independent eigenvectors associated with this eigenvalue.

. The geometric multiplicity is dim[ker( λi. A )], i.e. the number of linearly independent eigenvectors associated with this eigenvalue.

Density Matrix Description of NMR BCMB/CHEM 8190

( ) [ ] MAP Decision Rule

Variants of Pegasos. December 11, 2009

Lecture 5. Plane Wave Reflection and Transmission

Solution in semi infinite diffusion couples (error function analysis)

PHYS 1443 Section 001 Lecture #4

Motion in Two Dimensions

Variational method to the second-order impulsive partial differential equations with inconstant coefficients (I)

Research Article Cubic B-spline for the Numerical Solution of Parabolic Integro-differential Equation with a Weakly Singular Kernel

Chapter 5 Mobile Radio Propagation: Small-Scale Scale Fading and Multipath

Study on the Unexpected Wave of the Experiment of Ultrasonic Guided Wave in Square Steel Bar based on 2D Equivalent Simulation

Online Appendix for. Strategic safety stocks in supply chains with evolving forecasts

Observer Design for Nonlinear Systems using Linear Approximations

Velocity Modeling in a Vertical Transversely Isotropic Medium Using Zelt Method

Motion of Wavepackets in Non-Hermitian. Quantum Mechanics

TSS = SST + SSE An orthogonal partition of the total SS

10. A.C CIRCUITS. Theoretically current grows to maximum value after infinite time. But practically it grows to maximum after 5τ. Decay of current :

Density Matrix Description of NMR BCMB/CHEM 8190

J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes.

Outline. Review Solution Approaches. Review Basic Equations. Nature of Turbulence. Review Fluent Exercise. Turbulence Models

Electromagnetic waves in vacuum.

Turbulence Modelling (CFD course)

Bayes rule for a classification problem INF Discriminant functions for the normal density. Euclidean distance. Mahalanobis distance

XIII International PhD Workshop OWD 2011, October Three Phase DC/DC Boost Converter With High Energy Efficiency

Introduction to elastic wave equation. Salam Alnabulsi University of Calgary Department of Mathematics and Statistics October 15,2012

MEEN Handout 4a ELEMENTS OF ANALYTICAL MECHANICS

First-order piecewise-linear dynamic circuits

Chapter Finite Difference Method for Ordinary Differential Equations

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS

Robustness Experiments with Two Variance Components

Numerical simulation of flow reattachment length in a stilling basin with a step-down floor

Comb Filters. Comb Filters

M. Y. Adamu Mathematical Sciences Programme, AbubakarTafawaBalewa University, Bauchi, Nigeria

ELASTIC MODULUS ESTIMATION OF CHOPPED CARBON FIBER TAPE REINFORCED THERMOPLASTICS USING THE MONTE CARLO SIMULATION

Anisotropic Behaviors and Its Application on Sheet Metal Stamping Processes

Chapter 6 Conservation Equations for Multiphase- Multicomponent Flow Through Porous Media

Lecture 6: Learning for Control (Generalised Linear Regression)

Ordinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s

FI 3103 Quantum Physics

Chapter 6 DETECTION AND ESTIMATION: Model of digital communication system. Fundamental issues in digital communications are

Chapter 6 DETECTION AND ESTIMATION: Model of digital communication system. Fundamental issues in digital communications are

Approximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy

THERMODYNAMICS 1. The First Law and Other Basic Concepts (part 2)

Response of MDOF systems

Lect. 12: Oblique Incidence at Dielectric Interface

Notes on the stability of dynamic systems and the use of Eigen Values.

Displacement, Velocity, and Acceleration. (WHERE and WHEN?)

FTCS Solution to the Heat Equation

Diffraction by dielectric wedges: high frequency and time domain solutions

Real Time Hybrid Simulation using Shaking Tabels

Real-Time Trajectory Generation and Tracking for Cooperative Control Systems

Chapter Lagrangian Interpolation

Transcription: Messenger RNA, mrna, is produced and transported to Ribosomes

Incident ray Reflected ray θ i. θ r

On One Analytic Method of. Constructing Program Controls

3.2 Models for technical systems

Born Oppenheimer Approximation and Beyond

Lecture VI Regression

Dynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005

CHAPTER 10: LINEAR DISCRIMINATION

Lect. 13: Oblique Incidence at Dielectric Interface

WebAssign HW Due 11:59PM Tuesday Clicker Information

A symmetric boundary integral approach to transient poroelastic analysis

Lesson 2 Transmission Lines Fundamentals

Wronskian Determinant Solutions for the (3 + 1)-Dimensional Boiti-Leon-Manna-Pempinelli Equation

Transcription:

The Elasc Wave Eqaon Elasc waves n nfne homogeneos soropc meda Nmercal smlaons for smple sorces Plane wave propagaon n nfne meda Freqency, wavenmber, wavelengh Condons a maeral dsconnes nell s Law Reflecon coeffcens Free srface Reflecon esmology: an example from he Glf of Mexco

Eqaons of moon ρ f + j σ j Wha are he solons o hs eqaon? A frs we look a nfne homogeneos soropc meda, hen: σ j λθδ j + µε j σ j λ δ + µ k k j ( j + j ) ρ f + j k k j µ ( λ δ + ( + )) ρ f + λ + µ + µ k k j j j j j

Eqaons of moon homogeneos meda ρ f + λ + µ + µ k k j j j We can now smplfy hs eqaon sng he crl and dv operaors and x + - ρ f+ ( λ + µ ) - µ y + z hs holds n any coordnae sysem Ths eqaon can be frher smplfed, separang he wavefeld no crl free and dv free pars

Eqaons of moon P waves ρ ( λ + µ ) - µ Le s apply he dv operaor o hs eqaon, we oban ρ θ ( λ + µ ) θ Acosc wave eqaon where θ P wave velocy or θ α 1 θ α λ + µ ρ

Eqaons of moon shear waves ρ ( λ + µ ) - µ Le s apply he crl operaor o hs eqaon, we oban ρ ( λ + µ ) θ + µ ( ) we now make se of Wave eqaon for shear waves θ ϕ and defne o oban hear wave velocy ϕ β 1 ϕ β µ ρ

Elasodynamc Poenals Any vecor may be separaed no scalar and vecor poenals Φ + Ψ where P s he poenal for Φ waves and Ψ he poenal for shear waves θ Φ ϕ Ψ Ψ P-waves have no roaon hear waves have no change n volme θ α θ ϕ β ϕ

esmc Veloces Maeral and orce P-wave velocy (m/s) shear wave velocy (m/s) Waer 15 Loose sand 18 5 Clay 11-5 andsone 14-43 Anhydre, Glf Coas 41 Conglomerae 4 Lmesone 63 33 Grane 564 87 Granodore 478 31 Dore 578 36 Basal 64 3 Dne 8 437 Gabbro 645 34

olons o he wave eqaon - general Le s consder a regon who sorces η c η Where n cold be eher dlaaon or he vecor poenal and c s eher P- or shear-wave velocy. The general solon o hs eqaon s: η( x, ) G( a x c) j j ± Le s ake a look a a 1-D example

olons o he wave eqaon - harmonc Le s consder a regon who sorces η c η The mos approprae choce for G s of corse he se of harmonc fncons: ( x, ) A exp[ k( a x c)] j j

olons o he wave eqaon - harmonc akng only he real par and consderng only 1D we oban ( x, ) Acos[ k( x c)] k x c kx kc π π π ( ) x ω x λ λ T c k λ T ω A wave speed wavenmber wavelengh perod freqency amplde

phercal Waves η c η Le s assme ha η s a fncon of he dsance from he sorce r η 1 rη + rη r c η where we sed he defnon of he Laplace operaor n sphercal coordnaes le s defne o oban η η r η c η wh he known solon η f ( r α)

Geomercal spreadng so a dsrbance propagang away wh sphercal wavefrons decays lke r 1 η f ( r α) η r 1 r... hs s he geomercal spreadng for sphercal waves, he amplde decays proporonal o 1/r. If we had looked a cylndrcal waves he resl wold have been ha he waves decay as (e.g. srface waves) η 1 r

Plane waves... wha can we say abo he drecon of dsplacemen, he polarzaon of sesmc waves? Φ + Ψ P Φ P + Ψ Φ... we now assme ha he poenals have he well known form of plane harmonc waves Aexp( k x ω) Ψ B exp( k x ω) P Akexp( k x ω) k B exp( k x ω) P waves are longdnal as P s parallel o k shear waves are ransverse becase s normal o he wave vecor k

Heerogenees.. Wha happens f we have heerogenees? Dependng on he knd of reflecon par or all of he sgnal s refleced or ransmed. Wha happens a a free srface? Can a P wave be convered n an wave or vce versa? How bg are he ampldes of he refleced waves?

Bondary Condons... wha happens when he maeral parameers change? ρ 1 v 1 welded nerface ρ v A a maeral nerface we reqre conny of dsplacemen and racon A specal case s he free srface condon, where he srface racons are zero.

Reflecon and Transmsson nell s Law Wha happens a a (fla) maeral dsconny? Medm 1: v 1 1 sn sn 1 v v 1 Medm : v B how mch s refleced, how mch ransmed?

Reflecon and Transmsson coeffcens Le s ake he mos smple example: P-waves wh normal ncdence on a maeral nerface Medm 1: r1,v1 A R R A ρα ρ α ρ1α 1 + ρ α 1 1 Medm : r,v T T A ρ1α 1 ρ α + ρ α 1 1 A oblqe angles conversons from -P, P- have o be consdered.

Reflecon and Transmsson Ansaz How can we calclae he amon of energy ha s ransmed or refleced a a maeral dsconny? We know ha n homogeneos meda he dsplacemen can be descrbed by he correspondng poenals esmology and he Earh s Deep Ineror n -D hs yelds Φ + Ψ x y z x z z Φ Ψ x Φ + y x Ψ x Ψ an ncomng P wave has he form Φ z Ψ ω A exp ( a j x j α) α y z

Reflecon and Transmsson Ansaz... here a are he componens of he vecor normal o he wavefron : a (cos e,, -sn e), where e s he angle beween srface and ray drecon, so ha for he free srface Φ Ψ A expk( x B expk'( x 1 1 + x x 3 3 an e c) + an f c' ) Aexpk( x 1 + x 3 an e c) where c k α cos e ω cos e α ω c c' k' β cos f ω cos β f P e f P r V r wha we know s ha σ σ xz zz

Reflecon and Transmsson Coeffcens... png he eqaons for he poenals (dsplacemens) no hese eqaons leads o a relaon beween ncden and refleced (ransmed) ampldes R PP A A 4 an e an 4 an e an f f (1 + (1 an an f f ) ) R P A A 4 an e (1 an f ) 4 an e an f + (1 an f ) These are he reflecon coeffcens for a plane P wave ncden on a free srface, and refleced P and V waves.

Case 1: Reflecons a a free srface A P wave s ncden a he free srface... j P V P The refleced ampldes can be descrbed by he scaerng marx P P P d d P d d

Case : H waves For layered meda H waves are compleely decopled from P and V waves H There s no converson only H waves are refleced or ransmed d d d d

Case 3: old-sold nerface P V r P r V P To accon for all possble reflecons and ransmssons we need 16 coeffcens, descrbed by a 4x4 scaerng marx.

Case 4: old-fld nerface P V r P r P A a sold-fld nerface here s no converson o V n he lower medm.

Reflecon coeffcens - example

Reflecon coeffcens - example

Refracons waveform effecs

caerng and Aenaon Propagang sesmc waves loose energy de o geomercal spreadng e.g. he energy of sphercal wavefron emanang from a pon sorce s dsrbed over a sphercal srface of ever ncreasng sze nrnsc aenaon elasc wave propagaon consss of a permanen exchange beween poenal (dsplacemen) and knec (velocy) energy. Ths process s no compleely reversble. There s energy loss de o shear heang a gran bondares, mneral dslocaons ec. scaerng aenaon whenever here are maeral changes he energy of a wavefeld s scaered n dfferen phases. Dependng on he maeral properes hs wll lead o amplde decay and dspersve effecs.

Inrnsc aenaon How can we descrbe nrnsc aenaon? Le s ry a sprng model: The eqaon of moon for a damped harmonc oscllaor s mx + γx γ x + x + m x + εϖ + k m kx x x + ω x ϖ ε γ mϖ k m 1/ where ε s he frcon coeffcen.

Q The solon o hs sysem s x( ) A e εϖ sn( ϖ 1 ε ) so we have a me-dependen amplde of A ( ) A e εϖ A e ϖ Q and defnng 1 A1 ε δ ln Q A Q π δ Q s he energy loss per cycle. Inrnsc aenaon n he Earh s n general descrbed by Q.

Dsperson effecs Wha happens f we have freqency ndependen Q,.e. each freqency looses he same amon of energy per cycle? A( x) A e ( fπqv) x hgh freqences more oscllaons more aenaon low freqences less oscllaons less aenaon Conseqences: - hgh freqences decay very rapdly - plse broadenng In he Earh we observe ha Q p s large han Q. Ths s de o he fac ha nrnsc aenaon s predomnanly cased by shear lace effecs a gran bondares.

Q n he Earh Rock Type Q p Q hale 3 1 andsone 58 31 Grane 5 7-5 Perdoe Mdmanle Lowermanle Oer Core 65 36 1 8 8 5

caerng