lca_carbon_workflow/lca_carbon_workflow.py
2024-07-19 14:08:32 -04:00

98 lines
3.9 KiB
Python

"""
lca_carbon_workflow module
Returns the summarize of envelope and energy systems
SPDX - License - Identifier: LGPL - 3.0 - or -later
Copyright © 2024 Concordia CERC group
Code contributors: Alireza Adli alireza.adli@concordia.ca
Mohammad Reza Seyedabadi mohammad.seyedabadi@mail.concordia.ca
"""
from pathlib import Path
from hub.imports.geometry_factory import GeometryFactory
from hub.imports.construction_factory import ConstructionFactory
from hub.helpers.dictionaries import Dictionaries
from city_model_structure.life_cycle_assessment.access_nrcan_catalogue \
import AccessNrcanCatalog
from city_model_structure.life_cycle_assessment.opening_emission \
import OpeningEmission
class LCACarbonWorkflow:
def __init__(
self,
city_path,
archetypes_catalog_file_name,
constructions_catalog_file,
catalog='nrcan',
building_parameters=('height', 'year_of_construction', 'function')):
self.file_path = (Path(__file__).parent / 'input_files' / city_path)
self.catalogs_path = (Path(__file__).parent / 'input_files')
self.archetypes_catalog_file_name = archetypes_catalog_file_name
self.constructions_catalog_file = constructions_catalog_file
self.nrcan_catalogs = AccessNrcanCatalog(
self.catalogs_path,
archetypes=self.archetypes_catalog_file_name,
constructions=self.constructions_catalog_file)
self.out_path = (Path(__file__).parent / 'out_files')
self.handler = catalog
self.height, self.year_of_construction, self.function = \
building_parameters
print('[simulation start]')
self.city = GeometryFactory(
'geojson',
path=self.file_path,
height_field=self.height,
year_of_construction_field=self.year_of_construction,
function_field=self.function,
function_to_hub=Dictionaries().montreal_function_to_hub_function).city
print(f'city created from {self.file_path}')
ConstructionFactory(self.handler, self.city).enrich()
def calculate_opening_emission(self, building):
surface_opening_emission = []
period_of_construction = self.nrcan_catalogs.year_to_period_of_construction(building.year_of_construction)
print(period_of_construction)
print(building.year_of_construction)
opaque_surface_code = self.nrcan_catalogs.find_opaque_surface(
self.nrcan_catalogs.hub_to_nrcan_function(building.function),
self.nrcan_catalogs.year_to_period_of_construction(
building.year_of_construction),
'6')
for surface in building.surfaces:
boundary_opening_emission = []
for boundary in surface.associated_thermal_boundaries:
opening_emission = []
if boundary.window_ratio == 0:
opening_surface = 0
opening_material_emission = 0
# print('window = 0')
else:
print('window != 0')
for opening in boundary.thermal_openings:
transparent_surface_type = 'Window'
if building.year_of_construction >= 2020 and \
surface.type == 'Roof':
transparent_surface_type = 'Skylight'
opening_emission.append(OpeningEmission(
self.nrcan_catalogs.search_transparent_surfaces(
transparent_surface_type, opaque_surface_code)['embodied_carbon'],
opening.area).calculate_opening_emission())
print('\t\t\t', opening_emission)
print('opening_emission', opening_emission)
if opening_emission:
boundary_opening_emission += opening_emission
print('boundary_opening_emission', boundary_opening_emission)
if boundary_opening_emission:
surface_opening_emission += boundary_opening_emission
print('surface_opening_emission', surface_opening_emission, '\n', sum(surface_opening_emission))
def calculate_emission(self):
emitted_carbon = []
opening_emission = []
for building in self.city.buildings:
self.calculate_opening_emission(building)