We study the effects associated with nonlinearity of $f(R)$ gravity and of the background perfect fluid manifested in the Kaluza-Klein model with spherical compactification. The background space-time is perturbed by a massive gravitating source which is pressureless in the external space but has an arbitrary equation of state (EoS) parameter in the internal space. As characteristics of a nonlinear perfect fluid, the squared speeds of sound are not equal to the background EoS parameters in the external and internal spaces. In this setting, we find exact solutions to the linearized Einstein equations for the perturbed metric coefficients. For nonlinear models with $f^{\prime\prime}(R_0)\neq0$, we show that these coefficients acquire correction terms in the form of two summed Yukawa potentials and that in the degenerated case, the solutions are reduced to a single Yukawa potential with some "corrupted" prefactor (in front of the exponential function), which, in addition to the standard $1/r$ term, contains a contribution independent of the three-dimensional distance $r$. In the linear $f''(R)=0$ model, we generalize the previous studies to the case of an arbitrary nonlinear perfect fluid. We also investigate the particular case of the nonlinear background perfect fluid with zero speed of sound in the external space and demonstrate that a non-trivial solution exists only in the case of $f''(R_0)=0$.

A five-dimensional cosmological model including a single perfect fluid is studied in the framework of dynamical system analysis. All the critical points of the system with their stability properties are listed and some representative phase diagrams are explicitly shown. It is found that the stabilization of the extra dimension is possible and the observed flatness of the three-dimensional space is provided for certain ranges of the equation of state parameter of the fluid during the evolution of the universe. The model suggested here can be considered as a simplified model for examining the possible effects of the extra dimensions in the early universe.

Two distinct perturbative approaches have been recently formulated within General Relativity, arguing for the screening of gravity in the $\Lambda$CDM Universe. In this paper we compare them and show that the offered screening concepts, each characterized by its own interaction range, can peacefully coexist. Accordingly, we advance a united scheme, determining the gravitational potential at all scales, including regions of nonlinear density contrasts, by means of a simple Helmholtz equation with the effective cosmological screening length. In addition, we claim that cosmic structures may not grow at distances above this Yukawa range and confront its current value with dimensions of the largest known objects in the Universe.