The intention of this paper is mainly two-fold. \textit{First}, we point out a striking numerical agreement between the bulk viscosity in the lepton era calculated by Husdal (2016) and our own calculations of the present-day bulk viscosity when the functional form is $\zeta\,\sim\sqrt{\rho}$. From a phenomenological point of view, we thus seem to have an ansatz for the viscosity which bridges the infancy of the Universe ($\sim 1$ s) with the present. This can also be looked upon as a kind of symmetry between the early-time cosmology and the present-day cosmology: it is quite remarkable that the kinetic theory-based bulk viscosity in the early universe and the experimentally-based bulk viscosity in the present universe can be covered by the same simple analytical formula. \textit{Second}, we consider the Kasner universe as a typical anisotropic model of Bianchi-type I, investigating whether this geometrical model is compatible with constant viscosity coefficients in the fluid. Perhaps surprisingly, the existence of a shear viscosity turns out to be incompatible with the Kasner model. By contrast, a bulk viscosity is non-problematic in the {\it isotropic} version of the model. In the special case of a Zel'dovich (stiff) fluid, the three equal exponents in the Kasner metric are even determined by the bulk viscosity alone, independent of the value of the fluid energy density. We also give a brief comparison with some other recent approaches to viscous cosmology.

We analyze gravitational waves propagating in an isotropic cosmic fluid endowed with a bulk viscosity $\zeta$ and a shear viscosity $\eta$, assuming these coefficients to vary with fluid density $\rho$ as $\rho^\lambda$, with $\lambda=1/2$ favored by experimental evidence. We give the general governing equation for the gravitational waves, and focus thereafter on two examples. The first concerns waves in the very late universe, close to the Big Rip, where the fate of the comic fluid is dependent highly on the values of the parameters. Our second example considers the very early universe, the lepton era; the motivation for this choice being that the microscopical bulk viscosity as calculated from statistical mechanics is then at maximum. We find that the gravitational waves on such an underlying medium are damped, having a decay constant equal to the inverse of the conformal Hubble parameter. Our results turn out to be in good agreement with other viscosity-based approaches.

Gas hydrates formed in oceans and permafrost occur in vast quantities on Earth representing both a massive potential fuel source and a large threat in climate forecasts. They have been predicted to be important on other bodies in our solar systems such as Enceladus, a moon of Saturn. CO$_2$-hydrates likely drive the massive gas-rich water plumes seen and sampled by the spacecraft Cassini, and the source of these hydrates is thought to be due to buoyant gas hydrate particles. Dispersion forces cause gas hydrates to be coated in a 3-4 nm thick film of ice, or to contact water directly, depending on which gas they contain. These films are shown to significantly alter the properties of the gas hydrate clusters, for example, whether they float or sink. It is also expected to influence gas hydrate growth and gas leakage.

We look at the entropy production during the lepton era in the early Universe by using a model where we exclude all particles except the leptons and photons. We assume a temperature dependent viscosity as calculated recently by one of us (Husdal 2016) with use of relativistic kinetic theory. We consider only the bulk viscosity, the shear viscosity being omitted because of spatial isotropy. The rate of entropy production is highest just before the neutrinos decouple. Our results show that the increase in entropy during the lepton era is quite small, about 0.071 % at a decoupling temperature of $T=10^{10}~\mathrm{K}$. This result is slightly smaller than that obtained earlier by Caderni and Fabbri (1977). After the neutrino decoupling, when the Universe has entered the photon era, kinetic-theory arguments no longer support the appearance of a bulk viscosity. At high temperatures and a stable particle ratio, entropy production (d$\sigma$/d$T$) goes as $T^{-8}$, with the total entropy ($\Delta \sigma$) increasing as $T^{-7}$. These rates go slightly down just before the neutrinos decouple, where $\Delta \sigma \propto T^{-6.2}$.

We investigate the early-time accelerated universe after the Big Bang. We pay attention to the dissipative properties of the inflationary universe in the presence of a soft type singularity, making use of the parameters of the generalized equation of state of the fluid. Flat Friedmann-Robertson-Walker metric is being used. We consider cosmological models leading to the so-called type IV singular inflation. Our obtained theoretical results are compared with observational data from the Planck satellite. The theoretical predictions for the spectral index turn out to be in agreement with the data, while for the scalar-to tensor ratio there are minor deviations.

We analyze characteristic properties of two different cosmological models: (i) a one-component dark energy model where the bulk viscosity $\zeta$ is associated with the fluid as a whole, and (ii) a two-component model where $\zeta$ is associated with a dark matter component $\rho_{\rm m}$ only, the dark energy component considered inviscid. Shear viscosity is omitted. We assume throughout the simple equation of state $p=w\rho$, with $w$ a constant. In the one-component model we consider two possibilities, either to take $\zeta$ proportional to the scalar expansion (equivalent to the Hubble parameter), in which case the evolution becomes critically dependent on the value of the small constant $\alpha=1+w$ and the magnitude of $\zeta$. Second, we consider the case $\zeta=~$const., where a de Sitter final stage is reached in the future. In the two-component model we consider only the case where the dark matter viscosity $\zeta_{\rm m}$ is proportional to the square of $\rho_{\rm m}$, where again a de Sitter form is found in the future. In this latter case the formalism is supplemented by a phase space analysis. As a general result of our considerations we suggest that a value $\zeta_0\sim 10^6~$Pa s for the present viscosity is reasonable, and that the two-component model seems to be favored.

We investigate cosmological models with a linear inhomogeneous time-dependent equation of state for the dark energy, coupled with dark matter, leading to a bounce cosmology. Equivalent descriptions in terms of the equation-of-state parameters for an exponential, a power-law, or a double exponential form of the scale factor $a$ is obtained. The stability of the solutions is explored, by considering small perturbations around the critical points for the bounce in the early and in the late-time universe.

Dark energy of phantom or quintessence nature with an equation of state parameter $w$ almost equal to -1 often leads the universe evolution to a finite-time future singularity. An elegant solution to this problem has been recently proposed \cite{frampton11} under the form of the so-called Little Rip cosmology which appears to be a realistic alternative to the $\Lambda$CDM model. A viscous Little Rip cosmology is here proposed. Whereas generically bulk viscosity tends to promote the Big Rip, we find that there are a number of situations where this is not the case and where the formalism nicely adjusts itself to the Little Rip scenario. We prove, in particular, that a viscous fluid (or, equivalently, one with an inhomogeneous (imperfect) equation of state) is perfectly able to produce a Little Rip cosmology as a purely viscosity effect. The possibility of its induction as a combined result of viscosity and a general (power-like) equation of state is also investigated in detail. To finish, a physical, inertial force interpretation of the dissolution of bound structures in the Little Rip cosmology is presented.

As first part of this work, experimental information about the decay of isotropic turbulence in ordinary hydrodynamics, u^2(t) proportional to t^{-6/5}, is used as input in FRW equations in order to investigate how an initial fraction f of turbulent kinetic energy in the cosmic fluid influences the cosmological development in the late, quintessence/phantom, universe. First order perturbative theory to the first order in f is employed. It turns out that both in the Hubble factor, and in the energy density, the influence from the turbulence fades away at late times. The divergences in these quantities near the Big Rip behave essentially as in a non-turbulent fluid. However, for the scale factor, the turbulence modification turns out to diverge logarithmically. As second part of our work, we consider the full FRW equation in which the turbulent part of the dark energy is accounted for by a separate term. It is demonstrated that turbulence occurrence may change the future universe evolution due to dissipation of dark energy. For instance, phantom-dominated universe becomes asymptotically a de Sitter one in the future, thus avoiding the Big Rip singularity.

We derive a Cardy-Verlinde-like formula which relates the entropy of the closed FRW universe to its energy, and Casimir energy, for a multicomponent coupled fluid. The generalized fluid obeys an inhomogeneous equation of state. A viscous dark fluid is included, and we include also modified gravity using its fluid representation. It is demonstrated how such an expression reduces to the standard Cardy-Verlinde formula corresponding to the 2d CFT entropy only in some special cases. The dynamical entropy bound for a closed FRW universe with dark components is obtained. The universality of the dynamical entropy bound near a future singularity (of all four types), as well as near the Big Bang singularity, is investigated. It is demonstrated that except from some special cases of Type II and Type IV singularities the dynamical entropy bound is violated near the singularity even if quantum effects are taken into account. The dynamical entropy bound seems to be universal for the case of a regular universe, including the asymptotic de Sitter one.

Singularities in the dark energy late universe are discussed, under the assumption that the Lagrangian contains the Einstein term R plus a modified gravity term of the form R^\alpha, where \alpha is a constant. It is found, similarly as in the case of pure Einstein gravity [I. Brevik and O. Gorbunova, Gen. Rel. Grav. 37 (2005), 2039], that the fluid can pass from the quintessence region (w>-1) into the phantom region (w<-1) as a consequence of a bulk viscosity varying with time. It becomes necessary now, however, to allow for a two-fluid model, since the viscosities for the two components vary differently with time. No scalar fields are needed for the description of the passage through the phantom barrier.

We study the entropy of a FRW universe filled with dark energy (cosmological constant, quintessence or phantom). For general or time-dependent equation of state $p=w\rho$ the entropy is expressed in terms of energy, Casimir energy, and $w$. The correspondent expression reminds one about 2d CFT entropy only for conformal matter. At the same time, the cosmological Cardy-Verlinde formula relating three typical FRW universe entropies remains to be universal for any type of matter. The same conclusions hold in modified gravity which represents gravitational alternative for dark energy and which contains terms growing at low curvature. It is interesting that BHs in modified gravity are more entropic than in Einstein gravity. Finally, some hydrodynamical examples testing new shear viscosity bound, which is expected to be the consequence of the holographic entropy bound, are presented for the early universe in the plasma era and for the Kasner metric. It seems that the Kasner metric provides a counterexample to the new shear viscosity bound.

We consider the gravitationally induced particle production from the quantum vacuum which is defined by a free, massless and minimally coupled scalar field during the formation of a gauge cosmic string. Previous discussions of this topic estimate the power output per unit length along the string to be of the order of $10^{68}$ ergs/sec/cm in the s-channel. We find that this production may be completely suppressed. A similar result is also expected to hold for the number of produced photons.