Co-rotating binaries disk formation sims No .1

6 June '12, 10 AU case, Mdotprimary=10-5Mo yr-1

http://www.pas.rochester.edu/~martinhe/2011/binary/6jun12.png Disk mass evolution. Compare with pk model 2, fig. 1 top, but note that qpk=3 while qus=1.5. I'll soon have a similar plot for the 20 AU case which should be compared with M&M figure 5 bottom; M&M see an increasing disk mass with a final value of 5x10-6 Mo yr-1.
http://www.pas.rochester.edu/~martinhe/2011/binary/7jun12.png Evolution of the mass accretion rate onto the particle. The converging value of ~5x10-6 seems a factor of 50 higher than expected. Comparing with table 3 M&M, models 1 and 2: MdotM&M~.9-3x10-6 ; Mdotus~5x10-6, so we seem to be off for a factor of order 1.6-5.6. Comparing with pk model 2, fig. 1 bottom, (but note that qpk=3 while qus=1.5) we are off by a factor of ~30. Investigating further.

-20au preliminary early evo test with more resolution than before, ​http://www.pas.rochester.edu/~martinhe/2011/binary/20au-1000k-64x64x32-4amr-correctWINDdens.gif. This sim has a resolution of 64x64x32+3amr. AGB wind inflow enters the -x, +y and +-z boundaries. Improved versions of this sim are running in ranger and queued in kraken.

authors q=m1/m2 a [AU] tempw [k] velw [km/s] total run time [yr] resolution
us 1.5/1=1.5 10 1000 10 40 (2 orb) soon
us 1.5/1=1.5 20 1000 10 57 (2 orb) soon
us 1.5/1=1.5 30 1000 10 104 (2 orb) soon
us 1.5/1=1.5 40 1000 10 160 (2 orb) soon
pk 3/1=3 3 soon 10-20 104-6 soon
pk 3/1=3 10 soon 10-20 104-6 soon
pk 3/1=3 30 soon 10-20 104-6 soon
pk 3/1=3 100 soon 10-20 104-6 soon
pk 1.8/.6=3 3 soon 10-20 104-6 soon
pk 1.8/.6=3 10 soon 10-20 104-6 soon
pk 1.8/.6=3 30 soon 10-20 104-6 soon
pk 1.8/.6=3 100 soon 10-20 104-6 soon
vb 1.2/.6=2 10-100
  • velw=10km/s
  • tempw=1000K
  • mprim=1 Mo
  • msec = .5 Mo
  • resolution= 633 cells + 2 particle refinements
  • rBondi= G msec / (velw2 + velorb-second2)
  • lgrid = 2.5rBondi
  • rsoft = 4 cells

Inflow form the top, bottom, lower Y and larger X is killed.

date separation [AU] res rmax-ref-region rBondi [cells] rBondi/rsoft velorbit-sec [km/s] run disk moviesa
14may (7 orbits) 10 643 + 2amr rBondi/4 100 25 10. ~11hr/orbit, 32p, Bhive y ​http://www.pas.rochester.edu/~martinhe/2011/binary/4panels-10may12B.gif
19may 10 643 0 100 25 10. ~5hr/16p, ranger y ​http://www.pas.rochester.edu/~martinhe/2011/binary/4panels-20may12.gif
18may 10 643 + 2amr rBondi/2 100 25 10. ~52.5hr/64p, ranger y ​http://www.pas.rochester.edu/~martinhe/2011/binary/4panels-18may12.gif
14may (6 orbits) 20 643 + 2amr rBondi/4 100 25 6.4 ~17hr/orbit, 40p, Bhive y ​http://www.pas.rochester.edu/~martinhe/2011/binary/4panels-14may12a.gif
14may 30 643 + 2amr rBondi/4 100 25 5.2 ~42.5hr/128p, ranger (long q) small? ​http://www.pas.rochester.edu/~martinhe/2011/binary/4panels-14may12b.gif
18may 30 1283 + 2amr rBondi/2 100 25 5.2 running (1.6 orbits) 64p Bhive … ​http://www.pas.rochester.edu/~martinhe/2011/binary/4panels-21may12.gif
20may (stand by) 30 1283 + 2amr rBondi/4 100 25 5.2 running INTERMITTENTLYc since 2pm, 1536proc, kraken … …
20may (stand by) 30 643 + 3amr rBondi/4 100 25 5.2 running INTERMITTENTLYc since noon, 1536proc, kraken … …
14may 40 643 + 2amr rBondi/4 100 25 4.4 ~48hr/128p, ranger (long q) n ​http://www.pas.rochester.edu/~martinhe/2011/binary/4panels-14may12c.gif
23may (stand by) 40_2D 128x128x2 + 3amr rBondi/8 100 25 4.4 …/32p Bhive (q @ krakrn, Bgene) ? soon

a No. density log grey scale [cu] + color velocity [Mach]. Top left, top right, bottom left and bottom right show the entire orbital plane, a zoom into the orbital plane, the entire longitudinal plane for a phi angle parallel to the wind at the origin, and a zoom into the longitudinal plane, respectively.

c ​https://clover.pas.rochester.edu/trac/astrobear/ticket/209#trac-add-comment

Comments

No comments.