energyVector edit page

Calculates Energy velocity vector (km/s)

Description

Energy velocity for lossless elastic medium (i.e. no attenuation) Good proxy for group velocity, which typically has some energy loss The formula is given by F.I. Fedorov(1968)Theory of Elastic Waves in Crystals, 375 pp. New York: Penum Press.

Ve_i = C_ijkl P_j P_l X_k / rho*V

N.B. E_magnitude should be equal or more than plane wave velocity vp, vs1 or vs2

David Mainprice 6/02/2018

Syntax

E = energyVector(C,x,v,p)
E = energyVector(C,x,v,p,rho)
E = energyVector(C,[],vFun,pFun)
E = energyVector(C,[],vFun,pFun,rho)

Input

C stiffnessTensor (units GPa)
x vector3d propagation direction
v plane wave velocity (unit km/s) e.g. vp,vs1 or vs2
p vector3d plane wave polarization vector e.g. pp,ps1 or ps2
vFun S2Fun plane wave velocity (unit km/s) e.g. vp,vs1 or vs2
pFun S2AxisField plane wave polarization vector e.g. pp,ps1 or ps2
rho density in g/cm3

Output

E Energy velocity vector (units km/s)