Parallel Distributed Simulation of Elastic Waves

Page 56

48 - 2 !h [Tzy(i+ 1l2,j+ 1l2,k+312) -Tzy(i + 1l2,j + 112, k- 312)) (A.2) 82 P w(i+112,j,k+112,n) P [ (. I . I ) (. 1 ·k 1 ) 8 2 = - 2 . t w z + 1 2, J, k + 1 2, n + 1 - 2w z + 1 2, J, + 1 2, n ~t +w(i+ 1l2,j,k+ 1l2,n-1)) = Bgh [Tzz(i + 112,}, k + 1)- Tzz(i + 1l2,j, k)] 1 - 24h [Tzz(i+ 1l2,j,k+2) -Tzz(i+ 1l2,j,k-1)) + Bgh [rzx(i + 1,j, k + 112)- Tzx(i,j, k + 112)) - 2 !h [Tzx(i + 2,j, k + 112)- Tzx(i- 1,j, k + 112)) + B~l [Tzy(i + 1l2,j + 112, k + 112) -Tzy(i + 1l2,j- 112, k + 112)) - 2!h [rzy(i + 1l2,j + 312, k + 112) -Tzy(i + 1l2,j- 312, k + 112)) (A.3) Txx(i + 1/2, j, k) = {2p +A) ( 8 9 h (u{i +I, j, k)- u(i, j, k)] - 2~h [u(i + 2,j, k)- u(i- I,j, k)J) +A ( S~t [v(i + 1/2, j+ 1/2, k)- v(i + 1/2, j- 1/2, k)] 1 - 24h [v(i + 1l2,j + 312, k)- v(i + 1l2,j- 312, k)] 9 + Bh [w(i + 112, j, k + 112)- w(i + 112, j, k- 112)) - 2!h [w(i + l/2,j, k + 3/2)- w(i + I/2,j, k- 3/2)]) (A.4)