We investigate cosmological perturbations for nonlinear $f(R)$ models within the cosmic screening approach. Matter is considered both in the form of a set of discrete point-like massive bodies and in the form of a continuous pressureless perfect fluid. We perform full relativistic analysis of the first-order theory of scalar perturbations for arbitrary nonlinear $f(R)$ models and demonstrate that scalar potentials $\Phi(t,\mathbf{r})$ and $\Psi(t,\mathbf{r})$ are determined by a system of only two master equations. Our equations are applicable at all spatial scales as long as the approximation $\delta R/\bar R \ll 1$ (which is usually assumed in studies devoted to cosmological perturbations in $f(R)$ models) works.

We investigate dark and visible matter distribution in the Coma cluster in the case of the Navarro-Frenk-White (NFW) profile. A toy model where all galaxies in the cluster are concentrated inside a sphere of an effective radius $R_{\mathrm{eff}}$ is considered. It enables us to obtain the mean velocity dispersion as a function of $R_{\mathrm{eff}}$. We show that, within the observation accuracy of the NFW parameters, the calculated value of $R_{\mathrm{eff}}$ can be rather close to the observable cutoff of the galaxy distribution. Moreover, the comparison of our toy model with the observable data and simulations leads to the following preferable NFW parameters for the Coma cluster: $R_{200} \approx 1.77\,h^{-1} \, \mathrm{Mpc} = 2.61\, \mathrm{Mpc}$, $c=3\div 4$ and $M_{200}= 1.29h^{-1}\times10^{15}M_{\odot}$. In the Coma cluster the most of galaxies are concentrated inside a sphere of the effective radius $R_{\mathrm{eff}}\sim 3.7$ Mpc and the line-of-sight velocity dispersion is $1004\, \mathrm{km}\, \mathrm{s}^{-1}$.

We make a generalization of a self-consistent first-order perturbation scheme, being suitable for all (sub-horizon and super-horizon) scales, which has been recently constructed for the concordance cosmological model and discrete presentation of matter sources, to the case of extended models with extra perfect fluids and continuous presentation. Namely, we derive a single equation determining the scalar perturbation and covering the whole space as well as define the corresponding universal Yukawa interaction range. We also demonstrate explicitly that the structure growth is suppressed at distances exceeding this fundamental range.

In the $\Lambda$CDM framework, presenting nonrelativistic matter inhomogeneities as discrete massive particles, we develop the second-order cosmological perturbation theory. Our approach relies on the weak gravitational field limit. The derived equations for the second-order scalar, vector and tensor metric corrections are suitable at arbitrary distances including regions with nonlinear contrasts of the matter density. We thoroughly verify fulfilment of all Einstein equations as well as self-consistency of order assignments. In addition, we achieve logical positive results in the Minkowski background limit. Feasible investigations of the cosmological backreaction manifestations by means of relativistic simulations are also outlined.

We study the late-time evolution of the Universe where dark energy (DE) is parametrized by a modified generalized Chaplygin gas (mGCG) on top of cold dark matter (CDM). We also take into account the radiation content of the Universe. In this context, the late stage of the evolution of the universe refers to the epoch where CDM is already clustered into inhomogeneously distributed discrete structures (galaxies, groups and clusters of galaxies). Under these conditions, the mechanical approach is an adequate tool to study the Universe deep inside the cell of uniformity. To be more accurate, we study scalar perturbations of the Friedmann-Lema\^itre-Robertson-Walker metric due to inhomogeneities of CDM as well as fluctuations of radiation and mGCG, the later driving the late-time acceleration of the universe. Our analysis applies as well to the case where mGCG plays the role of DM and DE. We select the sets of parameters of the mGCG that are compatible with the mechanical approach. These sets define prospective mGCG models. By comparing the selected sets of models with some of the latest observational data results, we conclude that the mGCG is in tight agreement with those observations particularly for a mGCG playing the role of DE and DM.

We study the late-time evolution of the Universe where dark energy (DE) is presented by a barotropic fluid on top of cold dark matter (CDM). We also take into account the radiation content of the Universe. Here by the late stage of the evolution we refer to the epoch where CDM is already clustered into inhomogeneously distributed discrete structures (galaxies, groups and clusters of galaxies). Under this condition the mechanical approach is an adequate tool to study the Universe deep inside the cell of uniformity. More precisely, we study scalar perturbations of the FLRW metric due to inhomogeneities of CDM as well as fluctuations of radiation and DE. For an arbitrary equation of state for DE we obtain a system of equations for the scalar perturbations within the mechanical approach. First, in the case of a constant DE equation of state parameter $w$, we demonstrate that our method singles out the cosmological constant as the only viable dark energy candidate. Then, we apply our approach to variable equation of state parameters in the form of three different linear parametrizations of $w$, e.g., the Chevallier-Polarski-Linder perfect fluid model. We conclude that all these models are incompatible with the theory of scalar perturbations in the late Universe.