Stable Keplerian disks TABLE

Cylindrical disks (no grav support height), gamma=1.001

Time increases downwards and the images (click to enlarge) are logarithmic edge-on views of density. Note that to have pressure equilibrium I've set tempdisk=tempamb/100, whereas in the binary simulation these temperatures are the same. Also, in order to show as as much info about the structure of disks as possible, the color scale limits are not the same for all maps.

TABLE 1

rsoft=rd/2; rth=2rd rsoft=rd/2; rth=8rd rsoft=rd/2; rth=16rd rsoft=rd/4; rth=8rd rsoft=rd/8; rth=2rd rsoft=rd/8; rth=8rd rsoft=rd/16; rth=16rd
same but 2amr →same but 2 amr →same but 2 amr →same but 3 amr →same →http://www.pas.rochester.edu/~martinhe/2011/disk/ini-4amr.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-7.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm2-17.pnghttp://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-20.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm2-28.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft2-therm2-100.pnghttp://www.pas.rochester.edu/~martinhe/2011/disk/soft2-therm8-96.pnghttp://www.pas.rochester.edu/~martinhe/2011/disk/soft2-therm16-100.pnghttp://www.pas.rochester.edu/~martinhe/2011/disk/soft4-therm8-100.pnghttp://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm2-100.pnghttp://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-100.png

Flared disks (support height), gamma=1.001

TABLE 2

rsoft=rd/2; rth=2rd rsoft=rd/2; rth=8rd rsoft=rd/2; rth=16rd rsoft=rd/4; rth=8rd rsoft=rd/8; rth=2rd rsoft=rd/8; rth=8rd rsoft=rd/16; rth=16rd
http://www.pas.rochester.edu/~martinhe/2011/disk/soft4-therm8-f0.pnghttp://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm2-f0.pnghttp://www.pas.rochester.edu/~martinhe/2011/disk/ini-4amrf.pnghttp://www.pas.rochester.edu/~martinhe/2011/disk/ini-5amrf.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft16-therm16-1.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft16-therm16-f5.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft4-therm8-f20.pnghttp://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm2-f20.pnghttp://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f20.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft4-therm8-f40.pnghttp://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm2-f40.pnghttp://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f40.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft4-therm8-f60.pnghttp://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm2-f60.pnghttp://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f60.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft4-therm8-f80.pnghttp://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm2-f80.pnghttp://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f80.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft4-therm8-f100.pnghttp://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm2-f100.pnghttp://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f99.png

Flared disks (support height), gamma=1.001, tempd=tempamb (instead of pressd=pressamb):

TABLE 3

rsoft=rd/2; rth=2rd rsoft=rd/2; rth=8rd rsoft=rd/2; rth=16rd rsoft=rd/4; rth=8rd rsoft=rd/8; rth=2rd rsoft=rd/8; rth=8rd rsoft=rd/16; rth=16rd
http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f0-iso.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f20-iso.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f40-iso.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f60-iso.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f80-iso.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f100-iso.png

Same as table 3 (rsoft=rd/8; rth=8rd) but for a Plummer gravity soft profile, instead of a spline one:

TABLE 4

http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f0-iso.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f20-isoP.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f40-isoP.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f60-isoP.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f80-isoP.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f100-isoP.png

Here's a movie (clock to enlarge) of the run in table 4 but with the same color bar than the corresponding column in table 3.

This is zoomed out. Note that BC are periodic:

http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f-iso.gif


Same as table 3 (rsoft=rd/8; rth=8rd) but for gamma=1.667, instead of 1.001:

TABLE 5

http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f0-iso.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f20-isoG.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f40-isoG.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f60-isoG.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f80-isoG.png
http://www.pas.rochester.edu/~martinhe/2011/disk/soft8-therm8-f100-isoG.png

Comments

1. Jonathan -- 13 years ago

Shoudn't the disks be more flared in the initial conditions for various values of rth/r_d ?

2. martinhe -- 13 years ago

Yes Jonathan, the smaller the value of rd, the more flared the disks are.