Changes between Version 8 and Version 9 of u/erica/scratch4

Timestamp:

02/16/16 17:49:17 (9 years ago)

Author:

Erica Kaminski

Comment:

—

u/erica/scratch4

@@ -3 +3 @@
 = Radiation feedback from sink particles =
 
-The amount of thermal radiation produced in the grid is a function of temperature. Since sinks are a subgrid model they themselves do not have temperature (we are not sure how big the forming star is, how fast it is growing by contraction, etc., so there isn't an easy way of assigning the sub-grid object a 'temperature'). Thus, we are left to estimate the amount of radiation produced by the sink by its accretion luminosity. 
-
-Since we do not track energy accretion onto the sink, we are left to assume gas that is accreted from the surrounding zones contribute to the accretion luminosity of the forming protostar, and thus we need a way to estimate this accretion luminosity. To do this, we will recycle some fraction of infalling energy back into the grid. This ''accretion energy'' ([[latex($E_{acc}$)]]) will be distributed smoothly in a kernel surrounding the sink every time step. It will then diffuse away from the sink via FLD radiative transfer. In this way, sinks will act as additional sources of radiation within the grid. 
+The amount of thermal radiation produced in the grid is a function of temperature. Since sinks are a subgrid model they themselves do not have temperature (we are not sure how big the forming star is, how fast it is growing by contraction, etc., so there isn't an easy way of assigning the sub-grid object a 'temperature'). Thus, we are left to estimate the amount of radiation produced by the sink. We envision the sink particle as a protostar, and not yet producing radiation through fusion. Thus, the radiation our sinks produce comes from accretion.
 
 
-= Accretion Luminosity = 
+== Accretion Luminosity == 
 
-As discussed [https://astrobear.pas.rochester.edu/trac/blog/erica01262016 here], the following equation for the accretion luminosity:
+For spherical symmetry, a gas parcel starting from rest and freely falling to the star from infinity will have its kinetic and gravitational energy balance at the stellar surface:
 
-[[latex($L=\frac{G M dm}{R dt}$)]]
+[[latex($\frac{1}{2} m v_{ff}^2 = \frac{GmM_{*}}{R_{*}}$)]]
 
-should be a very good approximation for Astrobear to use.  
+As material passes through the accretion shocks at the surface of the star, its kinetic energy is converted into heat that is then radiated away. For an accretion rate [[latex($\dot{m}$)]], the rate at which this heat is produced, or the luminosity L, is given by:
 
+[[latex($ L = \frac{1}{2}\dot{m} v_{ff}^2 = \frac{G\dot{m}M_{*}}{R_{*}}$)]]
+
+Since we do not track energy accretion onto the sink, we are left to assume gas that is accreted from the surrounding zones contribute to this accretion luminosity directly. Thus, the best we can do for tracking the energy release from infall is to calculate the RHS of this equation in the code and use it as an estimate of the true accretion luminosity. (By the way, this form of the accretion luminosity was shown to be a good approximation for our purposes [https://astrobear.pas.rochester.edu/trac/blog/erica01262016 here]). 
+
+The accretion energy ([[latex($L*dt$)]]) will be distributed smoothly in a kernel surrounding the sink every time step. It will then diffuse away from the sink via FLD radiative transfer. In this way, sinks will act as additional sources of radiation within the grid. 
 = Tracking accretion luminosity in the code =