Parallel Distributed Simulation of Elastic Waves

Page 57

49 Tyy(i + 1/2,j, k) = (2p +A) uh [v(i + 1/2,j + 1/2, k)- v(i + 1/2,j- 1/2, k)] - 2~h [V(i + 1/2, j + 3/2, k)- v(i + 1/2, j- 3/2, k)J) +A ( 8 9 h [u(i + I, j, k)- u(i, j, k)] - 2~h [u(i + 2,j, k)- u(i- 1,j, k)] + 8 9 h [w(i + 112, j, k + 112) - w(i + 112, j, k- 112)] - 2~h [w(i + 1/2, j, k + 3/2)- w(i + 1/2, j, k- 3/2)]) (A.5) Tzz ( i + 1/2, j, k) = (2p +A) uh [w(i + 1/2, j, k + 1/2)- w(i + 1/2, j, k- 1/2)] - 2~h [w(i + 1/2, j, k + 3/2)- w(i + 1/2, j, k- 3/2)]) +A ( Bgh [u(i + I, j, k)- u(i, j, k)] 1 - 24 h [u(i + 2,j, k)- u(i- 1,j, k)] + 8~1 [ v ( i + 1 I 2, i + 1 I 2, k) - v ( i + 1 12, i - 1 12, k) J - 2~h [v(i + l/2,j + 3/2, k)- v(i + l/2,j- 3/2, k)J) (A.6) Try(i, j + 1/2, k) = p ( 8 9 h [u(i, j + I, k)- u(i, j, k)] 1 - 24 h [u(i,j + 2, k)- u(i,j- 1, k)] 9 + Bh [v(i + 1l2,j + 112, k)- v(i- 1l2,j + 112, k)] - 2~h [v(i + 3/2, j + 1/2, k)- v(i- 3/2, j + 1/2, k)J) (A. 7)