This paper presents a numerical study of the sensitivity of a fluid model known as time relaxation model with respect to variations of the time relaxation coefficient χ. The sensitivity analysis of this model is utilized by the sensitivity equation method and uses the finite element method along with Crank Nicolson method in the fully discretization of the partial differential equations. We present a test case in support of the sensitivity convergence and also provide a numerical comparison between two different strategies of computing the sensitivity, sensitivity equation method and forward finite differences.
