I am trying to model the seismic wave equation and have therefore been reading about discretization schemes and their stability. I recently came across an insightful paper on 'Galerkin FEM methods for wave equation' http://www.math.tamu.edu/~bangerth/publications/2010-waves.pdf. In this document and many other documents, for the newmark beta scheme, the final discretised form of wave equation is in terms of acceleration, velocity and displacement. That is how I have seen it in other documents as well.

enter image description here

However, I have also seen another form of Newmark scheme of wave equation, i.e., the following: (Pg 130)

enter image description here


Now, I am really confused as to how a wave equation can have 2 different forms of Newmark discretisation. I' ll be very grateful if you could shed some light on this or point me in the right direction.

The end purpose is to model it in fenics and DEAL 2

  • $\begingroup$ The $\beta \in [0, 0.5]$ and $\gamma \in [0, 1]$ values in Newmark's method interpolate linearly between $t_n$ and $t_{n+1}$ quantities. If you have a $t_{n-1}$ quantity in your algorithm it's not a Newmark-$\beta$ scheme. $\endgroup$ Nov 21, 2015 at 1:32


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.