Abstract
Using numerical methods, we examine, under a Gowdy symmetry assumption, the dynamics of nonlinearly perturbed Friedmann-Lemaître-Robertson-Walker (FLRW) fluid solutions of the Einstein-Euler-scalar field equations in the contracting direction for linear equations of state 𝑝=𝐾𝜌 and sound speeds 0≤𝐾<1/3. This article builds upon the numerical work from Beyer et al. [Phys. Rev. D 107, 104030 (2023)] in which perturbations of FLRW solutions to the Einstein-Euler equations with positive cosmological constant in the expanding time direction were studied. The numerical results presented here confirm that the instabilities observed in previous work [Phys. Rev. D 107, 104030 (2023), Lett. Math. Phys. 113, 102 (2023).] for 1/3<𝐾<1, first conjectured to occur in the expanding direction by Rendall [Ann. Henri Poincaré 5, 1041 (2004)], are also present in the contracting direction over the complementary parameter range 0≤𝐾<1/3. Our numerical solutions show that the fractional density gradient of the nonlinear perturbations develop steep gradients near a finite number of spatial points and become unbounded toward the big bang. This behavior, and in particular the characteristic profile of the fractional density gradient near the big bang, is strikingly similar to what was observed in the expanding direction near timelike infinity by Beyer et al. [Phys. Rev. D 107, 104030 (2023)].