Version 19 (modified by 9 years ago) ( diff ) | ,
---|
The accretion luminosity is summed up for all cells surrounding the sink particle, and the total is given by 'E_acc'. This total accretion energy is then smoothed over a kernel of cells surrounding the sink, proportionally to a decreasing exponential, such that,
In other words, the exponential function is normalized.
We also want the exponential function to go to zero at the boundary of the kernel.
So given these 2 constraints, we have the following equation,
which gives the normalization constant,
Note that the sum runs over cells i=0 to 1, where the '0th' cell is the cell the sink is in, and the '4th' cell is the furthest cell from the sink in the kernel, that the exponential goes to zero at this boundary, and dx*i gives a position (which here assumes the sink is at the cell center — the actual function is described below).
Attachments (7)
- accretionluminosity1.png (37.8 KB ) - added by 9 years ago.
- accretionluminosity2.png (25.7 KB ) - added by 9 years ago.
- accretionluminosity.nb (45.8 KB ) - added by 9 years ago.
- 4cellkernel1.png (4.0 KB ) - added by 9 years ago.
- 40cellkernel1.png (4.6 KB ) - added by 9 years ago.
- dx1kernel.png (4.6 KB ) - added by 9 years ago.
- dxpt1kernel.png (3.9 KB ) - added by 9 years ago.
Download all attachments as: .zip