bluemira.balance_of_plant.calculations
======================================

.. py:module:: bluemira.balance_of_plant.calculations

.. autoapi-nested-parse::

   Simple relations for power.



Functions
---------

.. autoapisummary::

   bluemira.balance_of_plant.calculations.cryo_power
   bluemira.balance_of_plant.calculations.He_pumping
   bluemira.balance_of_plant.calculations.H2O_pumping
   bluemira.balance_of_plant.calculations.superheated_rankine


Module Contents
---------------

.. py:function:: cryo_power(s_tf: float, m_cold: float, nucl_heating: float, e_pf_max: float, t_pulse: float, tf_current: float, n_TF: int) -> float

   Calculates cryogenic loads (taken from PROCESS)

   :param s_tf: TF coil total surface area [m^2]
   :param m_cold: Total cold mass [kg]
   :param nucl_heating: Total coil nuclear heating [W]
   :param e_pf_max: Maximum stored energy in the PF coils [J]
   :param t_pulse: Pulse length [s]
   :param tf_current: TF coil current per turn [A]
   :param n_TF: Number of TF coils

   :rtype: Total power required to cool cryogenic components

   .. note::

      Author: P J Knight, CCFE, Culham Science Centre
      D. Slack memo SCMDG 88-5-1-059, LLNL ITER-88-054, Aug. 1988


.. py:function:: He_pumping(pressure_in: float, pressure_out: float, t_in: float, t_out: float, blanket_power: float, eta_isen: float, eta_el: float) -> tuple[float, float]

   Calculate the pumping power for helium-cooled blankets.

   :param pressure_in: Inlet pressure [Pa]
   :param pressure_out: Pressure drop [Pa]
   :param t_in: Inlet temperature [K]
   :param t_out: Outlet temperature [K]
   :param blanket_power: Total blanket power excluding pumping power [W]
   :param eta_isen: Isentropic efficiency of the He compressors
   :param eta_el: Electrical efficiency of the He compressors

   :returns: * *P_pump_is* -- The isentropic pumping power (added to the working fluid) [W]
             * *P_pump_el* -- The electrical pumping power (parasitic load) [W]

               :math:`T_{in_{comp}} = \dfrac{T_{in_{BB}}}{\dfrac{P}{P-dP}^{\dfrac{\gamma-1}{\gamma}}}`

               :math:`f_{p} = \dfrac{T_{in_{comp}}}{\eta_{is}dT}\Bigg(\dfrac{P}{P-dP}^{\dfrac{\gamma-1}{\gamma}}-1\Bigg)`

   .. rubric:: Notes

   :math:`f_{p} = \dfrac{T_{in_{BB}}}{\eta_{is}dT}\Bigg(1-\dfrac{P-dP}{P}^{\dfrac{\gamma-1}{\gamma}}\Bigg)`

   **Outputs:**

       :math:`P_{pump} = \dfrac{f_{p}P_{plasma}}{1-f_p}` [W]

       :math:`P_{pump,el} = \dfrac{P_{pump}}{\eta_{el}}` [W]

   **No longer in use:**
       :math:`f_{pump}=\dfrac{dP}{dTc_P\rho_{av}}`


.. py:function:: H2O_pumping(p_blanket: float, f_pump: float, eta_isen: float, eta_el: float) -> tuple[float, float]

   H20-cooling pumping power calculation strategy

   :param f_pump: Fraction of thermal power required to pump
   :param eta_isen: Isentropic efficiency of the water pumps
   :param eta_el: Electrical efficiency of the water pumps
   :type eta_el: t

   :returns: * *P_pump_is* -- The isentropic pumping power (added to the working fluid)
             * *P_pump_el* -- The eletrical pumping power (parasitic load)

   .. rubric:: Notes

   The equations used in this function are:

   .. math::
       f_{pump} = \frac{f_{pump}}{\eta_{isen}}

   .. math::
       P_{pump\_is} = \frac{f_{pump} \cdot p_{blanket}}{1 - f_{pump}}

   .. math::
       P_{pump\_el} = \frac{P_{pump\_is}}{\eta_{el}}


.. py:function:: superheated_rankine(blanket_power: float, div_power: float, bb_outlet_temp: float, delta_t_turbine: float) -> float

   PROCESS C. Harrington correlation. Accounts for low-grade heat penalty.
   Used for He-cooled blankets. Not applicable to H2O temperatures.

   :param blanket_power: Blanket thermal power [W]
   :param div_power: Divertor thermal power [W]
   :param bb_outlet_temp: Blanket outlet temperature [K]
   :param delta_t_turbine: Turbine inlet temperature drop [K]

   :rtype: Efficiency of a superheated Rankine cycle

   .. rubric:: Notes

   The equations used in this function are:

   .. math::
       t_{turb} = bb_{outlet\_temp} - \delta_{t\_turbine}

   .. math::
       f_{lgh} = \frac{div_{power}}{blanket_{power} + div_{power}}

   .. math::
       \delta_{\eta} = 0.339 \cdot f_{lgh}

   .. math::
       Efficiency = 0.1802 \cdot log(t_{turb}) - 0.7823 - \delta_{\eta}