“…Here the fractional differential operator ∂ γ t := 0 I 1−γ t ∂ t for 0 ≤ γ < 1 with the fractional integral operator 0 I 1−γ t g := (t −γ /Γ(1 − γ)) * g and ∂ γ t := ∂ t u for γ = 1 [30]. Hence, the viscoelastic damping term −κ ∫ ∆Ω ∂ γ t udx accurately describes the behavior of viscoelastic damping, and includes the elastic resistance and viscoelastic damping as special cases [3,4,8,25,26,30,37,38], and has attracted growing research activities [9,11,13,17,18,21,24,28,31,32,41,42,44]. Consequently, the modeling equation becomes ∂ 2 t u(x, t) + κ∂ γ t u − K∇ 2 u(x, t) = f (x, t), (x, t) ∈ Ω × (0, T ], u(x, t) = 0, (x, t) ∈ ∂Ω × [0, T ]; u(x, 0) = u 0 (x), ∂ t u(x, 0) =ǔ 0 (x), x ∈ Ω,…”