hub/city_model_structure/building_demand/thermal_zone.py

"""
ThermalZone module
SPDX - License - Identifier: LGPL - 3.0 - or -later
Copyright © 2022 Concordia CERC group
Project Coder Guille Gutierrez guillermo.gutierrezmorote@concordia.ca
Code contributors: Pilar Monsalvete Alvarez de Uribarri pilar.monsalvete@concordia.ca
"""

import uuid
import copy
import numpy
from typing import List, Union, TypeVar
from city_model_structure.building_demand.occupancy import Occupancy
from city_model_structure.building_demand.appliances import Appliances
from city_model_structure.building_demand.lighting import Lighting
from city_model_structure.building_demand.internal_gain import InternalGain
from city_model_structure.building_demand.thermal_control import ThermalControl
from city_model_structure.energy_systems.hvac_system import HvacSystem
from city_model_structure.attributes.schedule import Schedule
import helpers.constants as cte

ThermalBoundary = TypeVar('ThermalBoundary')
InternalZone = TypeVar('InternalZone')


class ThermalZone:
  """
  ThermalZone class
  """
  def __init__(self, thermal_boundaries, parent_internal_zone, volume, footprint_area, usage=None):
    self._id = None
    self._parent_internal_zone = parent_internal_zone
    self._footprint_area = footprint_area
    self._thermal_boundaries = thermal_boundaries
    self._additional_thermal_bridge_u_value = None
    self._effective_thermal_capacity = None
    self._indirectly_heated_area_ratio = None
    self._infiltration_rate_system_on = None
    self._infiltration_rate_system_off = None
    self._volume = volume
    self._ordinate_number = None
    self._view_factors_matrix = None
    self._total_floor_area = None
    self._usage = usage
    self._usage_from_parent = False
    if usage is None:
      self._usage_from_parent = True
    self._hours_day = None
    self._days_year = None
    self._mechanical_air_change = None
    self._occupancy = None
    self._lighting = None
    self._appliances = None
    self._internal_gains = None
    self._thermal_control = None
    self._usage_zones = None

  @property
  def usage_zones(self):
    # example 70-office_30-residential
    if self._usage_from_parent:
      self._usage_zones = copy.deepcopy(self._parent_internal_zone.usage_zones)
    else:
      values = self._usage.split('_')
      usages = []
      for value in values:
        usages.append(value.split('-'))
      self._usage_zones = []
      for parent_usage_zone in self._parent_internal_zone.usage_zones:
        for value in usages:
          if parent_usage_zone.usage == value[1]:
            new_usage_zone = copy.deepcopy(parent_usage_zone)
            new_usage_zone.percentage = float(value[0])/100
            self._usage_zones.append(new_usage_zone)
    return self._usage_zones

  @property
  def id(self):
    """
    Get thermal zone id, a universally unique identifier randomly generated
    :return: str
    """
    if self._id is None:
      self._id = uuid.uuid4()
    return self._id

  @property
  def footprint_area(self) -> float:
    """
    Get thermal zone footprint area in m2
    :return: float
    """
    return self._footprint_area

  @property
  def thermal_boundaries(self) -> [ThermalBoundary]:
    """
    Get thermal boundaries bounding with the thermal zone
    :return: [ThermalBoundary]
    """
    return self._thermal_boundaries

  @property
  def additional_thermal_bridge_u_value(self) -> Union[None, float]:
    """
    Get thermal zone additional thermal bridge u value W/m2K
    :return: None or float
    """
    return self._additional_thermal_bridge_u_value

  @additional_thermal_bridge_u_value.setter
  def additional_thermal_bridge_u_value(self, value):
    """
    Set thermal zone additional thermal bridge u value W/m2K
    :param value: float
    """
    if value is not None:
      self._additional_thermal_bridge_u_value = float(value)

  @property
  def effective_thermal_capacity(self) -> Union[None, float]:
    """
    Get thermal zone effective thermal capacity in J/m2K
    :return: None or float
    """
    return self._effective_thermal_capacity

  @effective_thermal_capacity.setter
  def effective_thermal_capacity(self, value):
    """
    Set thermal zone effective thermal capacity in J/m2K
    :param value: float
    """
    if value is not None:
      self._effective_thermal_capacity = float(value)

  @property
  def indirectly_heated_area_ratio(self) -> Union[None, float]:
    """
    Get thermal zone indirectly heated area ratio
    :return: None or float
    """
    return self._indirectly_heated_area_ratio

  @indirectly_heated_area_ratio.setter
  def indirectly_heated_area_ratio(self, value):
    """
    Set thermal zone indirectly heated area ratio
    :param value: float
    """
    if value is not None:
      self._indirectly_heated_area_ratio = float(value)

  @property
  def infiltration_rate_system_on(self):
    """
    Get thermal zone infiltration rate system on in air changes per hour (ACH)
    :return: None or float
    """
    return self._infiltration_rate_system_on

  @infiltration_rate_system_on.setter
  def infiltration_rate_system_on(self, value):
    """
    Set thermal zone infiltration rate system on in air changes per hour (ACH)
    :param value: float
    """
    self._infiltration_rate_system_on = value

  @property
  def infiltration_rate_system_off(self):
    """
    Get thermal zone infiltration rate system off in air changes per hour (ACH)
    :return: None or float
    """
    return self._infiltration_rate_system_off

  @infiltration_rate_system_off.setter
  def infiltration_rate_system_off(self, value):
    """
    Set thermal zone infiltration rate system on in air changes per hour (ACH)
    :param value: float
    """
    self._infiltration_rate_system_off = value

  @property
  def volume(self):
    """
    Get thermal zone volume
    :return: float
    """
    return self._volume

  @property
  def ordinate_number(self) -> Union[None, int]:
    """
    Get the order in which the thermal_zones need to be enumerated
    :return: None or int
    """
    return self._ordinate_number

  @ordinate_number.setter
  def ordinate_number(self, value):
    """
    Set a specific order of the zones to be called
    :param value: int
    """
    if value is not None:
      self._ordinate_number = int(value)

  @property
  def hvac_system(self) -> Union[None, HvacSystem]:
    """
    Get HVAC system installed for this thermal zone
    From internal_zone
    :return: None or HvacSystem
    """
    return self._parent_internal_zone.hvac_system

  @property
  def view_factors_matrix(self):
    """
    Get thermal zone view factors matrix
    :return: [[float]]
    """
    return self._view_factors_matrix

  @view_factors_matrix.setter
  def view_factors_matrix(self, value):
    """
    Set thermal zone view factors matrix
    :param value: [[float]]
    """
    self._view_factors_matrix = value

  @property
  def usage(self) -> Union[None, str]:
    """
    Get thermal zone usage
    :return: None or str
    """
    if self._usage_from_parent:
      if self._parent_internal_zone.usage_zones is None:
        return None
      self._usage = ''
      for usage_zone in self._parent_internal_zone.usage_zones:
        self._usage += str(round(usage_zone.percentage * 100)) + '-' + usage_zone.usage + '_'
      self._usage = self._usage[:-1]
    return self._usage

  @staticmethod
  def _get_schedule_of_day(requested_day_type, schedules):
    for schedule in schedules:
      for day_type in schedule.day_types:
        if day_type == requested_day_type:
          return schedule
      return None

  @property
  def hours_day(self) -> Union[None, float]:
    """
    Get thermal zone usage hours per day
    :return: None or float
    """
    if self.usage_zones is None:
      return None
    if self._hours_day is None:
      self._hours_day = 0
      for usage_zone in self.usage_zones:
        self._hours_day += usage_zone.percentage * usage_zone.hours_day
    return self._hours_day

  @property
  def days_year(self) -> Union[None, float]:
    """
    Get thermal zone usage days per year
    :return: None or float
    """
    if self.usage_zones is None:
      return None
    if self._days_year is None:
      self._days_year = 0
      for usage_zone in self.usage_zones:
        self._days_year += usage_zone.percentage * usage_zone.days_year
    return self._days_year

  @property
  def mechanical_air_change(self) -> Union[None, float]:
    """
    Get thermal zone mechanical air change in air change per hour (ACH)
    :return: None or float
    """
    if self.usage_zones is None:
      return None
    if self._mechanical_air_change is None:
      self._mechanical_air_change = 0
      for usage_zone in self.usage_zones:
        if usage_zone.mechanical_air_change is None:
          return None
        self._mechanical_air_change += usage_zone.percentage * usage_zone.mechanical_air_change
    return self._mechanical_air_change

  @property
  def occupancy(self) -> Union[None, Occupancy]:
    """
    Get occupancy in the thermal zone
    :return: None or Occupancy
    """
    if self.usage_zones is None:
      return None

    self._occupancy = Occupancy()
    _occupancy_density = 0
    _convective_part = 0
    _radiative_part = 0
    _latent_part = 0
    for usage_zone in self.usage_zones:
      if usage_zone.occupancy is None:
        return None
      _occupancy_density += usage_zone.percentage * usage_zone.occupancy.occupancy_density
      if usage_zone.occupancy.sensible_convective_internal_gain is not None:
        _convective_part += usage_zone.percentage * usage_zone.occupancy.sensible_convective_internal_gain
        _radiative_part += usage_zone.percentage * usage_zone.occupancy.sensible_radiative_internal_gain
        _latent_part += usage_zone.percentage * usage_zone.occupancy.latent_internal_gain
    self._occupancy.occupancy_density = _occupancy_density
    self._occupancy.sensible_convective_internal_gain = _convective_part
    self._occupancy.sensible_radiative_internal_gain = _radiative_part
    self._occupancy.latent_internal_gain = _latent_part

    _occupancy_reference = self.usage_zones[0].occupancy
    if _occupancy_reference.occupancy_schedules is not None:
      _schedules = []
      for i_schedule in range(0, len(_occupancy_reference.occupancy_schedules)):
        schedule = copy.deepcopy(_occupancy_reference.occupancy_schedules[i_schedule])
        new_values = []
        for i_value in range(0, len(_occupancy_reference.occupancy_schedules[i_schedule].values)):
          _new_value = 0
          for usage_zone in self.usage_zones:
            _new_value += usage_zone.percentage * usage_zone.occupancy.occupancy_schedules[i_schedule].values[i_value]
          new_values.append(_new_value)
        schedule.values = new_values
        _schedules.append(schedule)
      self._occupancy.occupancy_schedules = _schedules
    return self._occupancy

  @property
  def lighting(self) -> Union[None, Lighting]:
    """
    Get lighting information
    :return: None or Lighting
    """
    if self.usage_zones is None:
      return None

    self._lighting = Lighting()
    _lighting_density = 0
    _convective_part = 0
    _radiative_part = 0
    _latent_part = 0
    for usage_zone in self.usage_zones:
      if usage_zone.lighting is None:
        return None
      _lighting_density += usage_zone.percentage * usage_zone.lighting.density
      if usage_zone.lighting.convective_fraction is not None:
        _convective_part += usage_zone.percentage * usage_zone.lighting.density \
                             * usage_zone.lighting.convective_fraction
        _radiative_part += usage_zone.percentage * usage_zone.lighting.density \
                           * usage_zone.lighting.radiative_fraction
        _latent_part += usage_zone.percentage * usage_zone.lighting.density \
                        * usage_zone.lighting.latent_fraction
    self._lighting.density = _lighting_density
    if _lighting_density > 0:
      self._lighting.convective_fraction = _convective_part / _lighting_density
      self._lighting.radiative_fraction = _radiative_part / _lighting_density
      self._lighting.latent_fraction = _latent_part / _lighting_density
    else:
      self._lighting.convective_fraction = 0
      self._lighting.radiative_fraction = 0
      self._lighting.latent_fraction = 0

    _lighting_reference = self.usage_zones[0].lighting
    if _lighting_reference.schedules is not None:
      _schedules = []
      for i_schedule in range(0, len(_lighting_reference.schedules)):
        schedule = copy.deepcopy(_lighting_reference.schedules[i_schedule])
        new_values = []
        for i_value in range(0, len(_lighting_reference.schedules[i_schedule].values)):
          _new_value = 0
          for usage_zone in self.usage_zones:
            _new_value += usage_zone.percentage * usage_zone.lighting.schedules[i_schedule].values[i_value]
          new_values.append(_new_value)
        schedule.values = new_values
        _schedules.append(schedule)
      self._lighting.schedules = _schedules
    return self._lighting

  @property
  def appliances(self) -> Union[None, Appliances]:
    """
    Get appliances information
    :return: None or Appliances
    """
    if self.usage_zones is None:
      return None

    self._appliances = Appliances()
    _appliances_density = 0
    _convective_part = 0
    _radiative_part = 0
    _latent_part = 0
    for usage_zone in self.usage_zones:
      if usage_zone.appliances is None:
        return None
      _appliances_density += usage_zone.percentage * usage_zone.appliances.density
      if usage_zone.appliances.convective_fraction is not None:
        _convective_part += usage_zone.percentage * usage_zone.appliances.density \
                            * usage_zone.appliances.convective_fraction
        _radiative_part += usage_zone.percentage * usage_zone.appliances.density \
                           * usage_zone.appliances.radiative_fraction
        _latent_part += usage_zone.percentage * usage_zone.appliances.density \
                        * usage_zone.appliances.latent_fraction
    self._appliances.density = _appliances_density
    if _appliances_density > 0:
      self._appliances.convective_fraction = _convective_part / _appliances_density
      self._appliances.radiative_fraction = _radiative_part / _appliances_density
      self._appliances.latent_fraction = _latent_part / _appliances_density
    else:
      self._appliances.convective_fraction = 0
      self._appliances.radiative_fraction = 0
      self._appliances.latent_fraction = 0

    _appliances_reference = self.usage_zones[0].appliances
    if _appliances_reference.schedules is not None:
      _schedules = []
      for i_schedule in range(0, len(_appliances_reference.schedules)):
        schedule = copy.deepcopy(_appliances_reference.schedules[i_schedule])
        new_values = []
        for i_value in range(0, len(_appliances_reference.schedules[i_schedule].values)):
          _new_value = 0
          for usage_zone in self.usage_zones:
            _new_value += usage_zone.percentage * usage_zone.appliances.schedules[i_schedule].values[i_value]
          new_values.append(_new_value)
        schedule.values = new_values
        _schedules.append(schedule)
      self._appliances.schedules = _schedules
    return self._appliances

  @property
  def internal_gains(self) -> Union[None, List[InternalGain]]:
    """
    Calculates and returns the list of all internal gains defined
    :return: [InternalGain]
    """
    if self.usage_zones is None:
      return None

    _internal_gain = InternalGain()
    _days = [cte.MONDAY, cte.TUESDAY, cte.WEDNESDAY, cte.THURSDAY, cte.FRIDAY, cte.SATURDAY, cte.SUNDAY, cte.HOLIDAY]
    _average_internal_gain = 0
    _convective_fraction = 0
    _radiative_fraction = 0
    _latent_fraction = 0
    _schedules = None
    _base_schedule = Schedule()
    _base_schedule.type = cte.INTERNAL_GAINS
    _base_schedule.time_range = cte.DAY
    _base_schedule.time_step = cte.HOUR
    _base_schedule.data_type = cte.ANY_NUMBER
    _schedules_defined = True
    values = numpy.zeros(24, 8)
    for usage_zone in self.usage_zones:
      for internal_gain in usage_zone.internal_gains:
        _average_internal_gain += internal_gain.average_internal_gain * usage_zone.percentage
        _convective_fraction += internal_gain.average_internal_gain * usage_zone.percentage \
                                * internal_gain.convective_fraction
        _radiative_fraction += internal_gain.average_internal_gain * usage_zone.percentage \
                               * internal_gain.radiative_fraction
        _latent_fraction += internal_gain.average_internal_gain * usage_zone.percentage \
                            * internal_gain.latent_fraction
    for usage_zone in self.usage_zones:
      for internal_gain in usage_zone.internal_gains:
        if len(internal_gain.schedules) == 0:
          _schedules_defined = False
          break
        for day, _schedule in enumerate(internal_gain.schedules):
          for v, value in enumerate(_schedule.values):
            values[v, day] += value * usage_zone.percentage

    if _schedules_defined:
      _schedules = []
      for day in range(0, len(_days)):
        _schedule = copy.deepcopy(_base_schedule)
        _schedule.day_types = [_days[day]]
        _schedule.values = values[:day]
        _schedules.append(_schedule)

    _internal_gain.convective_fraction = _convective_fraction / _average_internal_gain
    _internal_gain.radiative_fraction = _radiative_fraction / _average_internal_gain
    _internal_gain.latent_fraction = _latent_fraction / _average_internal_gain
    _internal_gain.average_internal_gain = _average_internal_gain
    _internal_gain.type = 'mean_value'
    _internal_gain.schedules = _schedules
    self._internal_gains = [_internal_gain]
    return self._internal_gains

  @property
  def thermal_control(self) -> Union[None, ThermalControl]:
    """
    Get thermal control of this thermal zone
    :return: None or ThermalControl
    """
    if self.usage_zones is None:
      return None

    self._thermal_control = ThermalControl()
    _mean_heating_set_point = 0
    _heating_set_back = 0
    _mean_cooling_set_point = 0
    for usage_zone in self.usage_zones:
      _mean_heating_set_point += usage_zone.percentage * usage_zone.thermal_control.mean_heating_set_point
      _heating_set_back += usage_zone.percentage * usage_zone.thermal_control.heating_set_back
      _mean_cooling_set_point += usage_zone.percentage * usage_zone.thermal_control.mean_cooling_set_point
    self._thermal_control.mean_heating_set_point = _mean_heating_set_point
    self._thermal_control.heating_set_back = _heating_set_back
    self._thermal_control.mean_cooling_set_point = _mean_cooling_set_point

    _thermal_control_reference = self.usage_zones[0].thermal_control
    _types_reference = []
    if _thermal_control_reference.hvac_availability_schedules is not None:
      _types_reference.append([cte.HVAC_AVAILABILITY, _thermal_control_reference.hvac_availability_schedules])
    if _thermal_control_reference.heating_set_point_schedules is not None:
      _types_reference.append([cte.HEATING_SET_POINT, _thermal_control_reference.heating_set_point_schedules])
    if _thermal_control_reference.cooling_set_point_schedules is not None:
      _types_reference.append([cte.COOLING_SET_POINT, _thermal_control_reference.cooling_set_point_schedules])

    for i_type in range(0, len(_types_reference)):
      _schedules = []
      _schedule_type = _types_reference[i_type][1]
      for i_schedule in range(0, len(_schedule_type)):
        schedule = copy.deepcopy(_schedule_type[i_schedule])
        new_values = []
        for i_value in range(0, len(_schedule_type[i_schedule].values)):
          _new_value = 0
          for usage_zone in self.usage_zones:
            if _types_reference[i_type][0] == cte.HVAC_AVAILABILITY:
              _new_value += usage_zone.percentage * \
                            usage_zone.thermal_control.hvac_availability_schedules[i_schedule].values[i_value]
            elif _types_reference[i_type][0] == cte.HEATING_SET_POINT:
              _new_value += usage_zone.percentage * \
                            usage_zone.thermal_control.heating_set_point_schedules[i_schedule].values[i_value]
            elif _types_reference[i_type][0] == cte.COOLING_SET_POINT:
              _new_value += usage_zone.percentage * \
                            usage_zone.thermal_control.cooling_set_point_schedules[i_schedule].values[i_value]
          new_values.append(_new_value)
        schedule.values = new_values
        _schedules.append(schedule)
      if i_type == 0:
        self._thermal_control.hvac_availability_schedules = _schedules
      elif i_type == 1:
        self._thermal_control.heating_set_point_schedules = _schedules
      elif i_type == 2:
        self._thermal_control.cooling_set_point_schedules = _schedules

    return self._thermal_control

  @property
  def total_floor_area(self):
    """
    Get the total floor area of this thermal zone
    :return: float
    """
    return self._total_floor_area

  @total_floor_area.setter
  def total_floor_area(self, value):
    """
    Set the total floor area of this thermal zone
    :param value: float
    """
    self._total_floor_area = value