Meeting update

Here is the new Jeans instability calculation of the growth rate

I see that the code reproduces the growth rate well up until ~4 e-folding times. But by 5 e-folding times, the difference is ~10% in the perturbation. This is much better than it's been: ​https://astrobear.pas.rochester.edu/trac/astrobear/blog/erica06262013

I could fidget around with the code a bit more, if necessary - Is this good enough for a test of self-gravity for the paper?

Also, here are fun plots showing the velocity and density perturbations on the same graph

t = 0:

Vx is in red, rho is in green. The density peaks at the zero's of Cos(kx). E.g., a local max is at 1. The vx perturbation is a Sin(x) function, and it has minima that corresponds with the density maxima. This produces converging flows at the max in density, as you would expect you need to develop the Jeans instability.

t = tfinal:

By the end of the sim, the density perturbation has grown a little, but velocity has grown a lot.

Attachments (3)

Download all attachments as: .zip

Comments

No comments.