# Problem for simulation on recollapse model

1. Jan 4, 2010

### fab13

Hello,

I post another topic because my first main problem about the modelisation of scale factor by solving numerically friedmann equations is resolved.
However, It remains one problem, that of the recollapse when integration starts at a early age of the universe.

i have reformulated the differential system through the equation :

$$H^2(a) = \bigg(\frac{1}{a} \frac{da}{dt}\bigg)^{2} = H_0^2 \left( \frac{\Omega_m}{a^3} + \frac{\Omega_k}{a^2} + \Omega_\Lambda \right)$$

which gives equation n°1 :

$$\left( \frac{da}{dt} \right)^2 = H_0^2 \left( \frac{\Omega_m}{a} + \Omega_k + \Omega_\Lambda a^{2} \right)$$

then we have :

$$a = \Omega_{m} \bigg( \frac{1}{H_{0}^{2}} \bigg(\frac{da}{dt}\bigg)^{2} - \Omega_{k} - \Omega_{\Lambda} a^{2} \bigg)^{-1}$$ ​

So, the differential equation can be written in this way :

$$a''=\bigg(\frac{-4 \pi G}{3c^{2}}(\frac{\rho_{0} c^{2}}{a^{3}}+3 p)+\frac{\Lambda}{3}\bigg) \Omega_{m} \bigg( \frac{1}{H_{0}^{2}} \bigg(\frac{da}{dt}\bigg)^{2} - \Omega_{k} - \Omega_{\Lambda} a^{2} \bigg)^{-1}$$ ​

Finally, i managed to get the curve representing the "recollapse" with $$\Omega_{m} > 1$$, $$a_{0} = 1$$, and $$a'_{0} = H_{0}$$ and $$\Omega_{\Lambda}=0$$ ( https://www.physicsforums.com/attachment.php?attachmentid=22861&d=1262647111") . The integration begins at $$t_{0}=13.7 Gyr$$. The results are validated by the value of the maximum of scale factor before the recollapse :

$$a_{max} = a ( - \frac{\Omega_{m}}{\Omega_{k}} )$$

$$t_{max} = - \frac{5}{-4} = 1.25$$ on the figure 5 ​

I would like to get the first part of this curve ( for $$t0 << 13.7 Gyr$$ ), ie from different initial conditions. I tried with these conditions :

$$a_{0} = 0.001$$, and $$a'_{0} =H_{0} (5*10^{3}-4)^{1/2} = 70.6824 H_{0}$$

respecting so equation n°1 for the value of $$a'_{0}$$. For this case, the Matlab solver begins integration at $$t_{0}=380.000 years$$, so a redshift z=1000.

I get the https://www.physicsforums.com/attachment.php?attachmentid=22862&d=1262647111". There's well the recollapse design but the maximum is not equal to 1.25 (much more little, about 0.09 ) and the big crunch occurs fastier than on figure 5.

What's the problem, How it is done ? Why i don't get the same recollapse as in figure 5, i mean the same maximum and time "t_final" of big crunch ( $$a(t_{final})=0$$.

Thanks a lot.

Last edited by a moderator: Apr 24, 2017