“…Consider the third order in time nonlinear partial differential equation model in abstract form ∂τττ y(x, τ ) + α∂ττ y(x, τ ) − c 2 ∂xxy(x, τ ) − b∂xxτ y(x, τ ) = G(x, τ , y, yτ , yττ ) (1) where G(x, t, y, yτ , yττ ) is a non-linear or linear function and α, c, b > 0 are given parameters. This type of model is often called Moore-Gibson-Thompson equation and appeared in many scientific fields such as nonlinear acoustics, medical ultrasound, viscoelasticity and thermoelasticity, [4,6,[10][11][12]20].…”