diff --git a/armi/physics/neutronics/globalFlux/globalFluxInterface.py b/armi/physics/neutronics/globalFlux/globalFluxInterface.py index a85723a94..b34acbfc8 100644 --- a/armi/physics/neutronics/globalFlux/globalFluxInterface.py +++ b/armi/physics/neutronics/globalFlux/globalFluxInterface.py @@ -132,7 +132,7 @@ def checkEnergyBalance(self): global power from the flux solve is computed by summing the block-wise power in the core. This value is then compared to the user-specified power and raises an error if relative difference is - above 1e-5. + above :math:`10^{-5}`. """ powerGenerated = ( self.r.core.calcTotalParam( @@ -1244,35 +1244,30 @@ def computeDpaRate(mgFlux, dpaXs): :id: I_ARMI_FLUX_DPA :implements: R_ARMI_FLUX_DPA - This method calculates DPA rates using the inputted multigroup flux and - DPA cross sections. Is simply does the operation: ``dpaPerSecond = - \sum_{g} mgFlux[g] * dpaXs[g]``. Where ``dpaPerSecond`` is the returned - DPA rate, ``mgFlux`` is the inputted multigroup flux, and ``dpaXS`` is - the inputted DPA cross section. + This method calculates DPA rates using the inputted multigroup flux and DPA cross sections. + Displacements calculated by displacement XS: - Displacements calculated by displacement XS: + .. math:: - .. math:: + \text{Displacement rate} &= \phi N_{\text{HT9}} \sigma \\ + &= (\#/\text{cm}^2/s) \cdot (1/cm^3) \cdot (\text{barn})\\ + &= (\#/\text{cm}^5/s) \cdot \text{(barn)} * 10^{-24} \text{cm}^2/\text{barn} \\ + &= \#/\text{cm}^3/s - \text{Displacement rate} &= \phi N_{\text{HT9}} \sigma \\ - &= (\#/\text{cm}^2/s) \cdot (1/cm^3) \cdot (\text{barn})\\ - &= (\#/\text{cm}^5/s) \cdot \text{(barn)} * 10^{-24} \text{cm}^2/\text{barn} \\ - &= \#/\text{cm}^3/s + :: - :: + DPA rate = displacement density rate / (number of atoms/cc) + = dr [#/cm^3/s] / (nHT9) [1/cm^3] + = flux * barn * 1e-24 - DPA rate = displacement density rate / (number of atoms/cc) - = dr [#/cm^3/s] / (nHT9) [1/cm^3] - = flux * barn * 1e-24 + .. math:: - .. math:: + \frac{\text{dpa}}{s} = \frac{\phi N \sigma}{N} = \phi * \sigma - \frac{\text{dpa}}{s} = \frac{\phi N \sigma}{N} = \phi * \sigma - - the Number density of the structural material cancels out. It's in the macroscopic - XS and in the original number of atoms. + the Number density of the structural material cancels out. It's in the macroscopic + XS and in the original number of atoms. Raises ------ @@ -1334,7 +1329,7 @@ def calcReactionRates(obj, keff, lib): :implements: R_ARMI_FLUX_RX_RATES This method computes 1-group reaction rates for the inputted - :py:class:`ArmiObject `. These + :py:class:`ArmiObject ` These reaction rates include: * fission @@ -1342,26 +1337,26 @@ def calcReactionRates(obj, keff, lib): * n2n * absorption - Scatter could be added as well. This function is quite slow so it is - skipped for now as it is uncommonly needed. + Scatter could be added as well. This function is quite slow so it is + skipped for now as it is uncommonly needed. - Reaction rates are: + Reaction rates are: - .. math:: + .. math:: - \Sigma \phi = \sum_{\text{nuclides}} \sum_{\text{energy}} \Sigma - \phi + \Sigma \phi = \sum_{\text{nuclides}} \sum_{\text{energy}} \Sigma + \phi - The units of :math:`N \sigma \phi` are:: + The units of :math:`N \sigma \phi` are:: - [#/bn-cm] * [bn] * [#/cm^2/s] = [#/cm^3/s] + [#/bn-cm] * [bn] * [#/cm^2/s] = [#/cm^3/s] - The group-averaged microscopic cross section is: + The group-averaged microscopic cross section is: - .. math:: + .. math:: - \sigma_g = \frac{\int_{E g}^{E_{g+1}} \phi(E) \sigma(E) - dE}{\int_{E_g}^{E_{g+1}} \phi(E) dE} + \sigma_g = \frac{\int_{E g}^{E_{g+1}} \phi(E) \sigma(E) + dE}{\int_{E_g}^{E_{g+1}} \phi(E) dE} """ rate = {} for simple in RX_PARAM_NAMES: