Universe density in Hyperbolic Universe

[Quarlep]

I want to know Universe density according to this equation( $k=-1$) ?

$H^2(t)-8πρG/3=-k/a^2(t)$
$ρ_U=ρ_m+p_r$
$ρ_U$=Universe density
$ρ_m$=Matter density
$p_r$=Radiation density

 

[wabbit]

With three terms ($1/a^2, 1/a^3, 1/a^4$) the differential equation for $a(t)$ is not solvable analytically, you only have an implicit form in terms of an integral which cannot be expressed in terms of usual functions - you can obtain this form by substituting $\frac{\dot a}{a}$ for $H$, you get something of the form $\int F(a)da=t-t_0$

 

[Quarlep]

Then What should I do ? Is there's any way to escape this situation ?

In $k=-1$ → $Ω_0<1$ Is this can help us ?

 

[wabbit]

It depends what you're looking for. If you're brave enough I think you can express the integral in terms of standard elliptic integrals, which are available in most mathematical packages. But the easy route is just to integrate numerically, this works fine.

$\Omega_0<1$ doesn't help - but the issue is not a hard one, the only thing is, the solution does not have a nice form. It is perfectly computable numerically though.

 

[Quarlep]

If you're brave enough you can express tge integral in terms of standard elliptic integrals, wgich are available in most mathematical packaged. But the easy route is just to integrate numerically, tgis wirks fine.

I am in high school so I did not understand a single word what you're saying

 

[Quarlep]

$Ω_0<1$ doesn't help - but the issue is not a hard one, the only thing is, the solution does not have a nice form. It is perfectly computable numerically though.

I can't do that myself Should I open a new thread or somebody can do that ?

 

[wabbit]

Sorry about that. Very impressive, when I was in high school in didn't know anything at all about the Friedman equation or about General Relativity !

To find $a(t)$ (and therefore $\rho(t)$ ), you need to program two things :
1. The numerical calculation of the integral $G(a)=\int F(a)da$
2. The numerical solution of the equation $G(a)=t$

This is a rather tough with high school tools though, you definitely need calculus. And the answer is a computer program, not a formula.

What are you trying to achieve exactly, is your question about the qualitative behaviour of the solution ? There are some things you can tell without going through a complete solution.

 

[Chalnoth]

I want to know Universe density according to this equation( $k=-1$) ?

$H^2(t)-8πρG/3=-k/a^2(t)$
$ρ_U=ρ_m+p_r$
$ρ_U$=Universe density
$ρ_m$=Matter density
$p_r$=Radiation density

Too many parameters, too few equations. Also, it depends a bit upon your conventions. There are two commonly-used conventions:
1. $k = {-1, 0, 1}$. With this convention, the scale parameter $a$ becomes the radius of curvature.
2. $a(now) = 1$. This is usually the easiest. With this convention, $k$ can take on any number, and represents the amount of spatial curvature today (when $a = 1$).

One way to make things easier to deal with is to use the concept of density fractions. With density fractions, I can rewrite:

H^2(t) = {8\pi G \over 3} \left(\rho_m + \rho_r\right) - {k \over a^2}

as:

H^2(t) = const \left({\Omega_m \over a^3} + {\Omega_r \over a^4} + {\Omega_k \over a^2}\right)

Here I've introduced a constant on the left, and a series of $\Omega$ constants. The equation becomes simplest if we require than when $a = 1$, $\Omega_m + \Omega_r + \Omega_k = 1$. In this case, when $a = 1$, we get:

H^2(now) = const

Where this is the same constant that goes out in front, so our constant is just $H_0^2$, and the Friedmann equation becomes:

H^2(t) = H_0^2 \left({\Omega_m \over a^3} + {\Omega_r \over a^4} + {\Omega_k \over a^2}\right)

So now we have four unknowns, but only one equation. That's not enough to give an answer unless three of the unknowns are set. What's done in practice is to take some measurable quantity which is a function of $H(t)$:

F(data) = f(H(t))

We then take lots of data points, and use computer simulations to figure out which choices for the constants $H_0, \Omega_m, \Omega_r, \Omega_k,$ and $\Omega_\Lambda$ fit the data most closely (that last one stems from the cosmological constant).

 

[Chalnoth]

I can't do that myself Should I open a new thread or somebody can do that ?

I don't think this is exactly what you've asked for, but this website calculates a bunch of different quantities based upon a variety of different inputs:
http://www.astro.ucla.edu/~wright/CosmoCalc.html

 

[Quarlep]

$H_0, \Omega_m, \Omega_r, \Omega_k,$ and $\Omega_\Lambda$ fit the data most closely (that last one stems from the cosmological constant).

In our universe $H_0$ is known 70 and $Ω_m=0.05$ ... so we can change this things isn't it.I mean this quantites describes known universe.So Can I put $H_0=65$ to fit the data ?

 

[Chalnoth]

In our universe $H_0$ is known 70 and $Ω_m=0.05$ ... so we can change this things isn't it.I mean this quantites describes known universe.So Can I put $H_0=65$ to fit the data ?

$\Omega_m$ is not 0.05. This parameter includes both normal matter and dark matter.

Planck's 2015 results report the best-fit estimates of these parameters are:
$H_0 = 68$km/sec/Mpc
$\Omega_m = 0.31$

Anyway, if you tried to fit the data with these parameter values, but $\Omega_\Lambda = 0$, you'd find pretty quickly that the data doesn't fit at all.

 

[Quarlep]

I am trying to find density of universe ) If I put matter density 0.31 then what's the meaning to find density.We have already decided it.Is there any calculator which modified by 2015 plank results.

 

[ChrisVer]

You can't just compute it... you have to make experiments in order to decide its value and try to fit the data... It's a free parameter.

 

[Quarlep]

http://www.astro.ucla.edu/~wright/CosmoCalc.html

In this site How can I put data ?

 

[Quarlep]

Ok,I understand the idea.But I want to ask something.My question is simple so I thought that there's no need to open a new thread and its related this question.
$H^2(t)-8ρπG/3=-k/a^2(t)$ (Lets suppose a(t)=1) then
$H^2(t)-8ρπG/3=-k$ In equations $k$ can be -1,0,1.
Then If $k<0$ then $-k>0$ then
$H^2(t)-8ρπG/3>0$→Hyperbolic Universe. This means $Ω_k>0$ I mean If $k$ negative $Ω_k$ must be positive isn't it ? I am confused here.
In cosmology calculator says $Ω_k=1-Ω_m-Ω_Λ$

(Its look like hyperbolic universe I am trying to understand)

Here I've introduced a constant on the left, and a series of $\Omega$ constants. The equation becomes simplest if we require than when $a = 1$, $\Omega_m + \Omega_r + \Omega_k = 1$.

Thanks