Context Navigation

scratch

-              v17
+              v18
 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,
+Sum dE_acc | cells = E_acc
+[[latex($\Sigma ~dE_{acc} = E_{acc}$)]]
 In other words, the exponential function is 'normalized'.
+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 set of equations,
+So given these 2 constraints, we have the following equation,
+[[latex($E_{acc} = \sum_{i=0}^{4} k(~e^{-dx*i} - e^{-4dx} )~E_{acc}  $)]]
+which gives the normalization constant,
+[[latex($ \boxed{k = \frac{1}{~\Sigma ~e^{-dx*i}-e^{-4dx}}} $)]]
+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, but will be more general in the code).