Changes between Version 1 and Version 2 of 1DPulsedJets


Ignore:
Timestamp:
01/21/12 12:56:02 (13 years ago)
Author:
ehansen
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • 1DPulsedJets

    v1 v2  
    44
    55== The Physics of Z Cooling ==
    6 
     6Z cooling is an energy loss due to forbidden line emission.  The routines use Pat Hartigan's cooling table Zcooling.tab.  The table does not include recombination or permitted line emission.
    77[[BR]]
    88
    99== Zcooling.tab ==
     10Understanding the format of this table is crucial to interpolating correctly.  There are a few properties that this table must have:
     111.) It must be ordered as log(ne), T, log(X), (9 species cooling fractions), cooling rate. Where ne is the electron density, T is temperature, X is ionization fraction, and the 9 species cooling fractions are the fractions of how much each species contributes to the cooling rate. The 9 species are ordered as follows: OI, OII, NI, NII, SII, FeII, SIII, MgII, CII.
     122.) The first three parameters in the table must be "cycled" through in the following order: log(X), log(ne), T.  In other words, if you follow the table going down, the log(X) column changes first while log(ne) and T remain constant.  Then once all values of log(X) are cycled through, log(ne) changes to its next value and log(X) gets cycled through again.  The process repeats until all possible values are cycled through.  This is important so that all the cooling rates
     133.) The very first line of the table needs to have nDensities, nTemps, nXs.  As an example, the first line might read: 40 30 21.  This means there are 40 different values for log(ne), 30 for T, and 21 for log(X).  Thus, there are a total of 25,200 different values for the cooling rate (40*30*21 = 25,200).
    1014
     15Also, it is important to note that the cooling rate has units of erg*cm^-3*nH^-2.  So before the rate can be added to the change in energy, dqdt(iE), it must be multiplied by nH^2 where nH is the number density of hydrogen.
    1116[[BR]]
    1217