========== Decay Heat ========== During the operation of a nuclear reactor, the core's components are heated by interactions of particles (neutrons, photons, electrons, ions and atoms) with component materials. The origin of most of those particles can be traced back in time to an induced fission event in the reactor. So they depend on the reactor power level history. We can further separate these particles into 3 categories: **Immediate** sources Particles that are created immediately after an induced fission event in the reactor (created within a time scale that is much smaller than the generation time). For example: prompt neutrons, prompt photons and fission products. **Delayed** sources Particles that are created after an induced fission event in the reactor with a delay that is larger than the generation time. For example: particles that are created from decay of activation products or fission products. **External** sources Particles that their origin *cannot* be traced back to an induced fission event in the reactor, and so this type of source does not depend on the reactor power level history. For example: particles that are created from an external neutron source, particles that are created by cosmic radiation and particles that are created by spontaneous fission of U238. Most of the time, external sources have a negligible effect on radiation heating. They are not considered here. The total heating must be evaluated to ensure proper heat evacuation from the core (such evaluation is beyond the scope of system codes and particularly :ref:`STREAM `, and serves as input for this subpackage). When the reactor is shut down, the power does not turn off completely, as these delayed sources of heating have a lingering (sometimes called `residual`) effect - **Decay Heat**. ----- The Fission Rate ================ Since the common source of (the majority of) immediate and delayed heating contributions is the induced fission event, it is convenient to write contributions per fission event, and multiply by the rate of fission events. The total decay heat per fission event is given in units of :math:`[\text{MeV} / \text{fission}]`, and denoted :math:`F_\text{Decay Heat}(t, T)` where :math:`T` is the operation time (constant fission rate during this period is assumed), and :math:`t` is the time after shutdown. The decay heat power is then defined as: .. math:: P_\text{Decay Heat} = \Phi F_\text{Decay Heat}(t, T) :label: PDH Where :math:`\Phi` is the **fission rate**, and has units of :math:`[\text{W} / (\text{MeV} / \text{fission})]` or :math:`[\approx 1.602\cdot 10^{-13}\text{ fission}/\text{s}]`, and is commonly written as: .. math:: \Phi = \frac{P_0}{Q} :label: FR Where :math:`P_0` is the reactor operation power and :math:`Q` is the total **recoverable** energy associated with a fission event (conservatively, :math:`Q\approx 200 \text{MeV}`). In the method contained herein this total decay heat per fission event is the sum of several contributions: .. math:: F_\text{Decay Heat}(t, T) = F_\text{Fissions}(t, T) + F_\text{Fission Products}(t, T) + \sum_i F_{\text{Activation Products} i}(t, T) :label: FDH .. note:: The user is left with the responsibility of constructing their own decay heat function, including the fission rate and the weighting of each contribution which depends on the component in question. In the following segments, each term of :eq:`FDH`: :ref:`Fissions`, :ref:`Fission Products`, :ref:`Activation` is explained and the corresponding functions are linked. Fissions ======== With the shutdown of the reactor the neutron population plunges, and with it the fission rate in accordance with the generation time :math:`\Lambda` and the delayed neutron decay rates :math:`\lambda_k`. The fraction of delayed fission yield :math:`\beta_\text{eff}` is spectrum dependent and should, along with the other kinetic parameters, be calculated for a given reactor. The decay heat is calculated as follows: .. math:: F_\text{Fissions}(t, T) = Q_\text{prompt}P(t) :label: fis Where :math:`P(t)` is the profile generated by solving the Point-Kinetics equations, having arbitrary units: .. math:: \dot{P} &= \frac{\rho(t) - \beta}{\Lambda} P + \sum_k C_k\lambda_k \\ \dot{C}_k &= - C_k\lambda_k + \frac{\beta_k}{ \Lambda} P For :math:`P(t=0) = 1` normalized initial power and controlled by the reactivity insertion :math:`\rho(t); \rho(0)=0`. These equations are further described at :class:`~stream.calculations.point_kinetics.PointKinetics`. :math:`Q_\text{prompt}\equiv Q - Q_\text{delayed}` is the total energy deposited in a component directly related to fission events (including, for example, the immediate :math:`\gamma` radiation released `immediately` after fission). The user is provided with methods to generate such a profile: .. autofunction:: stream.physical_models.decay_heat.fissions.profile .. autofunction:: stream.physical_models.decay_heat.fissions.profile_from_pk .. admonition:: Implementation Detail The profiles are generated once and interpolated for use at STREAM simulation-time, in order to save the runtime cost of calculating the profile for each call. Fission Products ================ The primary contributor to reactor decay heat is the :math:`\beta` and :math:`\gamma` decays of fission products. Several standards are readily available, of which the most notable is the American National Standard **ANSI/ANS-5.1**. The supported standards are: * **ANSI/ANS-5.1-1973** [#ANS1973]_ * **ANSI/ANS-5.1-2014** [#ANS2014]_ * **JAERI-M-91-034** [#JAERI1991]_ .. [#ANS1973] American Nuclear Society Proposed Standard, ANS-5.1 "Decay Energy Release Rates Following Shutdown of Uranium-Fueled Thermal Reactors," October (1971), Revised October (1973). .. [#ANS2014] "American National Standard for Decay Heat Power in Light Water Reactors", ANSI/ANS-5.1-2014, American Nuclear Society, 2014. .. [#JAERI1991] K. Tasaka, "Recommended values of decay heat power and method to utilize the data", Japan Atomic Energy Research Institute, Japan, JAERI-M-91-034, 1991. The american standards do not distinguish between :math:`\beta` and :math:`\gamma` contributions and are only meant for fuel elements. Since the :math:`\beta` contribution is naturally local while the :math:`\gamma` contribution may not be, for non-fuel elements only the latter profile should be used, a distinction available through JAERI [#JAERI1991]_, included mainly for this reason. For these cases, the user must normalize the profile and multiply by the correct fission product contribution evaluated at shutdown. The user may choose a standard and a source (for the reason outlined above), or provide their own constants for the general equation: .. math:: F_\text{Fission Products}(t, T) = \sum_i\frac{\alpha_i}{\lambda_i} e^{-\lambda_i t}\left(1 -e^{-\lambda_i T}\right) :label: FP Where :math:`\alpha_i \space [\text{MeV}/(\text{fission}/s)]` and :math:`\lambda_i \space [1/s]` are given by the standards. .. autoclass:: stream.physical_models.decay_heat.fission_products.Standard :members: :undoc-members: .. autoclass:: stream.physical_models.decay_heat.fission_products.Source :members: :undoc-members: .. autofunction:: stream.physical_models.decay_heat.fission_products.contribution .. autofunction:: stream.physical_models.decay_heat.fission_products.fp_inner_ Neutron Captures in Fission Products ------------------------------------ .. TODO:: Not implemented yet! Activation ========== The decay of activated materials is yet another source of decay heat. For an activated material whose decay rate is :math:`\lambda_1`, the decay heat :math:`F(t, T)` is given by .. math:: F_{\text{Activation }}(t, T) = Q_1 e^{-\lambda_1 t}\left(1 -e^{-\lambda_1 T}\right) :label: act1 .. autofunction:: stream.physical_models.decay_heat.activation.profile And if the resultant isotope also decays with decay rate :math:`\lambda_2`,(see :ref:`Actinides`), the solution is .. math:: F_\text{Activation 2}(t, T) = Q_2\left[ \frac{\lambda_1}{\lambda_1 - \lambda_2}e^{-\lambda_2 t} \left(1 -e^{-\lambda_2 T}\right) - \frac{\lambda_2}{\lambda_1 - \lambda_2}e^{-\lambda_1 t} \left(1 -e^{-\lambda_1 T}\right)\right] :label: act2 .. autofunction:: stream.physical_models.decay_heat.activation.double_decay_profile The profiles provided here should be factored by the appropriate :math:`Q` values, as seen in equations :eq:`act1`, :eq:`act2`. One way of doing that is by writing :math:`Q = E_d R` where :math:`E_d \ [\text{MeV}]` is the deposited energy per decay event, and :math:`R \ [1/\text{fission}]` is the number of capture events per fission event which lead to the creation of the decaying isotope. Actinides --------- The aforementioned ANS standard [#ANS1973]_ specifically mentions the contribution of :math:`^{239} U` and :math:`^{239} Np` which are created and decay by the following process: .. math:: {}^{238} U \xrightarrow{n \text{ capture}} {}^{239} U \xrightarrow{\beta^-} {}^{239} Np \xrightarrow{\beta^-} {}^{239} Pu Their contributions are given by: .. math:: F_{U}(t, T) &= E_{U}R_{U}e^{-\lambda_1 t}\left(1 -e^{-\lambda_1 T}\right) \\ F_{Np}(t, T) &= E_{Np}R_{U}\left[ \frac{\lambda_1}{\lambda_1 - \lambda_2}e^{-\lambda_2 t} \left(1 -e^{-\lambda_2 T}\right) - \frac{\lambda_2}{\lambda_1 - \lambda_2}e^{-\lambda_1 t} \left(1 -e^{-\lambda_1 T}\right)\right] :label: actinides_contribution .. automodule:: stream.physical_models.decay_heat.actinides :members: