Uniform Collapse of a Sphere
The free fall time should be 1.5323
A single sink particle appears at approximately that time but then drifts from the center of the grid. This may be due to poor resolution (only based on the Jeans length)
Artificial fragmentation will occur whenever the density violates the TrueLove criterion. However, the Federrath criteria for forming sink particles requires more then just a single cell above the truelove criterion. The cell must be in the minimum of a local potential well, have converging flow from each direction, and the surrounding sphere with a radius of 4*dx where dx is the finest level resolution must also be a bound system. If the sphere of 4*dx is much larger then the collapsing condensation, then it can lead to problems since particles will likely not form although the jeans criterion will lead to artificial fragmentation. The only solution is to add more resolution…
Resolution study
The same problem was run with a maxlevel of 0 through 5. Note that when the MaxLevel is 0, 4 particles are created instantly (although this is puzzling since they are not all local potential minima) Also note that for MaxLevel = 5, the particle wanders.
Choice of PhiDot
To investigate the wandering we modified the program to allow for calculating movie shows it seems to have fixed the wandering problem.
for a coarse step and then prolongating this to the finer level to use for updating the potential in ghost zones. This required advancing the coarse step before the finer steps and disables level threading. However as theAttachments (7)
 rho.gif (1.7 MB )  added by 13 years ago.
 rho0016.jpeg (38.4 KB )  added by 13 years ago.

sink_comparison0030.png
(10.5 KB
)  added by 13 years ago.
Frame of movie showing truelove uniform collapse problem with different Maxlevels (0 to 5)

sink_comparison.gif
(1.1 MB
)  added by 13 years ago.
Movie showing truelove uniform collapse problem with different Maxlevels (0 to 5)

nonthreaded0030.png
(44.6 KB
)  added by 13 years ago.
Movie showing the MaxLevel=5 case run but without using old time derivatives to prolongate/interpolate phi at coarsefine boundaries

nonthreaded.gif
(1.2 MB
)  added by 13 years ago.
Movie showing the MaxLevel=5 case run but without using old time derivatives to prolongate/interpolate phi at coarsefine boundaries

sink_comparison.flv
(363.7 KB
)  added by 13 years ago.
Sink Comparison Movie