It is found that our universe expands with acceleration [1-5]. The accelerating expansion of the universe may be explained in context of the dark energy [6]. Due to negative pressure, the simplest way for modeling the dark energy is the Einstein’s cosmological constant. On the other hand, the study of the cosmological constant is one of the important subjects in the theoretical and experimental physics [7-10]. Another candidate for the dark energy is scalarfield dark energy model [11-19]. However, presence of a scalar field is not only requirement of the transition from a universe filled with matter to an exponentially expanding universe. Therefore, Chaplygin gas is used as an exotic type of fluid, which is based on the recent observational fact that the equation of state parameter for dark energy can be less than −1.

On the other hand, we know that the viscosity plays an important role in the cosmology [20]. In another word, the presence of viscosity in the fluid introduces many interesting pictures in the dynamics of homogeneous cosmological models, which is used to study the evolution of universe. In Ref. [21], the exact solutions of the field equations for a five-dimensional space-time with viscous fluid were obtained. Also in Ref. [22] a cosmological model with viscous fluid in higher-dimensional spacetime was constructed. Then, in Ref. [23] the exact solutions of the field equations for a five-dimensional cosmological model with variable bulk viscosity were obtained. The isotropic homogeneous spatially flat cosmological model with bulk viscous fluid was constructed in Ref. [24]. The bulk viscous cosmological models with constant bulk viscosity coefficient were constructed in Ref. [25]. In the recent work [26] the FRW bulk viscous cosmology was considered and bulk viscous coefficient was obtained in the flatspace, and then extended to non-flat space [27]. In this work, we consider both bulk viscous effect and Chaplygin gas in FRW cosmology in flat space.

The Friedmann-Robertson-Walker (FRW) universe in four-dimensional space-time is described by the following metric [28,29],

where, and represents the scale factor. The and are the usual azimuthal and polar angles of spherical coordinates. Also, constant denotes the curvature of the space. In this paper we consider the case of only, which is corresponding to flat space. In that case the Einstein equation is given by,

where we assumed and. Also the energy-momentum tensor corresponding to the bulk viscous fluid and modified Chaplygin gas [30-35] is given by the following relation,

where is the energy density and is the velocity vector with normalization condition. Also, the total pressure and the proper pressure involve bulk viscosity coefficient and Hubble expansion parameter are given by the following equations [36-42],

and,

with and. The equation of state is one of the most important quantity to describe the features of dark energy models. It is clear that the parameter shows bulk viscosity and B shows effect of Chaplygin gas. In the Ref. [43] the dynamics of FRW cosmology with modified Chaplygin gas as the matter formulated. Then the nature of the critical points are studied by evaluating the eigenvalues of the linearized Jacobi matrix for the special case of. In this paper we consider special case with and extend the Ref. [43] to including bulk viscous coefficient.

In that case the Friedmann equations are given by,

and,

where dot denotes derivative with respect to cosmic time. The energy-momentum conservation law obtained as the following,

In the next section we try to obtain time-dependent density by using above equations.

Using the Equations (4)-(6) in the conservation relation (8) we have,

If we set, then one can extract energy density depend on scale factor [43],

where c is an integration constant. Here we also consider bulk viscous coefficient and would like to obtain energy density depend on time. In order to solve Equation (9) we use the following ansatz,

where constants, , , and should be determined. Substituting relation (11) in the Equation (9) gives us the following coefficients,

If we neglect both bulk viscosity and presence of Chaplygin gas then,

which is agree with results of the Refs. [27,43] where established. On the other hand for the large bulk viscosity coefficient one can find that and hence obtained. Also for the case of infinitesimal

where,

one can obtain constant negative energy density. In the general case, Equation (11) with coefficients (12)-(16) tells us that the energy density is decreasing function of time. Such behavior happen for the Hubble expansion parameter which is discussed below.

By using time-dependent density in the relation (6) one can obtain Hubble expansion parameter. In that case we draw plot of Hubble expansion parameter in the Figure 1 for.

In that case the modified Chaplygin gas model describes the evolution of the universe from the radiation regime to the Λ-cold dark matter scenario, where the fluid behaves as a cosmological constant, so there is an accelerated expansion of the universe.

It is possible to study deceleration parameter of this theory which obtained by the following relation,

Numerically we draw deceleration parameter in terms of time in the Figure 2.

In this work, we studied the FRW bulk viscous cosmology with modified Chaplygin gas as the matter contained. We obtained the modified Friedmann equations due to bulk viscous and Chaplygin gas coefficients. Then tried to solve equations and found time-dependent energy density. Therefore, we could extract Hubble expansion and deceleration parameters.

For the future work, it is possible to repeat calculation of this paper for the case of arbitrary α or non-flat universe where. In that case one deals with the following equation,

where.