| 13 | Checking the sink file data, I see that the first sink formed at r~(0.004, 0.004, 0.004). The smallest dx=0.0075, so this is approximately at the center of the center smallest cube on the grid. Ideally the sink would form at the origin, corresponding with the peak in phi. The sink then begins to migrate to the center... but then off the domain. I checked this in the sink data file, and indeed the sink acquired a position with negative coords. At frame 103 or so, a second sink was at position r=(0.001) and the first was at r=(-0.02, -0.02, -0.02). |

| 15 | The initial sink being attracted to the origin makes sense. It does not feel gravitational attraction from material at r>rsink by Gauss's law, and so is attracted to the center of the sphere as it is the location of the deepest part of the potential well. However, how is the sink position mapped to the domain, such that it can be off the grid? Does the sink always form at the center of a cell, and then can travel around subsequently non-restricted to cell centers? If so, why can't it form at any given r, not just (n-1/2)dx? And lastly, if we want to calculate the change in rho(R) over time after sink formation, what is best way to model the problem? Fixed grid, entire box, odd number of cells would insure a sink to form in the center (and remain in center), but the effective resolution would be extremely computationally heavy. Should this be a 1D calculation? |