Asymptotic behavior and finite element error estimates of Kelvin-Voigt viscoelastic fluid flow model

In this article, the convergence of the solution of the Kelvin-Voigt viscoelastic fluid flow model to its steady state solution with exponential rate is established under the uniqueness assumption. Then, a semidiscrete Galerkin method for spatial direction keeping time variable continuous is considered and asymptotic behavior of the semidiscrete solution is derived. Moreover, optimal error estimates are achieved for large time using steady state error estimates. Based on linearized backward Euler method, asymptotic behavior for the fully discrete solution is studied and optimal error estimates are derived for large time. All the results are even valid for κ→0, that is, when the Kelvin-Voigt model converges to the Navier-Stokes system. Finally, some numerical experiments are conducted to confirm our theoretical findings.