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) |