| | 3 | Fixed the AMR criteria so that just the inner regions of the sphere are refined and moved to modeling the collapse in 3D, beginning with just an octant. Running the code for a freefall time using the density of the inner region (~100 milliseconds) required extremely small initial time steps to avoid high cfls and nans. Kept increasing the initial max speed to get around this. Once the simulation got past a few frames, noticed strange behavior along the boundary, which looked like a pile up of mass there. The boundary conditions are currently set to outflow only on the faces of the cube, and the Poisson BC's to multipole... |
