Changes between Version 4 and Version 5 of u/erica/scratch


Ignore:
Timestamp:
01/27/16 16:05:54 (9 years ago)
Author:
Erica Kaminski
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • u/erica/scratch

    v4 v5  
     1From these plots (described in blog post),
     2
    13[[Image(accretionluminosity1.png, 50%)]]
    24[[Image(accretionluminosity2.png, 50%)]]
    35
     6I can read off the percent error for using various equations to calculate the accretion energy. To get a sense of the error, lets consider a few different initial kinetic energies, ke(r), -- some less than and some greater than the kinetic energy which would have been acquired from freefall alone (i.e. GM/r).
    47
    5 Having a kinetic energy (per unit mass) at r that is higher than would have been acquired from freefall alone (which would just be GM/r), i.e. solving the equations:
     8First, I will consider cases where the kinetic energy was greater than GM/r. This corresponds to the lower plot. I can solve the following equations for r, given various ke(r),
     9
    610
    711[[latex($ke(r)-\frac{G*M}{r}+\frac{G*M}{R} = 1.01 (\frac{GM}{R})$)]]
     
    913[[latex($ke(r)-\frac{G*M}{r}+\frac{G*M}{R} = 1.3 (\frac{GM}{R})$)]]
    1014
    11 again where the specific kinetic energy at the surface of the star (R) after having fallen from a distance r away from the center of the star is given by,
     15The LHS of these equations give the (specific) accretion energy at the surface of the star, ke(R), from gas that fell from a distance r away, i.e.,
    1216
    1317[[latex($ke(R)=ke(r)-\frac{G*M}{r}+\frac{G*M}{R} $)]]
    1418
    1519
     20Now, solving for r given various starting ke(r) and percents of error, we have,
    1621
    1722|| '''ke(r)''' || '''Error (%)''' || '''Distance (pc)'''||
     
    2429
    2530
     31This shows that you get the greatest error closer into the surface of the star, and even when you are going very fast starting from a distance r = .002 you still only have a .01% error. This means that for any cells larger than this r, wouldn't expect great deviation, even in the extreme 1000*freefall case.
     32
    2633Now, what if a gas parcel started from rest, a distance r away from the star surface? Now we are solving,
    2734
     
    3037[[latex($-\frac{G*M}{r}+\frac{G*M}{R} = .70 (\frac{GM}{R})$)]]
    3138
     39for a .01 and 30% error. This gives the following distances,
     40
    3241|| '''ke(r)''' || '''Error (%)''' || '''Distance (pc)'''||
    3342|| 0 || .01 || 2.2*10^-6^ ||
    3443|| 0 || 30 || 7.5*10^-8^ ||
    3544
     45Again, the closer into the star we get, the worse the approximation gets -- but that a cell size of r~ 2.2*10^-6^ shouldn't produce much deviation.
    3646
    37 Lastly, what if the parcel was moving, however, it was moving ''slower'' than freefall? Now ke(r) will be a fraction of the freefall energy in the table below. In particular, what if ke(r) was 1/2, 1/5, 1/10 the freefall kinetic energy... at what distance, r, would I see a 0.01% error? What about a 30% error?
     47
     48Lastly, what if the parcel was moving ''slower'' than freefall? Now ke(r) will be a fraction of the freefall energy in the table below. In particular, what if ke(r) was 1/2, 1/5, 1/10 the freefall kinetic energy... at what distance, r, would we see a 0.01% error? What about a 30% error?
    3849
    3950[[latex($ke(r)-\frac{G*M}{r}+\frac{G*M}{R} = .99 (\frac{GM}{R})$)]]
     
    4859|| 1/10*freefall || .01 || 2*10^-6 ||
    4960|| 1/10*freefall || 30 || 6.8*10^-8^ ||
     61
     62
     63This all shows us that '''the majority of energy gains from a stellar source occurs very close to the source''', and thus, using the equation of accretion luminosity for a particle starting from rest at infinity is a fine approximation.