The effective field theory for hydrodynamics allows to write the action functional for fluid. In this paper, some simplest possible higher derivative terms in the fluid action and the cosmological consequences of their presence are considered. Particular attention is given to dark matter, modelled as a dust with higher derivative corrections. We study the conditions of absence of singularities in the solutions of the background and perturbed equations and investigate the evolution of perturbations in two simple models of matter dominated Universe. There is a range of parameters describing the higher derivative terms, in which the short-wavelength perturbations of dark matter are suppressed and the dark matter can be seen as fairly homogeneous on a sufficiently small scale.

The spatial averaging used for the splitting of the local scale factor on the homogeneous background and small inhomogeneous perturbation leads to a non-local relationship between locally and globally defined comoving curvature perturbations. We study this relationship within a quasi-homogeneous, nearly spatially flat domain of the Universe. It is shown that, on scales larger than the size of the observed patch, the Fourier components of the locally defined comoving curvature perturbation are suppressed. We have also shown that the statistical properties of local and global comoving curvature perturbations are coincide on a small scale. Several examples are discussed in detail.

The general expression for the second order large scale curvature perturbation in the form of a functional over a background solution is derived. The explicit expressions was obtained for two special forms of the inflationary potential. In the considered cases, it is shown that a significant level of non-Gaussianity can be generated during the super-Hubble evolution only if nonadiabatic perturbations are non-negligible at the end of inflation.

The evolution of non-adiabatic perturbations in models with multiple coupled perfect fluids with non-adiabatic sound speed is considered. Instead of splitting the entropy perturbation into relative and intrinsic parts, we introduce a set of symmetric quantities, which also govern the non-adiabatic pressure perturbation in models with energy transfer. We write the gauge invariant equations for the variables that determine on a large scale the non-adiabatic pressure perturbation and the rate of changes of the comoving curvature perturbation. The analysis of evolution of the non-adiabatic pressure perturbation has been made for several particular models.

The scalar perturbations in inflationary models, based on a two-component diagonal non-linear sigma model, are considered. For inhomogeneities generated at an inflationary stage, the law of motion of the comoving curvature ${\cal R}$ is obtained (without using the slow roll approximations). A formal expressions for the power spectrum and its spectral index are obtained, which are valid and after the slow-roll stage. As an example, an inflationary model with a massive scalar field and quadratic curvature corrections is studied numerically.