Meeting update

Here are the Bannerjee plots for my sub-critical sphere (xi=4) that collapsed when the ambient medium was in pressure equilibrium with the edge of the sphere and at the same density. Entire box is isothermal.

Here is density: ​http://www.pas.rochester.edu/~erica/density.png

Here is radial velocity: ​http://www.pas.rochester.edu/~erica/radialVel.png

The scaled times from the Bannerjee paper run from ~3.7-7, so my time samples are slightly different. The radial velocity trend looks nice. As for the density, it seems that the ambient medium falls inward and piles up at the edge of the sphere at earlier times. By the last time curve, the curve has the right shape as the curves in Bannerjee's paper. One big difference I notice is that with my run the ambient seems to have set up some sort of equilibrium with the sphere. In contrast, in Bannerjee's sim the ambient stays uniformly flat throughout the collapse. I see that in that run, the ambient is 100 times lighter than that sphere's outer edge and the box's boundaries are very far from the sphere. There are other differences between the simulations as well, such as there the sphere was a critical sphere with a 10% density enhancement to initialize the collapse. Also, in my simul, a sink particle formed between t/to = 3.84 and 5.12. Is it no longer correct to compare sims after sink formation?

What next steps should be taken?

Comments

No comments.