Changes between Version 16 and Version 17 of u/erica/scratch4

Timestamp:

02/16/16 19:20:46 (9 years ago)

Author:

Erica Kaminski

Comment:

—

u/erica/scratch4

@@ -39 +39 @@
 [[latex($\Delta E_i= k E e^{-r_i/\sigma} ~~~~~~~(2)$)]]
 
-where [[latex($\sigma$)]] is a scaling factor. In the code this scaling factor covers 1 e-folding over the kernel. There is also an additional factor of e not shown here that  lets the exponential fall to zero at the boundary of the kernel, and a width that gives ... how to describe this. Now, to solve for the normalization constant, we insert (2) into (1):
+where [[latex($\sigma$)]] is a scaling factor. In the code the exponential function used is such that it falls smoothly to zero at the boundary of the kernel, and covers 4 e-foldings. Now, to solve for the normalization constant, we insert (2) into (1):
 
 [[latex($\sum k E e^{-r_i/\sigma}dV_i=E$)]]
@@ -55 +55 @@
 == Work arrays and feeding L into the radiative source function ==
 
+There is a subroutine in the code that is called by the radiative transfer module, 'apply kernel to work array'. This goes through and populates the work array index for each cell's accretion energy as:
+
+
+[[latex($ info\%workarray(i,j,k,iaccretion\_energy= k*E*e^{-r/\sigma}/dt $)]]
+
+This is an average luminosity over the hydro step. This is then fed into the source function as,
+
+[[latex($ source(i,j,k,iaccretion\_energy)= info\%workarray(i,j,k,iaccretion\_energy)*dt_{rad} $)]]
+
+where we now enter the total ''specific'' (recall k has units of 1/volume) accretion energy into the source function averaged over a radiative timestep. And now, we are ready to begin testing.
+
 = Tests =