Changes between Version 4 and Version 5 of UniformCollapse


Ignore:
Timestamp:
08/09/11 14:14:55 (13 years ago)
Author:
Jonathan
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • UniformCollapse

    v4 v5  
    1212[attachment:rho.gif movie]
    1313
    14 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...
     14Artificial 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...
    1515
    1616
     17== Resolution study ==
     18The 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. 
     19
     20[attachment:sink_comparison.gif movie]
     21
     22[[Image(sink_comparison0030.png, width=50%)]]
     23
     24== Choice of !PhiDot ==
     25
     26To investigate the wandering we modified the program to allow for calculating [[latex($\frac{\partial \phi}{\partial t}$)]] 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 the [attachment:nonthreaded.gif movie] shows it seems to have fixed the wondering problem.
     27
     28[[Image(nonthreaded0030.png, width=50%)]]