Frequency-domain ray series for viscoelastic waves with a non-symmetric stiffness matrix
|
|
- Darren Hardy
- 5 years ago
- Views:
Transcription
1 Frequency-domain ray series for viscoelastic waves with a non-symmetric stiffness matrix Ludě Klimeš Department of Geophysics, Faculty of Mathematics Physics, Charles University, Ke Karlovu 3, Praha 2, Czech Republic ( Received: September 8, 2017; Revised: November 29, 2017; Accepted: December 14, 2017 ABSTRACT In an elastic medium, it was proved that the stiffness tensor is symmetric with respect to the exchange of the first pair of indices the second pair of indices, but the proof does not apply to a viscoelastic medium. In order to indicate which phenomena could be observed in the wave field if the stiffness matrix were non symmetric, we propose the frequency domain ray series for viscoelastic waves with a non symmetric stiffness tensor in this paper. K e y wo r d s: viscoelastic media, stiffness tensor, wave propagation, ray theory, ray series, travel time, amplitude 1. INTRODUCTION The complex valued frequency domain stiffness tensor (elastic tensor, tensor of elastic moduli) c ijl = c ijl (x m, ω), projecting the strain tensor onto the stress tensor, depends on spatial coordinates x m circular frequency ω. It is symmetric with respect to the first pair of indices with respect to the second pair of indices c ijl = c jil (1) c ijl = c ijl. (2) It is thus frequently expressed in the form of the 6 6 stiffness matrix whose lines correspond to the first pair of indices columns to the second pair of indices according to the Voigt notation. In an elastic medium, it was proved that the stiffness tensor is symmetric with respect to the exchange of the first pair of indices the second pair of indices, c ijl = c lij. (3) The 6 6 stiffness matrix is thus symmetric in an elastic medium. The ray theory for a 3 D smoothly heterogeneous anisotropic elastic medium with a symmetric Stud. Geophys. Geod., 62 (2018), , DOI: /s c 2018 The Authors 261
2 L. Klimeš stiffness matrix is well developed (Babich, 1961; Červený, 1972; 2001; Chapman, 2004; Červený et al., 2007). However, the above mentioned proof does not apply to a viscoelastic medium, we do not now whether symmetry (3) holds in a viscoelastic medium. It is then reasonable to indicate which phenomena could be observed in the wave field if the stiffness matrix were non symmetric. Analogously to Thomson (1997), we thus do not assume symmetry (3) in this paper, propose the frequency domain ray series for viscoelastic waves with a non symmetric stiffness matrix, c ijl c lij. (4) The lower case Roman indices tae values 1, 2 3. The Einstein summation over repetitive lower case Roman indices is used throughout the paper. 2. STANDARD RAY SERIES In the frequency domain, the viscoelastodynamic equation for complex valued displacement u i = u i (x m, ω) reads (c ijl u l, ),j (iω) 2 u i = 0, (5) where lower case Roman subscript, following a comma denotes the partial derivative with respect to the corresponding spatial coordinate x. Here c ijl = c ijl (x m, ω) is the complex valued frequency domain stiffness tensor, = (x m ) is the density ω is the circular frequency. We express the displacement in terms of its complex valued vectorial amplitude U i = U i (x m, ω) complex valued travel time τ = τ(x m ) as u i = U i exp(iωτ). (6) We exp the frequency dependent amplitude into high frequency asymptotic series U i = (iω) n U [n] i, (7) where U [n] i n=0 = U [n] i (x m, ω) is the n th order vectorial amplitude. We insert displacement (6) into viscoelastodynamic equation (5) obtain equation (iω) 2 N i (U m, τ,u ) + iωm i (U m, τ,u ) + L i (U m ) = 0, (8) where the differential operators are defined as Here the Christoffel matrix, defined as N i (U m, τ,u ) = [Γ il (x m, τ,u )U l U i ], (9) M i (U m, τ,u ) = (c ijl τ, U l ),j + c ijl τ,j U l,, (10) L i (U m ) = (c ijl U l, ),j. (11) Γ il (x m, p u ) = a ijl (x m ) p j p, (12) 262 Stud. Geophys. Geod., 62 (2018)
3 Frequency-domain ray series for viscoelastic waves with a non-symmetric stiffness matrix is a function of six phase space coordinates x m, p u formed by three spatial coordinates x m three slowness vector components p u which represent the spatial gradient of travel time τ, p u = τ,u. In definition (12), a ijl (x m ) = c ijl (x m ) [ (x u ) ] 1 is the density reduced stiffness tensor. In order to obtain the stard ray series, we insert series (7) into viscoelastodynamic equation (8) sort the terms according to the order of iω, analogously to Babich (1961) or Červený (1972; 2001, Sec. 5.7). We then obtain the system of equations N i( U [n], τ,l ) + M i ( U [n 1], τ,l ) + L i ( U [n 2] (13) ) = 0 (14) for each order n = 0, 1, 2,... Here U [ 1] = 0 U [ 2] = 0, i.e. operator M i is missing in this equation for n = 0 operator L i is missing in this equation for n = 0, 1. Note that the system of equations (14) does not follow necessarily from viscoelastodynamic equation (8). The above sorting according to the order of iω is our assumption leading to the stard ray series. If we wish to obtain another ind of ray series, e.g., the coupling ray series, we may choose a solution of viscoelastodynamic equation (8) different from the system of equations (14), see Klimeš (2013). 3. CHRISTOFFEL EQUATION AND EIKONAL EQUATION Equation (14) for n = 0 constitutes the matrix Christoffel equation Γ il (x m, τ,u )U [0] l U [0] i = 0. (15) The three eigenvalues of Christoffel matrix (12) correspond to three waves: the P wave two S waves. Unlie as in the elastic case, the Christoffel matrix is not symmetric. We select one of three eigenvalues of Christoffel matrix (12) denote it as G = G(x m, τ,u ). We denote the corresponding right h eigenvector of the Christoffel matrix as g i = g i (x m, τ,u ), Γ il g l = Gg i, (16) the corresponding left h eigenvector of the Christoffel matrix as g i = g i (x m, τ,u ), g i Γ il = g l G. (17) The three right h eigenvectors of the Christoffel matrix the three left h eigenvectors of the Christoffel matrix are mutually biorthogonal. The lengths of the three right h eigenvectors of the Christoffel matrix are not determined, but we choose the lengths of the corresponding left h eigenvectors so that the three right h eigenvectors the three left h eigenvectors are mutually biorthonormal, g g = 1. (18) Stud. Geophys. Geod., 62 (2018) 263
4 L. Klimeš The zero order vectorial amplitude then reads U [0] i = U [0] g i, (19) where the zero order ray theory amplitude U [0] will be determined by the transport equation in Section 6. Whereas the two S waves, which propagate with different velocities, are linearly polarized in elastic media, they may be elliptically or even circularly polarized in viscoelastic media. Whereas the two elliptically polarized S waves always display equal hedness for a symmetric stiffness matrix, they display opposite hedness for a sufficiently non symmetric stiffness matrix. In order to satisfy Christoffel equation (15), selected eigenvalue G must be unit, G(x m, τ,u ) = 1. (20) Nonlinear first order partial differential equation (20) for travel time τ is called the Hamilton Jacobi equation. In wave propagation problems, it is also often referred to as the eional equation. The methods for solving the Hamilton Jacobi equation are already mostly developed (Hamilton, 1837; Červený, 1972; Klimeš, 2002; 2010; 2016). Hamilton Jacobi equation (20) generates the equations of rays (Hamilton equations, equations of geodesics) the related equations such as the Hamiltonian equations of geodesic deviation (dynamic ray tracing equations). The equations of rays may be expressed in various forms, e.g., as dx i dγ = 1 G, (21) 2 p i dp i dγ = 1 G 2 x i, (22) where parameter γ along rays coincides with the values of travel time τ. Refer to Section 4 for the case of complex valued eigenvalue G = G(x m, p u ) parameter γ. G Hereinafter, x G i p i denote the partial derivatives of function G(x m, p u ) of six phase space coordinates x m, p u. Using this notation, the partial derivatives of any function G ( x m, τ,u (x a ) ) of three spatial coordinates read Since G,j = G x j + G p τ,j. (23) G(x m, p u ) = g r (x m, p u )Γ rs (x a, p b )g s (x c, p d ), (24) the first order phase space derivatives of eigenvalue G(x m, p u ) can be expressed as G p i (x m, p u ) = g r (x m, p u ) Γrs p i (x a, p b )g s (x c, p d ) (25) G x i (xm, p u ) = g r (x m, p u ) Γrs x i (xa, p b )g s (x c, p d ). (26) 264 Stud. Geophys. Geod., 62 (2018)
5 Frequency-domain ray series for viscoelastic waves with a non-symmetric stiffness matrix Note that the phase space derivatives of the left h right h eigenvectors in the above derivatives (25) (26) of relation (24) vanish because of definitions (16) (17) normalization condition (18). We insert definition (12) into phase space derivatives (25) (26), obtain expressions Note that G p i (x m, p u ) = g r (x m, p u ) [ a ris (x a ) + a ris (x a ) ] p g s (x c, p d ) (27) G x i(xm, p u ) = g r (x m, p u ) arjs x i (x a )p j p g s (x c, p d ). (28) V i (x m, p u ) = 1 2 g r(x m, p u ) [ a ris (x a ) + a ris (x a ) ] p g s (x c, p d ) (29) is the ray velocity vector, which is sometimes also referred to as the group velocity vector, e.g., by Červený (2001). 4. COMPLEX VALUED AND REAL VALUED RAYS The above relations (20) (29) are applicable both to real valued travel time τ complex valued travel time τ. The methods for their solution are obvious if eigenvalue G = G(x m, p u ) in Hamilton Jacobi equation (20) is a real valued function of real valued slowness vector p u. However, eigenvalue G = G(x m, p u ) in Hamilton Jacobi equation (20) may be a complex valued function may generate complex valued rays. Since we usually do not now the values of the frequency domain stiffness tensor c ijl = c ijl (x m, ω) at imaginary coordinates x m, we can instead trace the reference real valued rays which are close to the complex valued rays defined by the equations (21) (22) of rays, to calculate travel time τ = τ(x m ) by means of the perturbation expansion by Klimeš Klimeš (2011) along the reference real valued rays. In this case, the slowness vector p u = τ,u in definition (29) is also approximated by the perturbation expansion. 5. PRINCIPAL AND ADDITIONAL AMPLITUDE COMPONENTS Analogously to the zero order vectorial amplitude (19), we define the amplitude components with respect to the three right h eigenvectors of Christoffel matrix (12): eigenvector g i corresponding to the selected eigenvalue G other two eigenvectors gi corresponding to other two eigenvalues G. We decompose each vectorial amplitude coefficient U [n] i in series (7) into the principal amplitude component U [n] two additional amplitude components U [n], U [n] i = U [n] g i + U [n] g i, (30) Stud. Geophys. Geod., 62 (2018) 265
6 L. Klimeš where the summation is performed over two superscripts. Considering expression (19), we assume that both U [0] = 0. Definition (9) with decomposition (30) yields N i( U m [n] ), τ,u = U [n] ( G 1 ) gi. (31) We multiply viscoelastodynamic equation (14) by left h eigenvector g i Christoffel matrix, consider relation (31), obtain expression of the U [n] = 1[ g i Mi( U [n 1] ), τ,u + g i L i( U [n 2] ) ] ( G 1 ) 1 (32) for the additional amplitude components in terms of lower order amplitudes. Expression (32) differs from the analogous expression for the symmetric stiffness matrix (Červený, 2001, Eq ) just by left h eigenvectors g i. 6. TRANSPORT EQUATION We multiply viscoelastodynamic equation (14) by left h eigenvector g i of the Christoffel matrix, consider relation (31), obtain transport equation g i M i( U [n], τ ),u + gi L i( U [n 1] ) = 0 (33) for the principal amplitude components. We separate the terms with higher order principal amplitude components from the terms containing higher order additional amplitude components lower order amplitude components, g i M i( U [n] ) g, τ,l = g i M i( U [n] g ), τ,l gi L i( U [n 1] ). (34) Relations (33) (34) differ from the analogous relations for the symmetric stiffness matrix (Červený, 2001, Eqs ) just by left h eigenvector g i. We express the left h side of transport equation (34) as g i M i (U [n] g m, τ,l ) = 2 V j U [n],j + ( V j ),j U [n] 2 S U [n], (35) where ray velocity vector V j is given by definition (29), which can also be expressed as (x a )V i (x m, p u ) = 1 2 g r(x m, p u ) [ c ris (x a ) + c ris (x a ) ] p g s (x c, p d ). (36) The first two terms on the right h side of relation (35) are well nown from the ray series with a symmetric stiffness matrix (Červený, 2001, Eq ). Quantity S in the rightmost term of transport equation (35) can be determined using definition (10) as S = 1 { 1 ( [ g i c ijl + c ijl) ] } τ, g l 2 2 g ( i c ijl ) τ, g l,j,j g ic ijl τ, g l,j, (37) see (36). Considering (35), we express transport equation (34) as V j U [n],j ),j U [n] = S U [n] + Z [n 1], (38) 266 Stud. Geophys. Geod., 62 (2018)
7 Frequency-domain ray series for viscoelastic waves with a non-symmetric stiffness matrix where Z [n 1] = 1 [ 2 g i M i( U [n] g, τ ),u + gi L i( U [n 1] ) ]. (39) The solution of transport equation (38) for n = 0 reads ( τ ) U [0] = U [0] 0 ( 0 J 0 ) 1 2 ( J) 1 2 exp dγ S τ 0, (40) where subscript 0 denotes the initial conditions. Squared geometrical spreading ( ) x i J = det γ a (41) (Babich, 1961, Eq. 3.7; Červený, 2001, Eq ) represents the Jacobian of transformation from ray coordinates γ 1, γ 2, γ 3 to spatial coordinates x i. Here γ 1 γ 2 are the ray parameters, γ 3 = γ. Factor exp ( τ τ 0 dγ S ) in (40) is present due to the sew part of the stiffness matrix, see the next section. The solution of transport equation (38) for n>0 reads (Červený, 2001, Eq ) [ U [n] = U [0] U [n] τ ] 0 Z [n 1] + dγ U [0] U [0] 0 τ 0. (42) 7. NON RECIPROCITY DUE TO A NON SYMMETRIC STIFFNESS MATRIX The only difference of expressions (40) (42) for the principal amplitudes from the analogous expressions derived for a symmetric stiffness matrix is exponential term exp ( τ τ 0 dγ S ) in expression (40). We shall now derive various expressions for quantity S = S(x m ). We express definition (37) as S = 1 { 1 ( [ g i c ijl +c ijl) ] } τ, g l 2 2 ( g i c ijl ) τ, g l,j,j + g i,jc ijl τ, g l g i c ijl τ, g l,j. (43) After summation, we obtain expression S = 1 2 ] [ g i,j a ijl τ, g l g i a ijl τ, g l,j 1 ( [ g i a ijl a ijl) ] τ, g l 4,j, (44) where a ijl is the density reduced stiffness tensor given by definition (13). The last term on the right h side differs from 1 j ( V ) 2,j = 1 ( [ g i a ijl + a ijl) ] τ, g l (45) 4,j Stud. Geophys. Geod., 62 (2018) 267
8 L. Klimeš just by the subtraction. Note also that [ ( a ijl a ijl) τ, ] τ,j = 0. (46) We differentiate the product in the rightmost term of expression (44) obtain expression S = 1 4 g i,j( a ijl +a ijl) τ, g l 1 4 g ( i a ijl +a ijl) τ, g l,j 1 [ 4 g i (a ijl a ijl)] τ,g l.,j (47) We differentiate characteristic equation (16) for the right h eigenvector with respect to spatial coordinates, consider the biorthonormality of the left h right h eigenvectors, obtain relation g i,j = g i g Γ l,jg l ( G G ) 1 + gi g g,j (48) for the spatial gradient of the right h eigenvector. The rightmost term in relation (48) accounts for the undefined changes of the length of the right h eigenvector g i. We differentiate characteristic equation (17) for the left h eigenvector with respect to spatial coordinates, obtain analogous relation g i,j = g Γ l,jg l g i ( G G ) 1 + gi g g,j (49) for the spatial gradient of the left h eigenvector. The rightmost term in relation (49) accounts for the undefined changes of the length of the left h eigenvector g i, satisfies identity g g,j = g g,j (50) obtained by differentiating normalization condition (18). We insert the gradients (48) (49) of the eigenvectors of the Christoffel matrix into expression (47), consider identity (50) arrive at expression S = 1 4 ( g Γ l,jg l g r Γ rs g s g r Γ rs g s g Γ l,jg l ) ( G G ) g ( i c ijl c ijl),j τ,g l g g,j V j, (51) where ray velocity vector V j is given by definition (29). We now express the spatial derivatives Γ l,j of the Christoffel matrix in terms of its phase space derivatives as Γ l,j = Γl x j + Γl p s τ,sj. (52) The partial derivatives of the Christoffel matrix (12) with respect to phase space coordinates x m p u read Γ l x j (xm, τ,u ) = a rsl,j τ,r τ,r (53) 268 Stud. Geophys. Geod., 62 (2018)
9 Frequency-domain ray series for viscoelastic waves with a non-symmetric stiffness matrix Γ l (x m, τ,u ) = ( a jrl + a rjl) τ,r. (54) Quantity (51) with identity (52) finally reads S = 1 Γ ( g l 4 x j g l g r Γ rs Γ l g s g gl g r 1 4 g ( i c ijl c ijl),j τ dg i,g l g i dγ The last term g i dg i dγ Γ rs ) ( x j g s G G ) 1. (55) in expression (55) represents just the correction of principal amplitude U [n] in decomposition (30) due to the undefined length of right h eigenvector g i, vanish if we put g i dg i dγ = 0 (56) along each ray. Each element of matrix (c ijl c ijl ),j τ, in the last but one term of relation (55) represents the divergence of a vector tangent to the wavefront, see identity (46). Note also that g Γ l g l τ,j = 0 (57) Vectors g Γ l g l τ,j = 0. (58) Γ l g gl (59) g Γ l g l (60) thus represent two sets of the contravariant basis vectors of the ray centred coordinate system. Expression (55) for quantity S may be singular at slowness surface singularities, but is regular at spatial caustics. Quantity S vanishes for a symmetric stiffness matrix. For a non symmetric stiffness matrix, quantity S vanishes in a homogeneous medium. Quantity S which expresses the difference between the ray series for symmetric non symmetric stiffness matrices is thus generated by a combination of a non symmetric stiffness matrix heterogeneities, is generally non zero in a heterogeneous medium with the non symmetric stiffness matrix. Expression (55) with identities (46), (57) (58) suggest that quantity S may mostly be influenced by the wavefront tangent component of the gradient of the sew part of the stiffness matrix. Stud. Geophys. Geod., 62 (2018) 269
10 L. Klimeš 8. CONCLUSIONS We have derived the anisotropic ray theory series for viscoelastic waves with a non-symmetric stiffness matrix. These ray series enable us to estimate which phenomena could be observed in the wave field if the stiffness matrix were non symmetric. Whereas the two S waves, which propagate with different velocities, are linearly polarized in elastic media, they may be elliptically or even circularly polarized in viscoelastic media. Whereas the two elliptically polarized S waves always display equal hedness for a symmetric stiffness matrix, they display opposite hedness for a sufficiently non symmetric stiffness matrix, similarly as electromagnetic waves in optically active media. The ray theory amplitudes corresponding to a non symmetric stiffness matrix are not reciprocal in the same way as the ray theory amplitudes corresponding to a symmetric stiffness matrix. This non reciprocity is expressed in terms of quantity (55) in the expression (40) for the zero order ray theory amplitude. Refer to Klimeš (2017, Eq. 18) for the sense in which the ray theory Green function corresponding to a non symmetric stiffness matrix is reciprocal. Acnowledgements: The suggestions by Einar Iversen, Edward S. Krebes two anonymous reviewers made it possible for me to improve the paper considerably. The research has been supported by the Grant Agency of the Czech Republic under contract S, by the members of the consortium Seismic Waves in Complex 3 D Structures (see ). References Babich V.M., Ray method of calculating the intensity of wavefronts in the case of a heterogeneous, anisotropic, elastic medium. In: Petrashen G.I. (Ed.), Problems of the Dynamic Theory of Propagation of Seismic Waves, Vol. 5. Leningrad Univ. Press, Leningrad, (in Russian, English translation: Geophys. J. Int., 118(1994), ). Červený V., Seismic rays ray intensities in inhomogeneous anisotropic media. Geophys. J. R. Astr. Soc., 29, Červený V., Seismic Ray Theory. Cambridge University Press, Cambridge, U.K. Červený V., Klimeš L. Pšenčí I., Seismic ray method: Recent developments. Adv. Geophys., 48, Chapman C.H., Fundamentals of Seismic Wave Propagation. Cambridge University Press, Cambridge, U.K. Hamilton W.R., Third supplement to an essay on the theory of systems of rays. Trans. Roy. Irish Acad., 17, Klimeš L., Second order higher order perturbations of travel time in isotropic anisotropic media. Stud. Geophys. Geod., 46, Stud. Geophys. Geod., 62 (2018)
11 Frequency-domain ray series for viscoelastic waves with a non-symmetric stiffness matrix Klimeš L., Transformation of spatial perturbation derivatives of travel time at a general interface between two general media. Seismic Waves in Complex 3 D Structures, 20, ( Klimeš L., Coupling ray series. Stud. Geophys. Geod., 57, Klimeš L., Transformation of spatial perturbation derivatives of travel time at a curved interface between two arbitrary media. Stud. Geophys. Geod., 60, Klimeš L., Representation theorem for viscoelastic waves with a non symmetric stiffness matrix. Seismic Waves in Complex 3 D Structures, 27, ( Klimeš M. Klimeš L., Perturbation expansions of complex valued traveltime along real valued reference rays. Geophys. J. Int., 186, Thomson C.J., Complex rays wave pacets for decaying signals in inhomogeneous, anisotropic anelastic media. Stud. Geophys. Geod., 41, Stud. Geophys. Geod., 62 (2018) 271
Superpositions of Gaussian beams and column Gaussian packets in heterogeneous anisotropic media
Superpositions of Gaussian beams and column Gaussian packets in heterogeneous anisotropic media Luděk Klimeš Department of Geophysics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu,
More informationPrevailing-frequency approximation of the coupling ray theory for electromagnetic waves or elastic S waves
Prevailing-frequency approximation of the coupling ray theory for electromagnetic waves or elastic S waves Luděk Klimeš and Petr Bulant Department of Geophysics, Faculty of Mathematics and Physics, Charles
More informationGaussian beams in inhomogeneous anisotropic layered structures
Gaussian beams in inhomogeneous anisotropic layered structures Vlastislav Červený 1 ) and Ivan Pšenčík 2 ) 1) Charles University, Faculty of Mathematics and Physics, Department of Geophysics, Ke Karlovu
More informationPrevailing-frequency approximation of the coupling ray theory for S waves
Prevailing-frequency approximation of the coupling ray theory for S waves Petr Bulant & Luděk Klimeš Department of Geophysics Faculty of Mathematics and Physics Charles University in Prague S EI S MIC
More informationRelation between the propagator matrix of geodesic deviation and the second-order derivatives of the characteristic function
Journal of Electromagnetic Waves an Applications 203 Vol. 27 No. 3 589 60 http://x.oi.org/0.080/0920507.203.808595 Relation between the propagator matrix of geoesic eviation an the secon-orer erivatives
More informationKirchhoff prestack depth migration in velocity models with and without rotation of the tensor of elastic moduli: Orthorhombic and triclinic anisotropy
Kirchhoff prestack depth migration in velocity models with and without rotation of the tensor of elastic moduli: Orthorhombic and triclinic anisotropy Václav Bucha Department of Geophysics, Faculty of
More informationSeismic ray theory. Summary. Introduction
Seismic ray theory Vlastislav Červený 1 ) and Ivan Pšenčík 2 ) 1 ) Charles University, Faculty of Mathematics and Physics, Department of Geophysics, Prague, Czech Republic. E-mail: vcerveny@seis.karlov.mff.cuni.cz.
More informationPARTICLE MOTION OF PLANE WAVES IN VISCOELASTIC ANISOTROPIC MEDIA
Geologiya and Geophysics i Geofizika Vol. 47, No. 5, pp. 557-567, 26 UDC 55.834 PARTICLE MOTION OF PLANE WAVES IN VISCOELASTIC ANISOTROPIC MEDIA V. C v erveny and I. Ps v enc v ík* Department of Geophysics,
More informationSeismic Waves in Complex 3 D Structures, 26 (2016), (ISSN , online at
Kirchhoff prestack depth migration in simple orthorhombic and triclinic models with differently rotated elasticity tensor: comparison with zero-offset travel-time perturbations Václav Bucha Department
More informationVáclav Bucha. Department of Geophysics Faculty of Mathematics and Physics Charles University in Prague. SW3D meeting June 6-7, 2016 C OM S TR 3 D
Kirchhoff prestack depth migration in simple orthorhombic and triclinic models with differently rotated elasticity tensor: comparison with zero-offset travel-time perturbations Václav Bucha Department
More informationMoveout approximation for P waves in a homogeneous VTI medium
Moveout approximation for P waves in a homogeneous VTI medium Véronique Farra 1 and Ivan Pšenčík 2 1 Institut de Physique du Globe de Paris, Sorbonne Paris Cité, Université Paris Diderot, UMR 7154 CNRS,
More informationRelationships between the velocities and the elastic constants of an anisotropic solid possessing orthorhombic symmetry
Relationships between the velocities and the elastic constants of an anisotropic solid possessing orthorhombic symmetry R. James Brown ABSTRACT This paper reviews the equations of body-wave propagation
More informationThe propagation of seismic body waves in complex, laterally varying 3-D layered
CHAPTER ONE Introduction The propagation of seismic body waves in complex, laterally varying 3-D layered structures is a complicated process. Analytical solutions of the elastodynamic equations for such
More informationFar-field radiation from seismic sources in 2D attenuative anisotropic media
CWP-535 Far-field radiation from seismic sources in 2D attenuative anisotropic media Yaping Zhu and Ilya Tsvankin Center for Wave Phenomena, Department of Geophysics, Colorado School of Mines, Golden,
More informationAcoustic axes in triclinic anisotropy
Acoustic axes in triclinic anisotropy Václav Vavryčuk a Geophysical Institute, Academy of Sciences, Boční II/40, 4 3 Praha 4, Czech Republic Received 2 November 2004; revised 9 May 2005; accepted 23 May
More informationReceiver. Johana Brokešová Charles University in Prague
Propagation of seismic waves - theoretical background Receiver Johana Brokešová Charles University in Prague Seismic waves = waves in elastic continuum a model of the medium through which the waves propagate
More informationKirchhoff prestack depth migration in simple models of various anisotropy
Kirchhoff prestack depth migration in simple models of various anisotropy Václav Bucha Department of Geophysics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu, 6 Praha, Czech Republic,
More informationLongitudinal waves of a finite amplitude in nonlinear elastic media
Longitudinal waves of a finite amplitude in nonlinear elastic media Ivan A. Molotkov 1), Ivan Pšenčík 2) 1) IZMIRAN, Troitsk, Moscow Region, 142190, Russia. E-mail molotkov@izmiran.rssi.ru 2) Geophysical
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 informationReal ray tracing in anisotropic viscoelastic media
Geophys. J. nt. (28) 175, 617 626 doi: 1.1111/j.1365-246X.28.3898.x eal ray tracing in anisotropic viscoelastic media Václav Vavryčuk nstitute of Geophysics, Academy of Sciences, Boční /141, Praha 4, Czech
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 informationPoint-source radiation in attenuative anisotropic media
GEOPHYSICS, VOL. 79, NO. 5 (SEPTEMBER-OCTOBER 4); P. WB5 WB34, 4 FIGS., TABLE..9/GEO3-47. Downloaded /7/4 to 38.67..93. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/
More informationIntroduction to Seismology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 1.510 Introduction to Seismology Spring 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 1.510 Introduction to
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationKirchhoff prestack depth migration in simple models with differently rotated elasticity tensor: orthorhombic and triclinic anisotropy
Kirchhoff prestack depth migration in simple models with differently rotated elasticity tensor: orthorhombic and triclinic anisotropy Václav Bucha Department of Geophysics, Faculty of Mathematics and Physics,
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 informationProperties of S waves near a kiss singularity: a comparison of exact and ray solutions
Geophys. J. Int. (1999) 138, 581 589 Properties of S waves near a kiss singularity: a comparison of exact and ray solutions Václav Vavryčuk Geophysical Institute, Academy of Sciences of the Czech Republic,
More informationSnell s law in transversely isotropic media using linearized group velocities and related quantities
Snell's law using group angles and velocities Snell s law in transversely isotropic media using linearized group velocities and related quantities P.F. Daley ABSTRACT Using a linearized approximation for
More informationRobert W. Vestrum and R. James Brown
Group versus phase velocity in laboratory measurements on orthorhombic samples Robert W. Vestrum and R. James Brown ABSTRACT In laboratory measurements of traveltimes across anisotropic materials, there
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 informationElements of Dynamic Ray tracing
Elements of Dynamic Ray tracing ErSE360 Contents List of Figures 1 1 Introduction Coordinates systems 1 Cartesian coordinates 3 Ray coordinates γ 1 γ s 3 3 Ray-centered coordinate system q 1 q s 3 3 Coordinate
More informationAcoustic axes in weak triclinic anisotropy
Geophys. J. Int. (2005) 163, 629 638 doi: 10.1111/j.1365-246X.2005.02762.x Acoustic axes in weak triclinic anisotropy Václav Vavryčuk Geophysical Institute, Academy of Sciences of the Czech Republic, Boční
More informationSymmetry and Properties of Crystals (MSE638) Stress and Strain Tensor
Symmetry and Properties of Crystals (MSE638) Stress and Strain Tensor Somnath Bhowmick Materials Science and Engineering, IIT Kanpur April 6, 2018 Tensile test and Hooke s Law Upto certain strain (0.75),
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 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 informationASYMPTOTIC RAY METHOD IN SEISMOLOGY A TUTORIAL
ASYMPTOTIC RAY METHOD IN SEISMOLOGY A TUTORIAL Johana Brokešová Matfyzpress 2006 Funded by the European Commission s Human Resources and Mobility Programme Marie Curie Research Training Network SPICE Contract
More informationDuality, Dual Variational Principles
Duality, Dual Variational Principles April 5, 2013 Contents 1 Duality 1 1.1 Legendre and Young-Fenchel transforms.............. 1 1.2 Second conjugate and convexification................ 4 1.3 Hamiltonian
More informationA Survey of Computational High Frequency Wave Propagation II. Olof Runborg NADA, KTH
A Survey of Computational High Frequency Wave Propagation II Olof Runborg NADA, KTH High Frequency Wave Propagation CSCAMM, September 19-22, 2005 Numerical methods Direct methods Wave equation (time domain)
More informationModern Optics Prof. Partha Roy Chaudhuri Department of Physics Indian Institute of Technology, Kharagpur
Modern Optics Prof. Partha Roy Chaudhuri Department of Physics Indian Institute of Technology, Kharagpur Lecture 09 Wave propagation in anisotropic media (Contd.) So, we have seen the various aspects of
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationSeismic ray theory. Chapter Introduction
Chapter 2 Seismic ray theory This chapter provides a background on seismic ray theory. The emphasis is placed on explanation and interpretation of the basic ingredients, rather than on the derivation of
More informationChap. 4. Electromagnetic Propagation in Anisotropic Media
Chap. 4. Electromagnetic Propagation in Anisotropic Media - Optical properties depend on the direction of propagation and the polarization of the light. - Crystals such as calcite, quartz, KDP, and liquid
More informationAdvanced Vibrations. Elements of Analytical Dynamics. By: H. Ahmadian Lecture One
Advanced Vibrations Lecture One Elements of Analytical Dynamics By: H. Ahmadian ahmadian@iust.ac.ir Elements of Analytical Dynamics Newton's laws were formulated for a single particle Can be extended to
More informationP137 Our Experiences of 3D Synthetic Seismic Modeling with Tip-wave Superposition Method and Effective Coefficients
P137 Our Experiences of 3D Synthetic Seismic Modeling with Tip-wave Superposition Method and Effective Coefficients M. Ayzenberg (StatoilHydro), A. Aizenberg (Institute of Petroleum Geology and Geophysics),
More informationChap. 1 Fundamental Concepts
NE 2 Chap. 1 Fundamental Concepts Important Laws in Electromagnetics Coulomb s Law (1785) Gauss s Law (1839) Ampere s Law (1827) Ohm s Law (1827) Kirchhoff s Law (1845) Biot-Savart Law (1820) Faradays
More informationEffects of Fracture Parameters in an Anisotropy Model on P-Wave Azimuthal Amplitude Responses
PROC. ITB Eng. Science Vol. 38 B, No. 2, 2006, 159-170 159 Effects of Fracture Parameters in an Anisotropy Model on P-Wave Azimuthal Amplitude Responses Fatkhan Program Studi Teknik Geofisika FIKTM-ITB
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationModern Optics Prof. Partha Roy Chaudhuri Department of Physics Indian Institute of Technology, Kharagpur
Modern Optics Prof. Partha Roy Chaudhuri Department of Physics Indian Institute of Technology, Kharagpur Lecture 08 Wave propagation in anisotropic media Now, we will discuss the propagation of electromagnetic
More informationINDEX 363. Cartesian coordinates 19,20,42, 67, 83 Cartesian tensors 84, 87, 226
INDEX 363 A Absolute differentiation 120 Absolute scalar field 43 Absolute tensor 45,46,47,48 Acceleration 121, 190, 192 Action integral 198 Addition of systems 6, 51 Addition of tensors 6, 51 Adherence
More informationWenyong Pan and Lianjie Huang. Los Alamos National Laboratory, Geophysics Group, MS D452, Los Alamos, NM 87545, USA
PROCEEDINGS, 44th Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 019 SGP-TR-14 Adaptive Viscoelastic-Waveform Inversion Using the Local Wavelet
More informationRay equations of a weak shock in a hyperbolic system of conservation laws in multi-dimensions
Proc. Indian Acad. Sci. (Math. Sci.) Vol. 126, No. 2, May 2016, pp. 199 206. c Indian Academy of Sciences Ray equations of a weak in a hyperbolic system of conservation laws in multi-dimensions PHOOLAN
More informationAzimuthal AVO and Curvature. Jesse Kolb* David Cho Kris Innanen
Azimuthal AVO and Curvature Jesse Kolb* David Cho Kris Innanen Azimuthal AVO What is azimuthal AVO?: Analysis of incidence angle and azimuthal amplitude variations of reflection coefficients; Measures
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 informationNumerical methods to compute optical errors due to stress birefringence
Numerical methods to compute optical errors due to stress birefringence Keith B. Doyle Optical Research Associates, 8 West Park Drive, Westborough, MA Victor L. Genberg & Gregory J. Michels Sigmadyne,
More informationChapter 0. Preliminaries. 0.1 Things you should already know
Chapter 0 Preliminaries These notes cover the course MATH45061 (Continuum Mechanics) and are intended to supplement the lectures. The course does not follow any particular text, so you do not need to buy
More informationULTRASONIC WAVE PROPAGATION IN DISSIMILAR METAL WELDS APPLICATION OF A RAY-BASED MODEL AND COMPARISON WITH EXPERIMENTAL RESULTS
ULTRASONIC WAVE PROPAGATION IN DISSIMILAR METAL WELDS APPLICATION OF A RAY-BASED MODEL AND COMPARISON WITH EXPERIMENTAL RESULTS Audrey GARDAHAUT 1, Hugues LOURME 1, Frédéric JENSON 1, Shan LIN 2, Masaki
More informationAmplitude calculations for 3-D Gaussian beam migration using complex-valued traveltimes
CWP-648 Amplitude calculations for 3-D Gaussian beam migration using complex-valued traveltimes Norman Bleistein 1 & S. H. Gray 2 1 Center for Wave Phenomena, Dept. of Geophysics, Colo. Sch. of Mines,
More informationPropagation of inhomogeneous plane waves in monoclinic crystals subject to initial fields
Annals of the University of Bucharest (mathematical series) 5 (LXIII) (014), 367 381 Propagation of inhomogeneous plane waves in monoclinic crystals subject to initial fields Olivian Simionescu-Panait
More informationExtended isochron rays in prestack depth (map) migration
Extended isochron rays in prestack depth (map) migration A.A. Duchkov and M.V. de Hoop Purdue University, 150 N.University st., West Lafayette, IN, 47907 e-mail: aduchkov@purdue.edu (December 15, 2008)
More information1. Reminder: E-Dynamics in homogenous media and at interfaces
0. Introduction 1. Reminder: E-Dynamics in homogenous media and at interfaces 2. Photonic Crystals 2.1 Introduction 2.2 1D Photonic Crystals 2.3 2D and 3D Photonic Crystals 2.4 Numerical Methods 2.4.1
More informationThe signature of attenuation and anisotropy on AVO and inversion sensitivities
The signature of attenuation and anisotropy on AVO and inversion sensitivities Shahpoor Moradi and Kristopher Innanen University of Calgary, Department of Geoscience, Calgary, Canada Summary We study the
More informationIranian Journal of Mathematical Sciences and Informatics Vol.2, No.2 (2007), pp 1-16
Iranian Journal of Mathematical Sciences and Informatics Vol.2, No.2 (2007), pp 1-16 THE EFFECT OF PURE SHEAR ON THE REFLECTION OF PLANE WAVES AT THE BOUNDARY OF AN ELASTIC HALF-SPACE W. HUSSAIN DEPARTMENT
More informationA direct proof of Malus theorem using the symplectic sructure of the set of oriented straight lines
A direct proof of Malus theorem using the symplectic sructure of the set of oriented straight lines Charles-Michel Marle Université Pierre et Marie Curie 58-th Souriau Colloquium Aix en Provence, July
More informationThe Basic Properties of Surface Waves
The Basic Properties of Surface Waves Lapo Boschi lapo@erdw.ethz.ch April 24, 202 Love and Rayleigh Waves Whenever an elastic medium is bounded by a free surface, coherent waves arise that travel along
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 informationElements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004
Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic
More informationWeek 6: Differential geometry I
Week 6: Differential geometry I Tensor algebra Covariant and contravariant tensors Consider two n dimensional coordinate systems x and x and assume that we can express the x i as functions of the x i,
More informationIntroduction to Polarization
Phone: Ext 659, E-mail: hcchui@mail.ncku.edu.tw Fall/007 Introduction to Polarization Text Book: A Yariv and P Yeh, Photonics, Oxford (007) 1.6 Polarization States and Representations (Stokes Parameters
More informationConstitutive Equations
Constitutive quations David Roylance Department of Materials Science and ngineering Massachusetts Institute of Technology Cambridge, MA 0239 October 4, 2000 Introduction The modules on kinematics (Module
More informationDistances, volumes, and integration
Distances, volumes, and integration Scribe: Aric Bartle 1 Local Shape of a Surface A question that we may ask ourselves is what significance does the second fundamental form play in the geometric characteristics
More informationIn a uniform 3D medium, we have seen that the acoustic Green s function (propagator) is
Chapter Geometrical optics The material in this chapter is not needed for SAR or CT, but it is foundational for seismic imaging. For simplicity, in this chapter we study the variable-wave speed wave equation
More informationLawrence Berkeley National Laboratory
Lawrence Berkeley National Laboratory Peer Reviewed Title: Fracture permeability and seismic wave scattering--poroelastic linear-slip interface model for heterogeneous fractures Author: Nakagawa, S. Publication
More informationLecture 4: Anisotropic Media. Dichroism. Optical Activity. Faraday Effect in Transparent Media. Stress Birefringence. Form Birefringence
Lecture 4: Anisotropic Media Outline Dichroism Optical Activity 3 Faraday Effect in Transparent Media 4 Stress Birefringence 5 Form Birefringence 6 Electro-Optics Dichroism some materials exhibit different
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationRadiation pattern in homogeneous and transversely isotropic attenuating media
Radiation pattern in homogeneous and transversely isotropic attenuating media Satish Sinha*, Sergey Abaseyev** and Evgeni Chesnokov** *Rajiv Gandhi Institute of Petroleum Technology, Rae Bareli, UP 229010
More informationOn divergence representations of the Gaussian and the mean curvature of surfaces and applications
Bull. Nov. Comp. Center, Math. Model. in Geoph., 17 (014), 35 45 c 014 NCC Publisher On divergence representations of the Gaussian and the mean curvature of surfaces and applications A.G. Megrabov Abstract.
More informationRecurrences and Full-revivals in Quantum Walks
Recurrences and Full-revivals in Quantum Walks M. Štefaňák (1), I. Jex (1), T. Kiss (2) (1) Department of Physics, Faculty of Nuclear Sciences and Physical Engineering, Czech Technical University in Prague,
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 informationOn Transversely Isotropic Elastic Media with Ellipsoidal Slowness Surfaces
On Transversely Isotropic Elastic Media with Ellipsoidal Slowness Surfaces Anna L. Mazzucato a,,1, Lizabeth V. Rachele b,2 a Department of Mathematics, Pennsylvania State University, University Par, PA
More informationEngineering Sciences 241 Advanced Elasticity, Spring Distributed Thursday 8 February.
Engineering Sciences 241 Advanced Elasticity, Spring 2001 J. R. Rice Homework Problems / Class Notes Mechanics of finite deformation (list of references at end) Distributed Thursday 8 February. Problems
More informationSHEAR-WAVE REFLECTION MOVEOUT FOR AZIMUTHALLY ANISOTROPIC MEDIA
SHEAR-WAVE REFLECTION MOVEOUT FOR AZIMUTHALLY ANISOTROPIC MEDIA AbdulFattah AI-Dajani and M. Nafi Toksoz Earth Resources Laboratory Department of Earth, Atmospheric, and Planetary Sciences Massachusetts
More informationTime-to-depth conversion and seismic velocity estimation using time-migration velocity a
Time-to-depth conversion and seismic velocity estimation using time-migration velocity a a Published in Geophysics, 73, no. 5, VE205-VE210, (2008) Maria Cameron, Sergey Fomel, and James Sethian ABSTRACT
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (04) 98 653 66 GJI Seismology doi: 0.093/gji/ggu9 Effective orthorhombic anisotropic models for wavefield extrapolation Wilson Ibanez-Jacome Tariq Alkhalifah
More informationarxiv: v1 [physics.geo-ph] 31 Dec 2018
Generalized Dix equation and analytic treatment of normal-moveout velocity for anisotropic media Vladimir Grechka 1,, Ilya Tsvankin 1, and Jack K Cohen 1 1 Center for Wave Phenomena, Colorado School of
More informationTable of Contents. Preface... 13
Table of Contents Preface... 13 Chapter 1. Vibrations of Continuous Elastic Solid Media... 17 1.1. Objective of the chapter... 17 1.2. Equations of motion and boundary conditions of continuous media...
More informationWave Propagation in Uniaxial Media. Reflection and Transmission at Interfaces
Lecture 5: Crystal Optics Outline 1 Homogeneous, Anisotropic Media 2 Crystals 3 Plane Waves in Anisotropic Media 4 Wave Propagation in Uniaxial Media 5 Reflection and Transmission at Interfaces Christoph
More informationChapter 3 Stress, Strain, Virtual Power and Conservation Principles
Chapter 3 Stress, Strain, irtual Power and Conservation Principles 1 Introduction Stress and strain are key concepts in the analytical characterization of the mechanical state of a solid body. While stress
More informationReflection of SV- Waves from the Free Surface of a. Magneto-Thermoelastic Isotropic Elastic. Half-Space under Initial Stress
Mathematica Aeterna, Vol. 4, 4, no. 8, 877-93 Reflection of SV- Waves from the Free Surface of a Magneto-Thermoelastic Isotropic Elastic Half-Space under Initial Stress Rajneesh Kakar Faculty of Engineering
More informationA Survey of Computational High Frequency Wave Propagation I
A Survey of Computational High Frequency Wave Propagation I Bjorn Engquist University of Texas at Austin CSCAMM Workshop on High Frequency Wave Propagation, University of Maryland, September 19-22, 2005
More informationTWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY. P. Yu 1,2 and M. Han 1
COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 3, Number 3, September 2004 pp. 515 526 TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY P. Yu 1,2
More informationGoing with the flow: A study of Lagrangian derivatives
1 Going with the flow: A study of Lagrangian derivatives Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc/ 12 February
More informationConstitutive model and wave equations for linear, viscoelastic, anisotropic media
GEOPHYSICS, VOL. 60, NO. 2 (MARCH-APRIL 995); P. 537-548, 5 FIGS., TABLE. Constitutive model and wave equations for linear, viscoelastic, anisotropic media Jose M. Carcione* ABSTRACT Rocks are far from
More informationArchivum Mathematicum
Archivum Mathematicum Zdeněk Dušek; Oldřich Kowalski How many are affine connections with torsion Archivum Mathematicum, Vol. 50 (2014), No. 5, 257 264 Persistent URL: http://dml.cz/dmlcz/144068 Terms
More informationPart III. Interaction with Single Electrons - Plane Wave Orbits
Part III - Orbits 52 / 115 3 Motion of an Electron in an Electromagnetic 53 / 115 Single electron motion in EM plane wave Electron momentum in electromagnetic wave with fields E and B given by Lorentz
More informationPhysics Letters A 374 (2010) Contents lists available at ScienceDirect. Physics Letters A.
Physics Letters A 374 (2010) 1063 1067 Contents lists available at ScienceDirect Physics Letters A www.elsevier.com/locate/pla Macroscopic far-field observation of the sub-wavelength near-field dipole
More informationMeasurement of deformation. Measurement of elastic force. Constitutive law. Finite element method
Deformable Bodies Deformation x p(x) Given a rest shape x and its deformed configuration p(x), how large is the internal restoring force f(p)? To answer this question, we need a way to measure deformation
More informationPart 5 ACOUSTIC WAVE PROPAGATION IN ANISOTROPIC MEDIA
Part 5 ACOUSTIC WAVE PROPAGATION IN ANISOTROPIC MEDIA Review of Fundamentals displacement-strain relation stress-strain relation balance of momentum (deformation) (constitutive equation) (Newton's Law)
More informationLecture 3: Asymptotic Methods for the Reduced Wave Equation
Lecture 3: Asymptotic Methods for the Reduced Wave Equation Joseph B. Keller 1 The Reduced Wave Equation Let v (t,) satisfy the wave equation v 1 2 v c 2 = 0, (1) (X) t2 where c(x) is the propagation speed
More informationSymmetry of the Dielectric Tensor
Symmetry of the Dielectric Tensor Curtis R. Menyuk June 11, 2010 In this note, I derive the symmetry of the dielectric tensor in two ways. The derivations are taken from Landau and Lifshitz s Statistical
More informationTitle: A note on perturbation formulae for the surface-wave speed due to perturbations in material properties
Editorial Manager(tm) for Journal of Elasticity Manuscript Draft Manuscript Number: ELAS248R2 Title: A note on perturbation formulae for the surface-wave speed due to perturbations in material properties
More information