#!/usr/bin/env python # coding: utf-8 # In[2]: import components as cp import numpy.random as nr import numpy as np import Model import pandas as pd import TruncatingFunctions as tf # In[5]: model = Model.Model('nano-TiO2 AUT from 2000 to 2016') # In[6]: # split of the total inflow to different compartments totalInflow = cp.FlowCompartment('Total Inflow', logInflows=True, logOutflows=True) #Import_Manuf = cp.FlowCompartment('Import to Manufacturing', logInflows=True, logOutflows=True, categories=['Import']) #Import_Cons = cp.FlowCompartment('Import to Consumption', logInflows=True, logOutflows=True, categories=['Import']) # In[7]: # Production of the nano-material and Manufacturing of the products it is used in Production = cp.Stock('Production', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['Production']) Manufacture = cp.Stock('Manufacture', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['Manufacture']) Consumption = cp.FlowCompartment('Consumption', logInflows=True, logOutflows=True, categories=['Consumption']) # In[8]: Prod_Air_N = cp.Sink('Pristine from Production to Air', logInflows=True, categories=['Air', 'Air_N', 'Pristine']) Manuf_Air_N = cp.Sink('Pristine from Manufacturing to Air', logInflows=True, categories=['Air', 'Air_N', 'Pristine']) # In[9]: # Manufacturing and waste Manuf_Text = cp.FlowCompartment('Manuf Textiles', logInflows=True, logOutflows=True, categories=['Manufacture']) Manuf_Plas = cp.FlowCompartment('Manuf Plastics', logInflows=True, logOutflows=True, categories=['Manufacture']) manuf_textiles_waste = cp.FlowCompartment('Waste from manuf Textiles', logInflows=True, logOutflows=True, categories=['Manufacture', 'Solid Waste']) manuf_plastics_waste = cp.FlowCompartment('Waste from manuf Plastics', logInflows=True, logOutflows=True, categories=['Manufacture', 'Solid Waste']) ManufTextW_Reuse_P = cp.Sink('Product-embedded from Textile manufacturing to reprocessing and reuse', logInflows=True, categories=['Reprocessing', 'Reuse', 'Reuse_P']) ManufPlasW_Reuse_P = cp.Sink('Product-embedded from Plastic manufacturing to reprocessing and reuse', logInflows=True, categories=['Reprocessing', 'Reuse', 'Reuse_P']) # In[10]: # Definition of Products categories (compartments) PersCare = cp.FlowCompartment('PersCare', logInflows=True, logOutflows=True, categories=['Products']) PersCare_Use = cp.Stock('PersCare_Use', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','Use','PMCrel']) PersCare_EoL = cp.Stock('PersCare_EoL', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','EoL']) Paints = cp.FlowCompartment('Paints', logInflows=True, logOutflows=True, categories='Products') Paints_Use = cp.Stock('Paints_Use', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','Use','PMCrel']) Paints_EoL = cp.Stock('Paints_EoL', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','EoL']) In = cp.FlowCompartment('In', logInflows=True, logOutflows=True, categories='Products') In_Use = cp.Stock('In_Use', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','Use','PMCrel']) In_EoL = cp.Stock('In_EoL', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','EoL']) Cement= cp.FlowCompartment('Cement', logInflows=True, logOutflows=True, categories='Products') Cement_Use =cp.Stock('Cement_Use', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','Use','PMCrel']) Cement_EoL =cp.Stock('Cement_EoL', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','EoL']) Glass = cp.FlowCompartment('Glass', logInflows=True, logOutflows=True, categories='Products') Glass_Use = cp.Stock('Glass_Use', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','Use','PMCrel']) Glass_EoL = cp.Stock('Glass_EoL', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','EoL']) Ceramics = cp.FlowCompartment('Ceramics', logInflows=True, logOutflows=True, categories='Products') Ceramics_Use = cp.Stock('Ceramics_Use', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','Use','PMCrel']) Ceramics_EoL = cp.Stock('Ceramics_EoL', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','EoL']) RubPlas = cp.FlowCompartment('RubPlas', logInflows=True, logOutflows=True, categories='Products') RubPlas_Use = cp.Stock('RubPlas_Use', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','Use','PMCrel']) RubPlas_EoL = cp.Stock('RubPlas_EoL', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','EoL']) Solar = cp.FlowCompartment('Solar', logInflows=True, logOutflows=True, categories='Products') Solar_EoL = cp.Stock('Solar_EoL', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','EoL']) Electronics= cp.FlowCompartment('Electronics', logInflows=True, logOutflows=True, categories='Products') Electronics_Use = cp.Stock('Electronics_Use', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','Use','PMCrel']) Electronics_EoL = cp.Stock('Electronics_EoL', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','EoL']) CleanAgents= cp.FlowCompartment('CleanAgents', logInflows=True, logOutflows=True, categories='Products') CleanAgents_Use = cp.Stock('CleanAgents_Use', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','Use','PMCrel']) CleanAgents_EoL = cp.Stock('CleanAgents_EoL', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','EoL']) Food = cp.FlowCompartment('Food', logInflows=True, logOutflows=True, categories='Products') Food_Use = cp.Stock('Food_Use', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','Use','PMCrel']) Food_EoL = cp.Stock('Food_EoL', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','EoL']) Textiles = cp.FlowCompartment('Textiles', logInflows=True, logOutflows=True, categories='Products') Textiles_Use = cp.Stock('Textiles_Use', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','Use','PMCrel']) Textiles_EoL = cp.Stock('Textiles_EoL', logInflows=True, logOutflows=True, logImmediateFlows=True, categories=['PMC','EoL']) # In[11]: # Definition of forms of release after use PersCare_WW = cp.FlowCompartment('From PersCare to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater']) PersCare_WW_N = cp.FlowCompartment('Pristine in Wastewater from PersCare', logInflows=True, logOutflows=True, categories=['Wastewater_N', 'Pristine']) PersCare_SW = cp.FlowCompartment('From PersCare to Surfacewater', logInflows=True, logOutflows=True, categories=['Surfacewater']) PersCare_SW_N = cp.Sink('Pristine in Surfacewater from PersCare', logInflows=True, categories=['Surfacewater_N', 'Pristine']) Paints_Air = cp.FlowCompartment('From Paints to Air', logInflows=True, logOutflows=True, categories=['Air']) Paints_Air_N = cp.Sink('Pristine in Air from Paints', logInflows=True, categories=['Air_N', 'Pristine']) Paints_Air_M = cp.Sink('Matrix-embedded in Air from Paints', logInflows=True, categories=['Air_M', 'Matrix-embedded']) Paints_WW = cp.FlowCompartment('From Paints to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater']) Paints_WW_N = cp.FlowCompartment('Pristine in Wastewater from Paints', logInflows=True, logOutflows=True, categories=['Wastewater_N', 'Pristine']) Paints_WW_M = cp.FlowCompartment('Matrix-embedded in Wastewater from Paints', logInflows=True, logOutflows=True, categories=['Wastewater_M', 'Matrix-embedded']) Paints_NUsoil = cp.FlowCompartment('From Paints to Soil', logInflows=True, logOutflows=True, categories=['Soil', 'NU Soil']) Paints_NUsoil_N = cp.Sink('Pristine in Soil from Paints', logInflows=True, categories=['Soil_N', 'NU Soil_N', 'Pristine']) Paints_NUsoil_M = cp.Sink('Matrix-embedded in Soil from Paints', logInflows=True, categories=['Soil_M', 'NU Soil_M', 'Matrix-embedded']) In_Air = cp.FlowCompartment('From In to Air', logInflows=True, logOutflows=True, categories=['Air']) In_Air_N = cp.Sink('Pristine in Air from In', logInflows=True, categories=['Air_N', 'Pristine']) In_Air_M = cp.Sink('Matrix-embedded in Air from In', logInflows=True, categories=['Air_M', 'Matrix-embedded']) In_WW = cp.FlowCompartment('From In to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater']) In_WW_N = cp.FlowCompartment('Pristine in Wastewater from In', logInflows=True, logOutflows=True, categories=['Wastewater_N', 'Pristine']) In_WW_M = cp.FlowCompartment('Matrix-embedded in Wastewater from In', logInflows=True, logOutflows=True, categories=['Wastewater_M', 'Matrix-embedded']) In_NUsoil = cp.FlowCompartment('From In to Soil', logInflows=True, logOutflows=True, categories=['Soil', 'NU Soil']) In_NUsoil_N = cp.Sink('Pristine in Soil from In', logInflows=True, categories=['Soil_N', 'NU Soil_N', 'Pristine']) In_NUsoil_M = cp.Sink('Matrix-embedded in Soil from In', logInflows=True, categories=['Soil_M', 'NU Soil_M', 'Matrix-embedded']) Cement_WW = cp.FlowCompartment('From Cement to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater']) Cement_WW_N = cp.FlowCompartment('Pristine from Cement to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater_N', 'Pristine']) Cement_WW_M = cp.FlowCompartment('Matrix-embedded from Cement to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater_M', 'Matrix-embedded']) Glass_WW = cp.FlowCompartment('From Glass coating to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater']) Glass_WW_N = cp.FlowCompartment('Pristine from Glass coating to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater_N', 'Pristine']) Glass_WW_M = cp.FlowCompartment('Matrix-embedded from Glass coating to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater_M', 'Matrix-embedded']) Glass_Air = cp.FlowCompartment('From Glass coating to Air', logInflows=True, logOutflows=True, categories=['Air']) Glass_Air_N = cp.Sink('Pristine from Glass coating to Air', logInflows=True, categories=['Air_N', 'Pristine']) Glass_Air_M = cp.Sink('Matrix-embedded from Glass coating to Air', logInflows=True, categories=['Air_M', 'Matrix-embedded']) Glass_NUsoil = cp.FlowCompartment('From Glass coating to Soil', logInflows=True, logOutflows=True, categories=['Soil', 'NU Soil']) Glass_NUsoil_N = cp.Sink('Pristine from Glass coating to NU Soil', logInflows=True, categories=['Soil_N', 'NU Soil_N', 'Pristine']) Glass_NUsoil_M = cp.Sink('Matrix-embedded from Glass coating to NU Soil', logInflows=True, categories=['Soil_M', 'NU Soil_M', 'Matrix-embedded']) Ceramics_WW = cp.FlowCompartment('From Ceramics coating to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater']) Ceramics_WW_N = cp.FlowCompartment('Pristine from Ceramics coating to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater_N', 'Pristine']) Ceramics_WW_M = cp.FlowCompartment('Matrix-embedded from Ceramics coating to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater_M', 'Matrix-embedded']) Ceramics_Air = cp.FlowCompartment('From Ceramics coating to Air', logInflows=True, logOutflows=True, categories=['Air']) Ceramics_Air_N = cp.Sink('Pristine from Ceramics coating to Air', logInflows=True, categories=['Air_N', 'Pristine']) Ceramics_Air_M = cp.Sink('Matrix-embedded from Ceramics coating to Air', logInflows=True, categories=['Air_M', 'Matrix-embedded']) Ceramics_NUsoil = cp.FlowCompartment('From Ceramics coating to Soil', logInflows=True, logOutflows=True, categories=['Soil', 'NU Soil']) Ceramics_NUsoil_N = cp.Sink('Pristine from Ceramics coating to NU Soil', logInflows=True, categories=['Soil_N', 'NU Soil_N', 'Pristine']) Ceramics_NUsoil_M = cp.Sink('Matrix-embedded from Ceramics coating to NU Soil', logInflows=True, categories=['Soil_M', 'NU Soil_M', 'Matrix-embedded']) RubPlas_Air = cp.FlowCompartment('From RubPlas to Air', logInflows=True, logOutflows=True, categories=['Air']) RubPlas_Air_N = cp.Sink('Pristine from RubPlas to Air', logInflows=True, categories=['Air_N', 'Pristine']) RubPlas_Air_M = cp.Sink('Matrix-embedded from RubPlas to Air', logInflows=True, categories=['Air_M', 'Matrix-embedded']) RubPlas_WW = cp.FlowCompartment('From RubPlas to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater']) RubPlas_WW_N = cp.FlowCompartment('Pristine from RubPlas to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater_N', 'Pristine']) RubPlas_WW_M = cp.FlowCompartment('Matrix-embedded from RubPlas to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater_M', 'Matrix-embedded']) Electronics_WW = cp.FlowCompartment('From Electronics to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater']) Electronics_WW_N = cp.FlowCompartment('Pristine from Electronics to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater_N', 'Pristine']) Electronics_WW_M = cp.FlowCompartment('Matrix-embedded from Electronics to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater_M', 'Matrix-embedded']) CleanAgents_WW = cp.FlowCompartment('From Cleaning agents to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater']) CleanAgents_WW_N = cp.FlowCompartment('Pristine from Cleaning agents to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater_N', 'Pristine']) Food_WW = cp.FlowCompartment('From Food to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater']) Food_WW_N = cp.FlowCompartment('Pristine from Food to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater_N', 'Pristine']) Textiles_Air = cp.FlowCompartment('From Textiles to Air', logInflows=True, logOutflows=True, categories=['Air']) Textiles_Air_N = cp.Sink('Pristine from Textiles to Air', logInflows=True, categories=['Air_N', 'Pristine']) Textiles_Air_M = cp.Sink('Matrix-embedded from Textiles to Air', logInflows=True, categories=['Air_M', 'Matrix-embedded']) Textiles_WW = cp.FlowCompartment('From Textiles to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater']) Textiles_WW_N = cp.FlowCompartment('Pristine from Textiles to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater_N', 'Pristine']) Textiles_WW_M = cp.FlowCompartment('Matrix-embedded from Textiles to Wastewater', logInflows=True, logOutflows=True, categories=['Wastewater_M', 'Matrix-embedded']) # In[12]: # Definition of WIP compartments WIP_P = cp.FlowCompartment('Product-embedded to Waste Incineration Plant', logInflows=True, logOutflows=True, categories=['WIP', 'WIP_P', 'Product-embedded']) WIP_N = cp.FlowCompartment('Pristine to Waste Incineration Plant', logInflows=True, logOutflows=True, categories=['WIP', 'WIP_N', 'Pristine']) WIP_M = cp.FlowCompartment('Matrix-embedded to Waste Incineration Plant', logInflows=True, logOutflows=True, categories=['WIP_M','Matrix-embedded']) WIP_T = cp.FlowCompartment('Transformed to Waste Incineration Plant', logInflows=True, logOutflows=True, categories=['WIP_T','Transformed']) Burning_N = cp.FlowCompartment('Pristine to Burning', logInflows=True, logOutflows=True, categories=['Burning', 'Burning_N', 'Pristine']) Burning_T = cp.FlowCompartment('Transformed to Burning', logInflows=True, logOutflows=True, categories=['Burning', 'Burning_T', 'Transformed']) Bottomash_N = cp.FlowCompartment('Pristine in Bottom ash', logInflows=True, logOutflows=True, categories=['Bottomash', 'Bottomash_N', 'Pristine']) Bottomash_T = cp.FlowCompartment('Transformed in Bottom ash', logInflows=True, logOutflows=True, categories=['Bottomash', 'Bottomash_M', 'Matrix-embedded']) Flyash_N = cp.FlowCompartment('Pristine in Flyash', logInflows=True, logOutflows=True, categories=['Flyash', 'Flyash_N', 'Pristine']) Flyash_T = cp.FlowCompartment('Transformed in Flyash', logInflows=True, logOutflows=True, categories=['Flyash', 'Flyash_T', 'Transformed']) Filterash_N = cp.FlowCompartment('Pristine in Filterash', logInflows=True, logOutflows=True, categories=['Filterash', 'Filterash_N', 'Pristine']) Filterash_T = cp.FlowCompartment('Transformed in Filterash', logInflows=True, logOutflows=True, categories=['Filterash', 'Filterash_T', 'Transformed']) WIP_Air_N = cp.Sink('From WIP to Air as Pristine', logInflows=True, categories=['Air', 'Air_N', 'Pristine']) WIP_Air_T = cp.Sink('From WIP to Air as Transformed', logInflows=True, categories=['Air', 'Air_T', 'Transformed']) # Definition of STP WW_N = cp.FlowCompartment('Pristine in Wastewater', logInflows=True, logOutflows=True, categories=['AllWastewater', 'AllWastewater_N', 'Pristine']) WW_M = cp.FlowCompartment('Matrix-embedded in Wastewater', logInflows=True, logOutflows=True, categories=['AllWastewater', 'AllWastewater_M', 'Matrix-embedded']) WW_T = cp.FlowCompartment('Transformed in Wastewater', logInflows=True, logOutflows=True, categories=['AllWastewater', 'AllWastewater_T', 'Transformed']) NoSewageSystem_N = cp.FlowCompartment('Pristine in NoSewageSystem', logInflows=True, logOutflows=True, categories=['NoSewageSystem', 'NoSewageSystem_N', 'Pristine']) NoSewageSystem_M = cp.FlowCompartment('Matrix-embedded in NoSewageSystem', logInflows=True, logOutflows=True, categories=['NoSewageSystem', 'NoSewageSystem_M', 'Matrix-embedded']) NoSewageSystem_T = cp.FlowCompartment('Transformed in NoSewageSystem', logInflows=True, logOutflows=True, categories=['NoSewageSystem', 'NoSewageSystem_T', 'Transformed']) NoSew_SW_N = cp.Sink('Pristine from no sewer to surface water', logInflows=True, categories=['Surfacewater', 'Surfacewater_N', 'Pristine']) NoSew_SW_M = cp.Sink('Matrix-embedded from no sewer to surface water', logInflows=True, categories=['Surfacewater', 'Surfacewater_M', 'Matrix-embedded']) NoSew_SW_T = cp.Sink('Transformed from no sewer to surface water', logInflows=True, categories=['Surfacewater', 'Surfacewater_T', 'Transformed']) SewageSystem_N = cp.FlowCompartment('Pristine in SewageSystem', logInflows=True, logOutflows=True, categories=['SewageSystem', 'SewageSystem_N', 'Pristine']) SewageSystem_M = cp.FlowCompartment('Matrix-embedded in SewageSystem', logInflows=True, logOutflows=True, categories=['SewageSystem', 'SewageSystem_M', 'Matrix-embedded']) SewageSystem_T = cp.FlowCompartment('Transformed in SewageSystem', logInflows=True, logOutflows=True, categories=['SewageSystem', 'SewageSystem_T', 'Transformed']) Sew_SW_N = cp.Sink('Pristine from sewer to surface water', logInflows=True, categories=['Surfacewater', 'Surfacewater_N', 'Pristine']) Sew_SW_M = cp.Sink('Matrix-embedded from sewer to surface water', logInflows=True, categories=['Surfacewater', 'Surfacewater_M', 'Matrix-embedded']) Sew_SW_T = cp.Sink('Transformed from sewer to surface water', logInflows=True, categories=['Surfacewater', 'Surfacewater_T', 'Transformed']) OnSiteTreat_N = cp.FlowCompartment('Pristine in OnsiteTreat', logInflows=True, logOutflows=True, categories = ['OnsiteTreat', 'OnsiteTreat_N', 'Pristine']) OnSiteTreat_M = cp.FlowCompartment('Matrix-embedded in OnsiteTreat', logInflows=True, logOutflows=True, categories = ['OnsiteTreat', 'OnsiteTreat_M', 'Matrix-embedded']) OnSiteTreat_T = cp.FlowCompartment('Transformed in OnsiteTreat', logInflows=True, logOutflows=True, categories = ['OnsiteTreat', 'OnsiteTreat_T', 'Transformed']) WWTP_N = cp.FlowCompartment('Pristine to WWTP', logInflows=True, logOutflows=True, categories=['WWTP', 'WWTP_N', 'Pristine']) WWTP_M = cp.FlowCompartment('Matrix-embedded to WWTP', logInflows=True, logOutflows=True, categories=['WWTP', 'WWTP_M', 'Matrix-embedded']) WWTP_T = cp.FlowCompartment('Transformed to WWTP', logInflows=True, logOutflows=True, categories=['WWTP', 'WWTP_T', 'Transformed']) TreatIOnSite_N = cp.FlowCompartment('Pristine in TreatIOnsite', logInflows=True, logOutflows=True, categories=['Treat', 'Treat_N', 'TreatIOnsite', 'TreatIOnsite_N', 'Pristine']) TreatIOnSite_M = cp.FlowCompartment('Matrix-embedded in TreatIOnSite', logInflows=True, logOutflows=True, categories=['Treat', 'Treat_M', 'TreatIOnsite', 'TreatIOnsite_M', 'Matrix-embedded']) TreatIOnSite_T = cp.FlowCompartment('Transformed in TreatIOnSite', logInflows=True, logOutflows=True, categories=['Treat', 'Treat_T', 'TreatIOnsite', 'TreatIOnsite_T', 'Transformed']) TreatI_N = cp.FlowCompartment('Pristine in TreatI', logInflows=True, logOutflows=True, categories=['STP', 'STP_N', 'Treat', 'Treat_N', 'TreatI', 'TreatI_N', 'Pristine']) TreatI_M = cp.FlowCompartment('Matrix-embedded in TreatI', logInflows=True, logOutflows=True, categories=['STP', 'STP_M', 'Treat', 'Treat_M', 'TreatI', 'TreatI_M', 'Matrix-embedded']) TreatI_T = cp.FlowCompartment('Transformed in TreatI', logInflows=True, logOutflows=True, categories=['STP', 'STP_T', 'Treat', 'Treat_T', 'TreatI', 'TreatI_T', 'Transformed']) TreatII_N = cp.FlowCompartment('Pristine in TreatII', logInflows=True, logOutflows=True, categories=['STP', 'STP_N', 'Treat', 'Treat_N', 'TreatII', 'TreatII_N', 'Pristine']) TreatII_M = cp.FlowCompartment('Matrix-embedded in TreatII', logInflows=True, logOutflows=True, categories=['STP', 'STP_M', 'Treat', 'Treat_M', 'TreatII', 'TreatII_M', 'Matrix-embedded']) TreatII_T = cp.FlowCompartment('Transformed in TreatII', logInflows=True, logOutflows=True, categories=['STP', 'STP_T', 'Treat', 'Treat_T', 'TreatII', 'TreatII_T,' 'Transformed']) TreatIII_N = cp.FlowCompartment('Pristine in TreatIII', logInflows=True, logOutflows=True, categories=['STP', 'STP_N', 'Treat', 'Treat_N', 'TreatIII', 'TreatIII_N', 'Pristine']) TreatIII_M = cp.FlowCompartment('Matrix-embedded in TreatIII', logInflows=True, logOutflows=True, categories=['STP', 'STP_M', 'Treat', 'Treat_M', 'TreatIII', 'TreatIII_M', 'Matrix-embedded']) TreatIII_T = cp.FlowCompartment('Transformed in TreatIII', logInflows=True, logOutflows=True, categories=['STP', 'STP_T', 'Treat', 'Treat_T', 'TreatIII', 'TreatIII_T', 'Transformed']) STPoverflow_N = cp.Sink('Pristine to STPoverflow and Surface water', logInflows=True, categories=['Overflow', 'Overflow_N', 'Surfacewater', 'Surfacewater_N' 'Pristine']) STPoverflow_M = cp.Sink('Matrix-embedded to STPoverflow and Surface water', logInflows=True, categories=['Overflow', 'Overflow_M', 'Surfacewater', 'Surfacewater_M', 'Matrix-embedded']) STPoverflow_T = cp.Sink('Transformed to STPoverflow and Surface water', logInflows=True, categories=['Overflow', 'Overflow_T', 'Surfacewater', 'Surfacewater_T', 'Transformed']) STPeffluent_N = cp.Sink('Pristine to STPeffluent and Surface water', logInflows=True, categories=['STPeffluent', 'STPeffluent_N', 'Surfacewater', 'Surfacewater_N', 'Pristine']) STPeffluent_M = cp.Sink('Matrix-embedded to STPeffluent and Surface water', logInflows=True, categories=['STPeffluent', 'STPeffluent_M', 'Surfacewater', 'Surfacewater_M', 'Matrix-embedded']) STPeffluent_T = cp.Sink('Transformed to STPeffluent and Surface water', logInflows=True, categories=['STPeffluent', 'STPeffluent_T', 'Surfacewater', 'Surfacewater_T', 'Transformed']) STPsludge_N = cp.FlowCompartment('Pristine in STPsludge', logInflows=True, logOutflows=True, categories=['STPSludge', 'STPSludge_N', 'Pristine']) STPsludge_M = cp.FlowCompartment('Matrix-embedded in STPsludge', logInflows=True, logOutflows=True, categories=['STPSludge', 'STPSludge_M', 'Matrix-embedded']) STPsludge_T = cp.FlowCompartment('Transformed in STPsludge', logInflows=True, logOutflows=True, categories=['STPSludge', 'STPSludge_T', 'Transformed']) # In[13]: #Definition of MMSW MMSW = cp.FlowCompartment('MMSW', logInflows=True, logOutflows=True, categories=['SolidWaste','PMCrel']) #Definition of PackW PackW = cp.FlowCompartment('PackW', logInflows=True, logOutflows=True, categories=['SolidWaste','PMCrel']) #Definition of WEEE WEEE = cp.FlowCompartment('WEEE', logInflows=True, logOutflows=True, categories=['SolidWaste','PMCrel']) #Definition of TextW TextW = cp.FlowCompartment('TextW', logInflows=True, logOutflows=True, categories=['SolidWaste','PMCrel']) #Definition of CDW CDW = cp.FlowCompartment('CDW', logInflows=True, logOutflows=True, categories=['SolidWaste', 'PMCrel']) #Definition of solar waste SolarWaste = cp.FlowCompartment('SolarWaste', logInflows=True, logOutflows=True, categories=['SolidWaste', 'PMCrel']) #Definition of PlasticW PlasticW = cp.FlowCompartment('PlasticW', logInflows=True, logOutflows=True, categories=['SolidWaste','PMCrel']) # In[14]: # Definition of Sorting_PlasticW Sorting_PlasticW = cp.FlowCompartment('Sorting_PlasticW', logInflows=True, logOutflows=True, categories=['Sorting', 'Sorting1']) # Definition of Sorting_PackW Sorting_PackW = cp.FlowCompartment('Sorting_PackW', logInflows=True, logOutflows=True, categories=['Sorting', 'Sorting1']) # Definition of Sorting_WEEE Sorting_WEEE = cp.FlowCompartment('Sorting_WEEE', logInflows=True, logOutflows=True, categories=['Sorting', 'Sorting1']) # Definition of WEEE sorting subcompartments WEEE_NotExported = cp.FlowCompartment('WEEE_NotExported', logInflows=True, logOutflows=True, categories=['Sorting', 'Sorting2']) OtherWEEE = cp.FlowCompartment('OtherWEEE', logInflows=True, logOutflows=True, categories=['Sorting', 'Sorting2']) WEEE_Plastics = cp.FlowCompartment('WEEE_Plastics', logInflows=True, logOutflows=True, categories=['Sorting', 'Sorting2']) WEEE_Plastics_OA = cp.FlowCompartment('WEEE_Plastics_OA', logInflows=True, logOutflows=True, categories=['Sorting', 'Sorting2']) WEEE_Resorting = cp.FlowCompartment('WEEE_Resorting', logInflows=True, logOutflows=True, categories=['Sorting', 'Sorting2']) # Definition of Sorting_TextW Sorting_TextW = cp.FlowCompartment('Sorting_TextW', logInflows=True, logOutflows=True, categories=['Sorting', 'Sorting1']) # Definition of TextW sorting TextW_Resorting = cp.FlowCompartment('TextW_Resorting', logInflows=True, logOutflows=True, categories = ['Sorting', 'Sorting2']) # Definition of Sorting_CDW Sorting_CDW = cp.FlowCompartment('Sorting_CDW', logInflows=True, logOutflows=True, categories=['Sorting', 'Sorting1']) #Definition of CDW sorting subcompartments Sorting_CDW_Mineral=cp.FlowCompartment('Sorting_CDW_Mineral', logInflows=True, logOutflows=True, categories=['Sorting', 'Sorting2']) Sorting_CDW_Glass=cp.FlowCompartment('Sorting_CDW_Glass', logInflows=True, logOutflows=True, categories=['Sorting', 'Sorting2']) # Disposal compartment with the same disposal rates as MMSW # - but actually it should be called 'Disposal' because it doesn't specifically needs to be disposed from sorting Sorting_Disposal = cp.FlowCompartment('Sorting_Disposal', logInflows=True, logOutflows=True, categories =['Sorting','Sorting_Disposal']) # In[15]: # Definition of Export SortWEEE_Export_P = cp.Sink('From Sorting WEEE to Export', logInflows=True, categories=['Export']) SortTextW_Export_P = cp.Sink('From Sorting Textile waste to Export', logInflows=True, categories=['Export']) SortPlasW_Export_P = cp.Sink('From Sorting Plastic waste to Export', logInflows=True, categories=['Export']) # Definition of Landfill ManufTextW_Landfill_P = cp.Sink('Product-embedded from Textile manufacturing waste to Landfill', logInflows=True, categories=['Landfill', 'Landfill_P', 'Product-embedded']) ManufPlasW_Landfill_P = cp.Sink('Product-embedded from Plastics manufacturing waste to Landfill', logInflows=True, categories=['Landfill', 'Landfill_P', 'Product-embedded']) STPSlud_Landfill_N = cp.Sink('Pristine from STP Sludge to Landfill', logInflows=True, categories=['Landfill', 'Landfill_N', 'Pristine']) STPSlud_Landfill_M = cp.Sink('Matrix-embedded from STP Sludge to Landfill', logInflows=True, categories=['Landfill', 'Landfill_M', 'Matrix-embedded']) STPSlud_Landfill_T = cp.Sink('Transformed from STP Sludge to Landfill', logInflows=True, categories=['Landfill', 'Landfill_T', 'Transformed']) MMSW_Landfill_P = cp.Sink('Product-embedded from MMSW to Landfill', logInflows=True, categories=['Landfill', 'Landfill_P', 'Product-embedded']) CDW_Landfill_P = cp.Sink('Product-embedded from CDW to Landfill', logInflows=True, categories=['Landfill', 'Landfill_P', 'Product-embedded']) SortCDWGlass_Landfill_P = cp.Sink('Product-embedded from Sorted CDW glass to Landfill', logInflows=True, categories=['Landfill', 'Landfill_P', 'Product-embedded']) SortCDWMiner_Landfill_P = cp.Sink('Product-embedded from Sorted CDW minerals to Landfill', logInflows=True, categories=['Landfill', 'Landfill_P', 'Product-embedded']) SortDisp_Landfill_P = cp.Sink('Product-embedded from other sorting to Landfill', logInflows=True, categories=['Landfill', 'Landfill_P', 'Product-embedded']) GranuPlas_Landfill_P = cp.Sink('Product-embedded from Plastic granulation to Landfill', logInflows=True, categories=['Landfill', 'Landfill_P', 'Product-embedded']) BaliText_Landfill_P = cp.Sink('Product-embedded from Textile baling to Landfill', logInflows=True, categories=['Landfill', 'Landfill_P', 'Product-embedded']) PuriSlag_Landfill_T = cp.Sink('Transformed from Purification slag to Landfill', logInflows=True, categories=['Landfill', 'Landfill_T', 'Transformed']) PuriSlag_Landfill_P = cp.Sink('Product-embedded from Purification slag to Landfill', logInflows=True, categories=['Landfill', 'Landfill_P', 'Product-embedded']) SortMiner_Landfill_P = cp.Sink('Product-embedded from Sorted minerals to Landfill', logInflows=True, categories=['Landfill', 'Landfill_P', 'Product-embedded']) Filterash_Landfill_N = cp.Sink('Pristine from Filter ash to Landfill', logInflows=True, categories=['Landfill', 'Landfill_N', 'Pristine']) Filterash_Landfill_T = cp.Sink('Transformed from Filter ash to Landfill', logInflows=True, categories=['Landfill', 'Landfill_T', 'Transformed']) Bottomash_Landfill_N = cp.Sink('Pristine from Bottom ash to Landfill', logInflows=True, categories=['Landfill', 'Landfill_N', 'Pristine']) Bottomash_Landfill_T = cp.Sink('Transformed from Bottom ash to Landfill', logInflows=True, categories=['Landfill', 'Landfill_T', 'Transformed']) # In[16]: # Definition of Sludge treated soil STsoil_N = cp.Sink('Pristine in Sludge treated soil', logInflows=True, categories=['Soil', 'Soil_N', 'ST Soil', 'ST Soil_N', 'Pristine']) STsoil_M = cp.Sink('Matrix-embedded in Sludge treated soil', logInflows=True, categories=['Soil', 'Soil_M','ST Soil', 'ST Soil_M', 'Matrix-embedded']) STsoil_T = cp.Sink('Transformed in Sludge treated soil', logInflows=True, categories=['Soil', 'Soil_T', 'ST Soil', 'ST Soil_T', 'Transformed']) #Definition of Subsurface Sew_Subsurface_N = cp.Sink('Pristine from Sewer to Subsurface', logInflows=True, categories=['Subsurface', 'Subsurface_N', 'Pristine']) Sew_Subsurface_M = cp.Sink('Matrix-embedded from Sewer to Subsurface', logInflows=True, categories=['Subsurface', 'Subsurface_M', 'Matrix-embedded']) Sew_Subsurface_T = cp.Sink('Transformed from Sewer to Subsurface', logInflows=True, categories=['Subsurface', 'Subsurface_T', 'Transformed']) OnSiteTreat_Subsurface_N = cp.Sink('Pristine from On-site treatment to Subsurface', logInflows=True, categories=['Subsurface', 'Subsurface_N', 'Pristine']) OnSiteTreat_Subsurface_M = cp.Sink('Matrix-embedded from On-site treatment to Subsurface', logInflows=True, categories=['Subsurface', 'Subsurface_M', 'Matrix-embedded']) OnSiteTreat_Subsurface_T = cp.Sink('Transformed from On-site treatment to Subsurface', logInflows=True, categories=['Subsurface', 'Subsurface_T', 'Transformed']) #Definition of the sludge that stays on site OnSiteSludge_N = cp.Sink('Pristine in OnsiteSludge', logInflows=True, categories=['OnsiteSludge', 'OnsiteSludge_N', 'Pristine']) OnSiteSludge_M = cp.Sink('Matrix-embedded in OnsiteSludge', logInflows=True, categories=['OnsiteSludge', 'OnsiteSludge_M', 'Matrix-embedded']) OnSiteSludge_T = cp.Sink('Transformed in OnsiteSludge', logInflows=True, categories=['OnsiteSludge', 'OnsiteSludge_T', 'Transformed']) # In[17]: ''' -----------------Define necessary commpartments for recycling systems--------------------------''' # Definition of Reprocessing #define the waste categories in reprocesseing systems as product-embedded forms Reprocess_WEEE_NotExW_P=cp.FlowCompartment('Light bulbs waste to reprocess as product-embedded forms',logInflows=True, logOutflows=True, categories=['Reprocessing', 'Reprocess_P']) Reprocess_Plas_P=cp.FlowCompartment('Plastics waste to reprocess as product-embedded forms',logInflows=True, logOutflows=True, categories=['Reprocessing','Reprocess_P','Reprocessed Plastics']) Reprocess_Glass_P=cp.FlowCompartment('Glass waste to reprocess as product-embedded forms',logInflows=True, logOutflows=True, categories=['Reprocessing','Reprocess_P','Reprocessed Glass']) Reprocess_Metal_P=cp.FlowCompartment('Metal waste to reprocess as product-embedded forms',logInflows=True, logOutflows=True, categories=['Reprocessing','Reprocess_P','Reprocessed Metals']) Reprocess_Mineral_P=cp.FlowCompartment('Mineral waste to reprocess as product-embedded forms',logInflows=True, logOutflows=True, categories=['Reprocessing','Reprocess_P','Reprocessed Minerals']) Reprocess_Text_P=cp.FlowCompartment('Textile waste to reprocess as product-embedded forms',logInflows=True, logOutflows=True, categories=['Reprocessing','Reprocess_P','Reprocessed Textiles']) Reprocess_Solar_P=cp.Sink('Solar waste to reprocess as product-embedded forms',logInflows=True, categories=['Reprocessing','Reprocess_P','Reprocessed Solar waste']) # define plastic recycling system Washing_Plas=cp.FlowCompartment('Washing of plastics',logInflows=True, logOutflows=True, categories=['ReprWashPlas']) Melting_Plas=cp.FlowCompartment('Melting of plastics',logInflows=True, logOutflows=True, categories=['ReprMeltPlas']) Granulation_Plas=cp.FlowCompartment('Granulation of plastics',logInflows=True, logOutflows=True, categories=['ReprGranPlas']) Air_Plas= cp.FlowCompartment('ENM particles in Air from shredding plastics', logInflows=True, logOutflows=True,categories=['Air']) # define textile recycling system Washing_Text=cp.FlowCompartment('Washing of textiles',logInflows=True, logOutflows=True, categories=['ReprWashText']) Baling_Text=cp.FlowCompartment('ENM particles in baling as product-embedded forms with textiles',logInflows=True, logOutflows=True, categories=['ReprWashBali']) Air_Text= cp.FlowCompartment('Textile particles in Air from shredding textiles', logInflows=True,logOutflows=True, categories=['Air']) WW_Text = cp.FlowCompartment('ENM particles in Wastewater from washing textiles',logInflows=True,logOutflows=True, categories=['Wastewater']) # define metal recycling system Melting_Metal=cp.FlowCompartment('Melting of metals',logInflows=True, logOutflows=True, categories=['ReprMeltMeta']) Purification_Metal_P=cp.FlowCompartment('ENM as product-embedded forms in flows to purifaction of metals',logInflows=True, logOutflows=True, categories=['ReprPuriMeta', 'ReprPuriMeta_P']) Purification_Metal_T=cp.FlowCompartment('ENM as transformed forms in flows to purifaction of metals',logInflows=True, logOutflows=True, categories=['ReprPuriMeta', 'ReprPuriMeta_T']) #Moulding_Metal_T=cp.FlowCompartment('ENM as transformed forms in flows to moulding of metals',logInflows=True, logOutflows=True, categories=['Procedures in reprocessing system']) #Moulding_Metal_P=cp.FlowCompartment('ENM as product-embedded forms in flows to moulding of metals',logInflows=True, logOutflows=True, categories=['Procedures in reprocessing system']) Puri_slag_P=cp.FlowCompartment('ENM as product-embedded forms in flows to slags of metals',logInflows=True, logOutflows=True, categories=['ReprPuriSlag', 'ReprPuriSlag_P']) Puri_slag_T=cp.FlowCompartment('ENM as transformed forms in flows to slags of metals',logInflows=True, logOutflows=True, categories=['ReprPuriSlag', 'ReprPuriSlag_T']) Air_Metal= cp.FlowCompartment('ENM particles in Air from shredding metals', logInflows=True, logOutflows=True,categories=['Air']) # define glass recycling system Melting_Glass=cp.FlowCompartment('Melting of glass',logInflows=True, logOutflows=True, categories=['ReprMeltGlas']) #Moulding_Glass_T=cp.FlowCompartment('Moulding of glass',logInflows=True, logOutflows=True, categories=['Procedures in reprocessing system']) Air_Glass= cp.FlowCompartment('ENM particles in Air from crushing glass', logInflows=True,logOutflows=True, categories=['Air']) #Landfill_Glass_P = cp.Sink('ENM partciles in Landfill as product-embedded forms with glass', logInflows=True, categories=['Landfill']) # define mineral recycling system Sorting_Mineral=cp.FlowCompartment('Sorting/Sieving of mineral',logInflows=True, logOutflows=True, categories=['ReprSortMine']) Air_Mineral= cp.FlowCompartment('ENM particles in Air from crushing minerals', logInflows=True,logOutflows=True, categories=['Air']) # define air NM Air_NM_Plas = cp.FlowCompartment('N and M particles in air from plastics',logInflows=True,logOutflows=True, categories=['Air_NM']) Air_NM = cp.FlowCompartment('N and M particles in air from metal/textile/mineral/glass',logInflows=True,logOutflows=True, categories=['Air_NM']) # defined RPAir here in order to sum up flows to air from reprocessing RPAir_N= cp.Sink('Pristine particles in air from reprocessing systems',logInflows=True, categories=['Air_N', 'Pristine']) RPAir_M= cp.Sink('Matrix-embedded particles in air from reprocessing systems',logInflows=True, categories=['Air_M','Matrix-embedded']) RPWIP_P = cp.FlowCompartment('Product-embedded to Waste Incineration Plant from reprocessing', logInflows=True, logOutflows=True, categories=['WIP','WIP_P', 'Product_embedded']) RPWW_N = cp.FlowCompartment('Pristine in Wastewater from reprocessing', logInflows=True, logOutflows=True, categories=['Wastewater', 'Wastewater_N', 'Pristine']) RPWW_M = cp.FlowCompartment('Matrix-embedded in Wastewater from reprocessing', logInflows=True, logOutflows=True, categories=['Wastewater', 'Wastewater_M', 'Matrix-embedded']) # define reuse compartment to store all final products from recycling factories GranuPlas_Reuse_P = cp.Sink('Product-embedded from Plastic granulation to Reuse', logInflows=True, categories=['Reuse', 'Reuse_P', 'Product-embedded']) BaliText_Reuse_P = cp.Sink('Product-embedded from Textile baling to Reuse', logInflows=True, categories=['Reuse', 'Reuse_P', 'Product-embedded']) PuriMetal_Reuse_P = cp.Sink('Product-embedded from Metal purification to Reuse', logInflows=True, categories=['Reuse', 'Reuse_P', 'Product-embedded']) PuriMetal_Reuse_T = cp.Sink('Transformed from Metal purification to Reuse', logInflows=True, categories=['Reuse', 'Reuse_T', 'Transformed']) PuriSlag_Reuse_P = cp.Sink('Product-embedded from Purification slag to Reuse', logInflows=True, categories=['Reuse', 'Reuse_P', 'Product-embedded']) PuriSlag_Reuse_T = cp.Sink('Transformed from Purification slag to Reuse', logInflows=True, categories=['Reuse', 'Reuse_T', 'Transformed']) MeltGlass_Reuse_T = cp.Sink('Transformed from Glass melting to Reuse', logInflows=True, categories=['Reuse', 'Reuse_T', 'Transformed']) SortMiner_Reuse_P = cp.Sink('Product-embedded from Sorted minerals to Reuse', logInflows=True, categories=['Reuse', 'Reuse_P', 'Product-embedded']) Filterash_Reuse_N = cp.Sink('Pristine from Filter ash to Reuse', logInflows=True, categories=['Reuse', 'Reuse_N', 'Pristine']) Filterash_Reuse_T = cp.Sink('Transformed from Filter ash to Reuse', logInflows=True, categories=['Reuse', 'Reuse_T', 'Transformed']) Bottomash_Reuse_N = cp.Sink('Pristine from Bottom ash to Reuse', logInflows=True, categories=['Reuse', 'Reuse_N', 'Pristine']) Bottomash_Reuse_T = cp.Sink('Transformed from Bottom ash to Reuse', logInflows=True, categories=['Reuse', 'Reuse_T', 'Transformed']) # In[19]: compartmentList=[totalInflow, Production, Prod_Air_N, Manuf_Air_N, Manufacture, Manuf_Plas, Manuf_Text, manuf_textiles_waste, manuf_plastics_waste, ManufTextW_Reuse_P, ManufPlasW_Reuse_P, Consumption, PersCare, PersCare_Use, PersCare_EoL, Paints, Paints_Use, Paints_EoL, In, In_Use, In_EoL, Cement, Cement_Use, Cement_EoL, Glass, Glass_Use, Glass_EoL, Ceramics, Ceramics_Use, Ceramics_EoL, RubPlas, RubPlas_Use, RubPlas_EoL, Solar, Solar_EoL, Electronics, Electronics_Use, Electronics_EoL, CleanAgents, CleanAgents_Use, CleanAgents_EoL, Food, Food_Use, Food_EoL, Textiles, Textiles_Use, Textiles_EoL, PersCare_WW, PersCare_WW_N, PersCare_SW, PersCare_SW_N, Paints_Air, Paints_Air_N, Paints_Air_M, Paints_WW, Paints_WW_N, Paints_WW_M, Paints_NUsoil, Paints_NUsoil_N, Paints_NUsoil_M, In_Air, In_Air_N, In_Air_M, In_WW, In_WW_N, In_WW_M, In_NUsoil, In_NUsoil_N, In_NUsoil_M, Cement_WW, Cement_WW_N, Cement_WW_M, Glass_WW, Glass_WW_N, Glass_WW_M, Glass_Air, Glass_Air_N, Glass_Air_M, Glass_NUsoil, Glass_NUsoil_N, Glass_NUsoil_M, Ceramics_WW, Ceramics_WW_N, Ceramics_WW_M, Ceramics_Air, Ceramics_Air_N, Ceramics_Air_M, Ceramics_NUsoil, Ceramics_NUsoil_N, Ceramics_NUsoil_M, RubPlas_Air, RubPlas_Air_N, RubPlas_Air_M, RubPlas_WW, RubPlas_WW_N, RubPlas_WW_M, Electronics_WW, Electronics_WW_N, Electronics_WW_M, CleanAgents_WW, CleanAgents_WW_N, Food_WW, Food_WW_N, Textiles_WW, Textiles_WW_N, Textiles_WW_M, Textiles_Air, Textiles_Air_N, Textiles_Air_M, WW_N, WW_M, WW_T, OnSiteTreat_N, OnSiteTreat_M, OnSiteTreat_T, SewageSystem_N, SewageSystem_M, SewageSystem_T, NoSewageSystem_N, NoSewageSystem_M, NoSewageSystem_T, TreatIOnSite_N, TreatIOnSite_M, TreatIOnSite_T, WWTP_N, WWTP_M, WWTP_T, STPoverflow_N, STPoverflow_M, STPoverflow_T, STPeffluent_N, STPeffluent_M, STPeffluent_T, STPsludge_N, STPsludge_M, STPsludge_T, OnSiteSludge_N, OnSiteSludge_M, OnSiteSludge_T, TreatI_N, TreatII_N, TreatIII_N, TreatI_M, TreatII_M, TreatIII_M, TreatI_T, TreatII_T, TreatIII_T, NoSew_SW_N, NoSew_SW_T, NoSew_SW_M, Sew_SW_N, Sew_SW_T, Sew_SW_M, MMSW, PackW, WEEE, TextW, CDW, SolarWaste, PlasticW, Sorting_PackW, Sorting_WEEE, Sorting_TextW, Sorting_CDW, Sorting_PlasticW, WEEE_NotExported, OtherWEEE, WEEE_Plastics, WEEE_Plastics_OA, WEEE_Resorting, TextW_Resorting, Sorting_CDW_Mineral, Sorting_CDW_Glass, Sorting_Disposal, GranuPlas_Reuse_P, BaliText_Reuse_P, PuriMetal_Reuse_P, PuriMetal_Reuse_T, PuriSlag_Reuse_P, PuriSlag_Reuse_T, MeltGlass_Reuse_T, SortMiner_Reuse_P, Filterash_Reuse_N, Filterash_Reuse_T, Bottomash_Reuse_N, Bottomash_Reuse_T, Air_NM, RPAir_N, RPAir_M, RPWIP_P, RPWW_N, RPWW_M, Washing_Plas, Melting_Plas, Granulation_Plas, Air_Plas, Air_NM_Plas, Washing_Text, Baling_Text, Air_Text, WW_Text, Melting_Metal, Purification_Metal_P, Purification_Metal_T, Puri_slag_P, Puri_slag_T,#Moulding_Metal_P,Moulding_Metal_T, Air_Metal, Melting_Glass, Air_Glass, Sorting_Mineral, Air_Mineral, Reprocess_WEEE_NotExW_P, Reprocess_Plas_P, Reprocess_Text_P, Reprocess_Glass_P, Reprocess_Metal_P, Reprocess_Mineral_P, Reprocess_Solar_P, WIP_N, WIP_M, WIP_P, WIP_T, Burning_N, Burning_T, Bottomash_N, Bottomash_T, Flyash_N, Flyash_T, Filterash_N, Filterash_T, WIP_Air_N, WIP_Air_T, SortWEEE_Export_P, SortTextW_Export_P, SortPlasW_Export_P, ManufTextW_Landfill_P, ManufPlasW_Landfill_P, STPSlud_Landfill_N, STPSlud_Landfill_M, STPSlud_Landfill_T, MMSW_Landfill_P, CDW_Landfill_P, SortCDWGlass_Landfill_P, SortCDWMiner_Landfill_P, SortDisp_Landfill_P, GranuPlas_Landfill_P, BaliText_Landfill_P, PuriSlag_Landfill_T, PuriSlag_Landfill_P, SortMiner_Landfill_P, Filterash_Landfill_N, Filterash_Landfill_T, Bottomash_Landfill_N, Bottomash_Landfill_T, STsoil_N, STsoil_M, STsoil_T, Sew_Subsurface_N, Sew_Subsurface_M, Sew_Subsurface_T, OnSiteTreat_Subsurface_N, OnSiteTreat_Subsurface_M, OnSiteTreat_Subsurface_T] model.setCompartments(compartmentList) # In[20]: Productiondata = pd.read_csv('nano_TiO2production2000-2016.csv', header=None) Productiondata.shape # In[21]: # Production for the whole Europe: Productiondata_EU = Productiondata.values #in earlier code as " Productiondata = Productiondata.as_matrix(Productiondata) # Scaling production for the UK only: Productiondata_anatase = Productiondata_EU*0.062*0.096 # Putting the numbers in the model: Productionvolume_EU = [] Productionvolume_anatase = [] for i in np.arange(0,17): Productionvolume_EU.append(Productiondata_EU[:,i]) Productionvolume_anatase.append(Productiondata_anatase[:,i]) #model.addInflow(cp.ExternalListInflow(totalInflow, inflowList=Productionvolume)) periodRange = np.arange(0,17) # Setting inflow to production model.addInflow(cp.ExternalListInflow(totalInflow, [cp.RandomChoiceInflow(Productionvolume_anatase[x]) for x in periodRange])) # Adding imports. (country-specific) #SF_Prod = 0 #SF_Manuf = 0 #SF_Cons = 0 #model.addInflow(cp.ExternalListInflow(Import_Manuf, [cp.RandomChoiceInflow(Productionvolume_anatase[x]) for x in periodRange])) #model.addInflow(cp.ExternalListInflow(Import_Cons, [cp.RandomChoiceInflow(Productionvolume_anatase[x]) for x in periodRange])) # In[22]: totalInflow.transfers = [cp.ConstTransfer(1, Production)] #totalInflow.transfers = [cp.ConstTransfer(1, Manufacture)] #totalInflow.transfers = [cp.ConstTransfer(1, Consumption)] # In[23]: ## Production and Manufacture Production.localRelease = cp.ListRelease([1.0]) Manufacture.localRelease = cp.ListRelease([1.0]) Production.transfers = [cp.StochasticTransfer(nr.triangular, [0.0002, 0.0125, 0.0467], WW_N, priority=2), cp.StochasticTransfer(nr.triangular, [0.0000008, 0.0000016, 0.0000024], Prod_Air_N, priority=2), cp.ConstTransfer(1, Manufacture, priority=1)] Manufacture.transfers = [cp.StochasticTransfer(nr.triangular, [0, 0.968, 1], Manuf_Plas, priority=2), cp.StochasticTransfer(nr.triangular, [0, 0.032, 0.04], Manuf_Text, priority=2) ] Manuf_Text.transfers = [cp.StochasticTransfer(nr.triangular, [0.00000085, 0.0000017, 0.00000255], Manuf_Air_N, priority=2), cp.StochasticTransfer(nr.uniform, [0.001, 0.0016], WW_N, priority=2), cp.StochasticTransfer(nr.uniform, [0, 0.018], manuf_textiles_waste, priority=2), cp.ConstTransfer(1, Consumption, priority=1)] Manuf_Plas.transfers = [cp.StochasticTransfer(nr.triangular, [0.00000085, 0.0000017, 0.00000255], Manuf_Air_N, priority=2), cp.StochasticTransfer(nr.uniform, [0.001, 0.0016], WW_N, priority=2), cp.StochasticTransfer(nr.uniform, [0.03, 0.08], manuf_plastics_waste, priority=2), cp.ConstTransfer(1, Consumption, priority=1)] manuf_textiles_waste.transfers = [cp.StochasticTransfer(nr.triangular, [0.25, 0.5, 0.75], ManufTextW_Reuse_P), cp.StochasticTransfer(nr.triangular, [0.3*0.25, 0.3*0.5, 0.3*0.75], WIP_P), cp.StochasticTransfer(nr.triangular, [0.7*0.25, 0.7*0.5, 0.7*0.75], ManufTextW_Landfill_P)] manuf_plastics_waste.transfers = [cp.RandomChoiceTransfer(tf.TriangTruncDet(0.455, 0.91, 1.365, 10000, 0, 1), ManufPlasW_Reuse_P), cp.StochasticTransfer(nr.triangular, [0.045, 0.09, 0.135], WIP_P), cp.StochasticTransfer(nr.triangular, [0.001, 0.002, 0.003], ManufPlasW_Landfill_P)] # In[24]: # ENM allocation to product categories Consumption.transfers = [cp.ConstTransfer(0, PersCare), cp.StochasticTransfer(nr.triangular,[0, 0.968, 1], RubPlas), cp.ConstTransfer(0, Electronics), cp.ConstTransfer(0, CleanAgents), cp.ConstTransfer(0, Food), cp.StochasticTransfer(nr.triangular,[0, 0.032, 0.04], Textiles), cp.ConstTransfer(0, Paints), cp.ConstTransfer(0, In), cp.ConstTransfer(0, Cement), cp.ConstTransfer(0, Glass), cp.ConstTransfer(0, Ceramics), cp.ConstTransfer(0, Solar)] # In[25]: # Products categories to their USE and EoL compartments PersCare.transfers = [cp.ConstTransfer(1, PersCare_Use, priority=1), cp.StochasticTransfer(nr.triangular, [0.025, 0.05, 0.075], PersCare_EoL, priority=2)] Paints.transfers = [cp.ConstTransfer(1, Paints_EoL, priority=1), cp.StochasticTransfer(nr.triangular, [0.005, 0.01, 0.015], Paints_Use, priority=2)] In.transfers = [cp.ConstTransfer(1, In_EoL, priority=1), cp.StochasticTransfer(tf.TriangTrunc, [0.001, 3.22, 1, 0, 1], In_Use, priority=2)] Cement.transfers = [cp.ConstTransfer(1, Cement_EoL, priority=1), cp.StochasticTransfer(nr.triangular, [0.005, 0.01, 0.015], Cement_Use, priority=2)] Glass.transfers = [cp.ConstTransfer(1, Glass_EoL, priority=1), cp.StochasticTransfer(nr.triangular, [0.175, 0.35, 0.525], Glass_Use, priority=2)] Ceramics.transfers = [cp.ConstTransfer(1, Ceramics_EoL, priority=1), cp.StochasticTransfer(nr.triangular, [0.175, 0.35, 0.525], Ceramics_Use, priority=2)] RubPlas.transfers = [cp.ConstTransfer(1, RubPlas_EoL, priority=1), cp.StochasticTransfer(nr.triangular, [0.015, 0.035, 0.045], RubPlas_Use, priority=2)] Solar.transfers = [cp.ConstTransfer(1, Solar_EoL, priority=1)] Electronics.transfers = [cp.ConstTransfer(1, Electronics_EoL, priority=1), cp.StochasticTransfer(nr.triangular, [0.15, 0.30, 0.45], Electronics_Use, priority=2)] CleanAgents.transfers = [cp.ConstTransfer(1, CleanAgents_Use, priority=1), cp.StochasticTransfer(nr.triangular, [0.025, 0.05, 0.075], CleanAgents_EoL, priority=2)] Food.transfers = [cp.ConstTransfer(1, Food_Use, priority=1), cp.StochasticTransfer(nr.triangular, [0.05, 0.10, 0.15], Food_EoL, priority=2)] # Determination of releases during use: # previous from Sun et al = 0.03 (0.015 the first year), # new from LEITAT = 0.093 after ten washings (i.e. assumed during the first year), # so in total (wastewater + air) = 0.186 released during use. Textiles_Use_data = np.concatenate([nr.triangular(0.015, 0.03, 0.045, 5000), nr.triangular(0.093, 0.186, 0.279, 5000)]) Textiles.transfers = [cp.ConstTransfer(1, Textiles_EoL, priority=1), cp.RandomChoiceTransfer(Textiles_Use_data, Textiles_Use, priority=2)] # In[26]: # Personal care PersCare_Use.localRelease = cp.ListRelease([0.90,0.10]) PersCare_EoL.localRelease = cp.ListRelease([0.90,0.10]) PersCare_Use.transfers = [cp.ConstTransfer(1, PersCare_WW, priority=1), cp.StochasticTransfer(nr.triangular, [0.05,0.10,0.15], PersCare_SW, priority=2)] PersCare_WW.transfers = [cp.ConstTransfer(1, PersCare_WW_N, priority=1)] PersCare_SW.transfers = [cp.ConstTransfer(1, PersCare_SW_N, priority=1)] PersCare_EoL.transfers = [cp.ConstTransfer(1, PackW, priority=1)] # In[27]: # Paints Paints_Use.localRelease = cp.ListRelease([0.9,0.01666667,0.01666667,0.01666667,0.01666667, 0.01666667,0.01666667]) Paints_EoL.localRelease = cp.ListRelease([0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.001722281, 0.001037804,0.001572363,0.00232936,0.003374171,0.004779078, 0.006618624,0.00896268,0.01186737,0.01536449,0.01945038, 0.02407602,0.0291399,0.03448561,0.03990567,0.04515209, 0.04995371,0.05403864,0.0571594,0.05911772,0.05978529, 0.05911772,0.0571594,0.05403864,0.04995371,0.04515209, 0.03990567,0.03448561,0.0291399,0.02407602,0.01945038, 0.01536449,0.01186737,0.00896268,0.006618624,0.004779078, 0.003374171,0.00232936,0.001572363,0.001037804,0.001722281]) Paints_Use.transfers = [cp.StochasticTransfer(nr.triangular, [0.25, 0.50, 0.75], Paints_WW, priority=1), cp.StochasticTransfer(nr.triangular, [0.125, 0.25,0.375], Paints_Air, priority=1), cp.StochasticTransfer(nr.triangular, [0.125, 0.25,0.375], Paints_NUsoil, priority=1)] Paints_WW.transfers = [cp.StochasticTransfer(nr.triangular, [0.15, 0.30, 0.45], Paints_WW_N, priority=1), cp.StochasticTransfer(tf.TriangTruncDet, [0.35, 0.7, 1.05, 1, 0, 1], Paints_WW_M, priority=1)] Paints_Air.transfers = [cp.StochasticTransfer(nr.uniform, [1e-6, 1], Paints_Air_N, priority=1), cp.StochasticTransfer(nr.uniform, [1e-6, 1], Paints_Air_M, priority=1)] Paints_NUsoil.transfers = [cp.StochasticTransfer(nr.triangular, [0.15, 0.30, 0.45], Paints_NUsoil_N, priority=1), cp.StochasticTransfer(tf.TriangTruncDet, [0.35, 0.70, 1.05, 1, 0, 1], Paints_NUsoil_M, priority=1)] Paints_EoL.transfers = [cp.StochasticTransfer(tf.TriangTruncDet, [0.475, 0.95, 1.9, 1, 0, 1], CDW, priority=1), cp.StochasticTransfer(nr.triangular, [0.025, 0.05, 0.075], MMSW, priority=1)] # In[27]: # In In_Use.localRelease = cp.FixedRateRelease(0.1) #Or cp.ListRelease([0.1, 0.1, 0,1, 0.1, 0.1, 0,1, 0.1, 0.1, 0,1, 0.1]) In_EoL.localRelease = cp.ListRelease([0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.001722281, 0.001037804,0.001572363,0.00232936,0.003374171,0.004779078, 0.006618624,0.00896268,0.01186737,0.01536449,0.01945038, 0.02407602,0.0291399,0.03448561,0.03990567,0.04515209, 0.04995371,0.05403864,0.0571594,0.05911772,0.05978529, 0.05911772,0.0571594,0.05403864,0.04995371,0.04515209, 0.03990567,0.03448561,0.0291399,0.02407602,0.01945038, 0.01536449,0.01186737,0.00896268,0.006618624,0.004779078, 0.003374171,0.00232936,0.001572363,0.001037804,0.001722281]) In_Use.transfers = [cp.StochasticTransfer(nr.triangular, [0.09, 0.11, 0.136], In_WW, priority=1), cp.StochasticTransfer(nr.triangular, [0.125, 0.25,0.375], In_Air, priority=1), cp.StochasticTransfer(nr.triangular, [0.512, 0.64,0.768], In_NUsoil, priority=1)] In_WW.transfers = [cp.StochasticTransfer(nr.triangular, [0.15, 0.30, 0.45], In_WW_N, priority=1), cp.StochasticTransfer(tf.TriangTruncDet, [0.35, 0.7, 1.05, 1, 0, 1], In_WW_M, priority=1)] In_Air.transfers = [cp.StochasticTransfer(nr.uniform, [1e-6, 1], In_Air_N, priority=1), cp.StochasticTransfer(nr.uniform, [1e-6, 1], In_Air_M, priority=1)] In_NUsoil.transfers = [cp.StochasticTransfer(nr.triangular, [0.15, 0.30, 0.45], In_NUsoil_N, priority=1), cp.StochasticTransfer(tf.TriangTruncDet, [0.35, 0.70, 1.05, 1, 0, 1], In_NUsoil_M, priority=1)] In_EoL.transfers = [cp.StochasticTransfer(tf.TriangTruncDet, [0.475, 0.95, 1.9, 1, 0, 1], CDW, priority=1), cp.StochasticTransfer(nr.triangular, [0.025, 0.05, 0.075], MMSW, priority=1)] # In[28]: # Cement Cement_Use.localRelease = cp.ListRelease([0.9,0.01666667,0.01666667,0.01666667,0.01666667, 0.01666667,0.01666667,0.01666667]) Cement_EoL.localRelease = cp.ListRelease([0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.001722281, 0.001037804,0.001572363,0.00232936,0.003374171,0.004779078, 0.006618624,0.00896268,0.01186737,0.01536449,0.01945038, 0.02407602,0.0291399,0.03448561,0.03990567,0.04515209, 0.04995371,0.05403864,0.0571594,0.05911772,0.05978529, 0.05911772,0.0571594,0.05403864,0.04995371,0.04515209, 0.03990567,0.03448561,0.0291399,0.02407602,0.01945038, 0.01536449,0.01186737,0.00896268,0.006618624,0.004779078, 0.003374171,0.00232936,0.001572363,0.001037804,0.001722281]) Cement_Use.transfers = [cp.ConstTransfer(1, Cement_WW, priority=1)] Cement_WW.transfers = [cp.StochasticTransfer(nr.uniform, [1e-6, 0.01], Cement_WW_N, priority=1), cp.StochasticTransfer(nr.uniform, [0.99, 1], Cement_WW_M, priority=1)] Cement_EoL.transfers = [cp.ConstTransfer(1, CDW, priority=1)] # In[29]: # Glass coatings Glass_Use.localRelease = cp.ListRelease([0.9,0.01111111,0.01111111,0.01111111,0.01111111, 0.01111111,0.01111111,0.01111111,0.01111111,0.01111111]) Glass_EoL.localRelease = cp.ListRelease([0,0,0,0,0.003466974, 0.01439745,0.04894278,0.1172529,0.1980285,0.2358228, 0.1980285,0.1172529,0.04894278,0.01439745,0.003466974]) Glass_Use.transfers = [cp.ConstTransfer(1, Glass_WW, priority=1), cp.StochasticTransfer(nr.triangular, [0.05,0.10,0.15], Glass_Air, priority=2), cp.StochasticTransfer(nr.triangular, [0.05,0.10,0.15], Glass_NUsoil, priority=2)] Glass_WW.transfers = [cp.StochasticTransfer(nr.triangular, [0.1, 0.3, 0.5], Glass_WW_N, priority=1), cp.StochasticTransfer(nr.triangular, [0.5, 0.7, 0.9], Glass_WW_M, priority=1)] Glass_Air.transfers = [cp.StochasticTransfer(nr.uniform, [1e-6, 1], Glass_Air_N, priority=1), cp.StochasticTransfer(nr.uniform, [1e-6, 1], Glass_Air_M, priority=1)] Glass_NUsoil.transfers = [cp.StochasticTransfer(nr.triangular, [0.1, 0.3, 0.5], Glass_NUsoil_N, priority=1), cp.StochasticTransfer(nr.triangular, [0.5, 0.7, 0.9], Glass_NUsoil_M, priority=1)] Glass_EoL.transfers = [cp.ConstTransfer(1, CDW, priority=1)] # In[30]: # Ceramic coatings Ceramics_Use.localRelease = cp.ListRelease([0.9,0.01111111,0.01111111,0.01111111,0.01111111, 0.01111111,0.01111111,0.01111111,0.01111111,0.01111111]) Ceramics_EoL.localRelease = cp.ListRelease([0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.001722281, 0.001037804,0.001572363,0.00232936,0.003374171,0.004779078, 0.006618624,0.00896268,0.01186737,0.01536449,0.01945038, 0.02407602,0.0291399,0.03448561,0.03990567,0.04515209, 0.04995371,0.05403864,0.0571594,0.05911772,0.05978529, 0.05911772,0.0571594,0.05403864,0.04995371,0.04515209, 0.03990567,0.03448561,0.0291399,0.02407602,0.01945038, 0.01536449,0.01186737,0.00896268,0.006618624,0.004779078, 0.003374171,0.00232936,0.001572363,0.001037804,0.001722281]) Ceramics_Use.transfers = [cp.ConstTransfer(1, Ceramics_WW, priority=1), cp.StochasticTransfer(nr.triangular, [0.05,0.10,0.15], Ceramics_Air, priority=2), cp.StochasticTransfer(nr.triangular, [0.05,0.10,0.15], Ceramics_NUsoil, priority=2)] Ceramics_WW.transfers = [cp.StochasticTransfer(nr.triangular, [0.1, 0.3, 0.5], Ceramics_WW_N, priority=1), cp.StochasticTransfer(nr.triangular, [0.5, 0.7, 0.9], Ceramics_WW_M, priority=1)] Ceramics_Air.transfers = [cp.StochasticTransfer(nr.uniform, [1e-6, 1], Ceramics_Air_N, priority=1), cp.StochasticTransfer(nr.uniform, [1e-6, 1], Ceramics_Air_M, priority=1)] Ceramics_NUsoil.transfers = [cp.StochasticTransfer(nr.triangular, [0.1, 0.3, 0.5], Ceramics_NUsoil_N, priority=1), cp.StochasticTransfer(nr.triangular, [0.5, 0.7, 0.9], Ceramics_NUsoil_M, priority=1)] Ceramics_EoL.transfers = [cp.ConstTransfer(1, CDW, priority=1)] # In[31]: # Rubber, ceramic and plastic additives #Plastics_Use.localRelease = cp.FixedRateRelease(0.125) #Sport_Use.localRelease = cp.FixedRateRelease(0.1428571) #Mean = 0.134, so: RubPlas_Use.localRelease = cp.FixedRateRelease(0.134) #Plastics_EoL.localRelease = cp.ListRelease([0, 0, 0.003466974, 0.01439745, 0.04894278, # 0.1172529, 0.1980285, 0.2358228, 0.1980285, 0.1172529, # 0.04894278,0.01439745,0.003466974]) #Sport_EoL.localRelease = cp.ListRelease([0, 0, 0, 0.006209665, 0.06059754, # 0.2417303, 0.3829249, 0.2417303, 0.06059754, 0.006209665]) #Mean = [0, 0, 0.001733487, 0.0103035575, 0.05477016, # 0.1794916, 0.2904767, 0.23877655, 0.12931302, 0.0617312825, # 0.02447139, 0.007198725, 0.001733487], #adds to 1, so: RubPlas_EoL.localRelease = cp.ListRelease([0, 0, 0.001733487, 0.0103035575, 0.05477016, 0.1794916, 0.2904767, 0.23877655, 0.12931302, 0.0617312825, 0.02447139, 0.007198725, 0.001733487]) RubPlas_Use.transfers = [cp.StochasticTransfer(nr.triangular, [0.075, 0.15, 0.225], RubPlas_Air, priority=1), cp.StochasticTransfer(tf.TriangTruncDet, [0.425, 0.85, 1.275, 1, 0, 1], RubPlas_WW, priority=1)] RubPlas_Air.transfers = [cp.StochasticTransfer(nr.uniform, [1e-6, 1], RubPlas_Air_N, priority=1), cp.StochasticTransfer(nr.uniform, [1e-6, 1], RubPlas_Air_M, priority=1)] TransfSportWWN = nr.uniform(1e-6, 1, 10000) TransfPlasWWN = nr.triangular(0.15, 0.30, 0.45, 10000) TransfRubPlasWWN = (TransfSportWWN + TransfPlasWWN)/2 TransfSportWWM = nr.uniform(1e-6, 1, 10000) TransfPlasWWM = tf.TriangTruncDet(0.35, 0.70, 1.05, 10000, 0, 1) TransfRubPlasWWM = (TransfSportWWM + TransfPlasWWM)/2 RubPlas_WW.transfers = [cp.TimeDependendListTransfer([nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN), nr.choice(TransfRubPlasWWN)], RubPlas_WW_N, RubPlas_WW, priority=2), cp.ConstTransfer(1, RubPlas_WW_M, priority=1)] # Data sources: Kawecki et al. 2018 # Take the average of 7 types of polymers RubPlas_EoL.transfers = [cp.ConstTransfer(1, PlasticW, priority=1), cp.StochasticTransfer(nr.triangular, [0.462, 0.556, 0.650], MMSW, priority=2)] Plas_Ex_y0 = cp.TransferDistribution(nr.triangular, [0.221,0.23,0.239]) Plas_Ex_y1 = cp.TransferDistribution(nr.triangular, [0.239,0.249,0.259]) Plas_Ex_y2 = cp.TransferDistribution(nr.triangular, [0.258,0.268,0.27899999999999997]) Plas_Ex_y3 = cp.TransferDistribution(nr.triangular, [0.276,0.28800000000000003,0.299]) Plas_Ex_y4 = cp.TransferDistribution(nr.triangular, [0.294,0.307,0.319]) Plas_Ex_y5 = cp.TransferDistribution(nr.triangular, [0.313,0.326,0.33899999999999997]) Plas_Ex_y6 = cp.TransferDistribution(nr.triangular, [0.331,0.345,0.359]) Plas_Ex_y7 = cp.TransferDistribution(nr.triangular, [0.35,0.364,0.379]) Plas_Ex_y8 = cp.TransferDistribution(nr.triangular, [0.368,0.38299999999999995,0.39899999999999997]) Plas_Ex_y9 = cp.TransferDistribution(nr.triangular, [0.386,0.40299999999999997,0.419]) Plas_Ex_y10 = cp.TransferDistribution(nr.triangular, [0.405,0.42200000000000004,0.439]) Plas_Ex_y11 = cp.TransferDistribution(nr.triangular, [0.423,0.441,0.45799999999999996]) Plas_Ex_y12 = cp.TransferDistribution(nr.triangular, [0.442,0.46,0.478]) Plas_Ex_y13 = cp.TransferDistribution(nr.triangular, [0.45399999999999996,0.473,0.49200000000000005]) Plas_Ex_y14 = cp.TransferDistribution(nr.triangular, [0.467,0.48700000000000004,0.506]) Plas_Ex_y15 = cp.TransferDistribution(nr.triangular, [0.48,0.5,0.52]) Plas_Ex_y16 = cp.TransferDistribution(nr.triangular, [0.48,0.5,0.52]) Plas_Ex_y17 = cp.TransferDistribution(nr.triangular, [0.48,0.5,0.52]) Plas_Ex_y18 = cp.TransferDistribution(nr.triangular, [0.48,0.5,0.52]) Plas_Ex_y19 = cp.TransferDistribution(nr.triangular, [0.48,0.5,0.52]) Plas_Ex_y20 = cp.TransferDistribution(nr.triangular, [0.48,0.5,0.52]) PlasticW.transfers = [cp.ConstTransfer(1, Sorting_PlasticW, priority=1)] Sorting_PlasticW.transfers = [cp.TimeDependendDistributionTransfer([Plas_Ex_y0, Plas_Ex_y1,Plas_Ex_y2,Plas_Ex_y3,Plas_Ex_y4,Plas_Ex_y5, Plas_Ex_y6,Plas_Ex_y7,Plas_Ex_y8,Plas_Ex_y9,Plas_Ex_y10, Plas_Ex_y11,Plas_Ex_y12,Plas_Ex_y13,Plas_Ex_y14,Plas_Ex_y15, Plas_Ex_y16,Plas_Ex_y17,Plas_Ex_y18,Plas_Ex_y19,Plas_Ex_y20], SortPlasW_Export_P, Sorting_PlasticW, priority=2), cp.ConstTransfer(1, Reprocess_Plas_P, priority=1)] # In[32]: # Solar # IRENA, 2016, End-of-life management - Solar Photovoltaic Panels Report: Average lifetime of a solar panel is 30 years. # Since our model covers only 20 years and we assume no TiO2 was included in solar panels before 2000, we don't need to # look at what happens to waste, since the nano will still be in stock in 2020/2030. # So I assume everything goes to waste at 31 years, and everything is reprocessed, only to give the model transfers. Solar_EoL.localRelease = cp.ListRelease([0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1]) Solar_EoL.transfers = [cp.ConstTransfer(1, SolarWaste, priority=1)] # To change if the model goes further than 2030. SolarWaste.transfers = [cp.ConstTransfer(1, Reprocess_Solar_P, priority=1)] # In[33]: # Electronics Electronics_Use.localRelease = cp.FixedRateRelease(0.125) Electronics_EoL.localRelease = cp.ListRelease([0.002457901,0.004936706,0.01218547,0.02617355,0.04892212, 0.07957497,0.112637,0.1387466,0.1487314,0.1387466, 0.112637,0.07957497,0.04892212,0.02617355,0.01218547, 0.004936706,0.002457901]) Electronics_Use.transfers = [cp.ConstTransfer(1, Electronics_WW, priority=1)] Electronics_WW.transfers = [cp.StochasticTransfer(tf.TriangTruncDet, [0.475, 0.95, 1.425, 1, 0, 1], Electronics_WW_N, priority=1), cp.StochasticTransfer(nr.triangular, [0.025, 0.05, 0.075], Electronics_WW_M, priority=1)] Electronics_EoL.transfers = [cp.ConstTransfer(1, WEEE, priority=1), cp.StochasticTransfer(nr.triangular, [0.085, 0.17, 0.255], MMSW, priority=2)] # In[34]: # Cleaning Agents CleanAgents_Use.localRelease = cp.ListRelease([1.0]) CleanAgents_EoL.localRelease = cp.ListRelease([1.0]) CleanAgents_Use.transfers = [cp.ConstTransfer(1, CleanAgents_WW, priority=1)] CleanAgents_WW.transfers = [cp.ConstTransfer(1, CleanAgents_WW_N, priority=1)] CleanAgents_EoL.transfers = [cp.ConstTransfer(1, PackW, priority=1)] # In[35]: # Food Food_Use.localRelease = cp.ListRelease([1.0]) Food_EoL.localRelease = cp.ListRelease([1.0]) Food_Use.transfers = [cp.ConstTransfer(1, Food_WW, priority=1)] Food_WW.transfers = [cp.ConstTransfer(1, Food_WW_N, priority=1)] Food_EoL.transfers = [cp.StochasticTransfer(tf.TriangTruncDet,[0.495, 0.99, 1.485, 1, 0, 1], MMSW, priority=1), cp.StochasticTransfer(nr.triangular, [0.005, 0.01, 0.015], PackW, priority=1) ] # In[36]: # Textiles Textiles_Use.localRelease = cp.ListRelease([0.5,0.3,0.2]) Textiles_EoL.localRelease = cp.ListRelease([0.01222447,0.2144029,0.5467453,0.2144029,0.01222447]) Textiles_Use.transfers = [cp.ConstTransfer(1, Textiles_WW, priority=1), cp.StochasticTransfer(nr.triangular, [0.10,0.20,0.30], Textiles_Air, priority=2)] Textiles_WW.transfers = [cp.StochasticTransfer(nr.uniform, [0.25, 0.75], Textiles_WW_N, priority=1), cp.StochasticTransfer(nr.uniform, [0.25, 0.75], Textiles_WW_M, priority=1)] Textiles_Air.transfers = [cp.StochasticTransfer(nr.uniform, [1e-6, 1], Textiles_Air_N, priority=1), cp.StochasticTransfer(nr.uniform, [1e-6, 1], Textiles_Air_M, priority=1)] Textiles_EoL.transfers = [cp.ConstTransfer(1, TextW, priority=1)] # In[37]: ############################################################################## ########################## WASTE MANAGEMENT ################################## ############################################################################## ################## sampling of the removal data from different sources ######## s = 10000 #sampling size used for sampling in uncertainty ranges of several sources/values # In[38]: ################################# WWTP #################################### # In[39]: # Gather all pristine flows to wastewater PersCare_WW_N.transfers = [cp.ConstTransfer(1, WW_N, priority=1)] Paints_WW_N.transfers = [cp.ConstTransfer(1, WW_N, priority=1)] In_WW_N.transfers = [cp.ConstTransfer(1, WW_N, priority=1)] Cement_WW_N.transfers = [cp.ConstTransfer(1, WW_N, priority=1)] Glass_WW_N.transfers = [cp.ConstTransfer(1, WW_N, priority=1)] Ceramics_WW_N.transfers = [cp.ConstTransfer(1, WW_N, priority=1)] RubPlas_WW_N.transfers = [cp.ConstTransfer(1, WW_N, priority=1)] Electronics_WW_N.transfers = [cp.ConstTransfer(1, WW_N, priority=1)] CleanAgents_WW_N.transfers = [cp.ConstTransfer(1, WW_N, priority=1)] Food_WW_N.transfers = [cp.ConstTransfer(1, WW_N, priority=1)] Textiles_WW_N.transfers = [cp.ConstTransfer(1, WW_N, priority=1)] # Gather all matrix-embedded flows to wastewater Paints_WW_M.transfers = [cp.ConstTransfer(1, WW_M, priority=1)] In_WW_M.transfers = [cp.ConstTransfer(1, WW_M, priority=1)] Cement_WW_M.transfers = [cp.ConstTransfer(1, WW_M, priority=1)] Glass_WW_M.transfers = [cp.ConstTransfer(1, WW_M, priority=1)] Ceramics_WW_M.transfers = [cp.ConstTransfer(1, WW_M, priority=1)] RubPlas_WW_M.transfers = [cp.ConstTransfer(1, WW_M, priority=1)] Electronics_WW_M.transfers = [cp.ConstTransfer(1, WW_M, priority=1)] Textiles_WW_M.transfers = [cp.ConstTransfer(1, WW_M, priority=1)] # In[40]: # From wastewater to non-sewer WW_NS_y0 = cp.TransferDistribution(nr.triangular, [0.0815, 0.163, 0.2445]) WW_NS_y1 = cp.TransferDistribution(nr.triangular, [0.075, 0.15, 0.225]) WW_NS_y2 = cp.TransferDistribution(nr.triangular, [0.0775, 0.155, 0.2325]) WW_NS_y3 = cp.TransferDistribution(nr.triangular, [0.0765, 0.153, 0.2295]) WW_NS_y4 = cp.TransferDistribution(nr.triangular, [0.076, 0.152, 0.228]) WW_NS_y5 = cp.TransferDistribution(nr.triangular, [0.0785, 0.157, 0.2355]) WW_NS_y6 = cp.TransferDistribution(nr.triangular, [0.067, 0.134, 0.201]) WW_NS_y7 = cp.TransferDistribution(nr.triangular, [0.0715, 0.143, 0.2145]) WW_NS_y8 = cp.TransferDistribution(nr.triangular, [0.071, 0.142, 0.213]) WW_NS_y9 = cp.TransferDistribution(nr.triangular, [0.072, 0.144, 0.216]) WW_NS_y10 = cp.TransferDistribution(nr.triangular, [0.0715, 0.143, 0.2145]) WW_NS_y11 = cp.TransferDistribution(nr.triangular, [0.074, 0.148, 0.222]) WW_NS_y12 = cp.TransferDistribution(nr.triangular, [0.0725, 0.145, 0.2175]) WW_NS_y13 = cp.TransferDistribution(nr.triangular, [0.072, 0.144, 0.216]) WW_NS_y14 = cp.TransferDistribution(nr.triangular, [0.071, 0.142, 0.213]) WW_NS_y15 = cp.TransferDistribution(nr.triangular, [0.0715, 0.143, 0.2145]) WW_NS_y16 = cp.TransferDistribution(nr.triangular, [0.0715, 0.143, 0.2145]) WW_NS_y17 = cp.TransferDistribution(nr.triangular, [0.0715, 0.143, 0.2145]) WW_NS_y18 = cp.TransferDistribution(nr.triangular, [0.0715, 0.143, 0.2145]) WW_NS_y19 = cp.TransferDistribution(nr.triangular, [0.0715, 0.143, 0.2145]) WW_NS_y20 = cp.TransferDistribution(nr.triangular, [0.072, 0.144, 0.216]) WW_N.transfers = [cp.TimeDependendDistributionTransfer([WW_NS_y0, WW_NS_y1,WW_NS_y2,WW_NS_y3,WW_NS_y4,WW_NS_y5, WW_NS_y6,WW_NS_y7,WW_NS_y8,WW_NS_y9,WW_NS_y10, WW_NS_y11,WW_NS_y12,WW_NS_y13,WW_NS_y14,WW_NS_y15, WW_NS_y16,WW_NS_y17,WW_NS_y18,WW_NS_y19,WW_NS_y20], NoSewageSystem_N, WW_N, priority=2), cp.ConstTransfer(1, SewageSystem_N, priority=1)] WW_M.transfers = [cp.TimeDependendDistributionTransfer([WW_NS_y0, WW_NS_y1,WW_NS_y2,WW_NS_y3,WW_NS_y4,WW_NS_y5, WW_NS_y6,WW_NS_y7,WW_NS_y8,WW_NS_y9,WW_NS_y10, WW_NS_y11,WW_NS_y12,WW_NS_y13,WW_NS_y14,WW_NS_y15, WW_NS_y16,WW_NS_y17,WW_NS_y18,WW_NS_y19,WW_NS_y20], NoSewageSystem_M, WW_M, priority=2), cp.ConstTransfer(1, SewageSystem_M, priority=1)] WW_T.transfers = [cp.TimeDependendDistributionTransfer([WW_NS_y0, WW_NS_y1,WW_NS_y2,WW_NS_y3,WW_NS_y4,WW_NS_y5, WW_NS_y6,WW_NS_y7,WW_NS_y8,WW_NS_y9,WW_NS_y10, WW_NS_y11,WW_NS_y12,WW_NS_y13,WW_NS_y14,WW_NS_y15, WW_NS_y16,WW_NS_y17,WW_NS_y18,WW_NS_y19,WW_NS_y20], NoSewageSystem_T, WW_T, priority=2), cp.ConstTransfer(1, SewageSystem_T, priority=1)] NoSewageSystem_N.transfers = [cp.ConstTransfer(0.03, NoSew_SW_N, priority=1), cp.ConstTransfer(0.11, OnSiteTreat_N, priority=1)] NoSewageSystem_M.transfers = [cp.ConstTransfer(0.03, NoSew_SW_M, priority=1), cp.ConstTransfer(0.11, OnSiteTreat_M, priority=1)] NoSewageSystem_T.transfers = [cp.ConstTransfer(0.03, NoSew_SW_T, priority=1), cp.ConstTransfer(0.11, OnSiteTreat_T, priority=1)] OnSiteTreat_N.transfers = [cp.ConstTransfer(1, TreatIOnSite_N, priority=1)] OnSiteTreat_M.transfers = [cp.ConstTransfer(1, TreatIOnSite_M, priority=1)] OnSiteTreat_T.transfers = [cp.ConstTransfer(1, TreatIOnSite_T, priority=1)] # In[41]: # From sewer system to surface water (non-treated, NT) SS_NT_y0 = cp.TransferDistribution(nr.triangular, [0.0325, 0.065, 0.0975]) SS_NT_y1 = cp.TransferDistribution(nr.triangular, [0.034, 0.068, 0.102]) SS_NT_y2 = cp.TransferDistribution(nr.triangular, [0.028, 0.056, 0.084]) SS_NT_y3 = cp.TransferDistribution(nr.triangular, [0.0245, 0.049, 0.0735]) SS_NT_y4 = cp.TransferDistribution(nr.triangular, [0.0225, 0.045, 0.0675]) SS_NT_y5 = cp.TransferDistribution(nr.triangular, [0.018, 0.036, 0.054]) SS_NT_y6 = cp.TransferDistribution(nr.triangular, [0.0285, 0.057, 0.0855]) SS_NT_y7 = cp.TransferDistribution(nr.triangular, [0.024, 0.048, 0.072]) SS_NT_y8 = cp.TransferDistribution(nr.triangular, [0.0235, 0.047, 0.0705]) SS_NT_y9 = cp.TransferDistribution(nr.triangular, [0.0205, 0.041, 0.0615]) SS_NT_y10 = cp.TransferDistribution(nr.triangular, [0.0185, 0.037, 0.0555]) SS_NT_y11 = cp.TransferDistribution(nr.triangular, [0.016, 0.032, 0.048]) SS_NT_y12 = cp.TransferDistribution(nr.triangular, [0.014, 0.028, 0.042]) SS_NT_y13 = cp.TransferDistribution(nr.triangular, [0.0125, 0.025, 0.0375]) SS_NT_y14 = cp.TransferDistribution(nr.triangular, [0.0105, 0.021, 0.0315]) SS_NT_y15 = cp.TransferDistribution(nr.triangular, [0.0085, 0.017, 0.0255]) SS_NT_y16 = cp.TransferDistribution(nr.triangular, [0.007, 0.014, 0.021]) SS_NT_y17 = cp.TransferDistribution(nr.triangular, [0.005, 0.01, 0.015]) SS_NT_y18 = cp.TransferDistribution(nr.triangular, [0.0035, 0.007, 0.0105]) SS_NT_y19 = cp.TransferDistribution(nr.triangular, [0.0015, 0.003, 0.0045]) SS_NT_y20 = cp.TransferDistribution(nr.triangular, [0, 0, 0.0001]) # changed by hand to be different from 0 SewageSystem_N.transfers= [cp.StochasticTransfer(nr.triangular,[0.016, 0.032, 0.048], STPoverflow_N, priority=2), # overflow rate cp.StochasticTransfer(nr.triangular, [0.0325, 0.065, 0.0975], Sew_Subsurface_N, priority=2), #Exfiltration cp.TimeDependendDistributionTransfer([SS_NT_y0, SS_NT_y1,SS_NT_y2,SS_NT_y3,SS_NT_y4,SS_NT_y5, SS_NT_y6,SS_NT_y7,SS_NT_y8,SS_NT_y9,SS_NT_y10, SS_NT_y11,SS_NT_y12,SS_NT_y13,SS_NT_y14,SS_NT_y15, SS_NT_y16,SS_NT_y17,SS_NT_y18,SS_NT_y19,SS_NT_y20], Sew_SW_N, SewageSystem_N, priority=2), #discharges of sewer to nature (when no wwtp is installed) cp.ConstTransfer(1, WWTP_N, priority=1)] SewageSystem_M.transfers= [cp.StochasticTransfer(nr.triangular,[0.016, 0.032, 0.048], STPoverflow_M, priority=2), # overflow rate cp.StochasticTransfer(nr.triangular, [0.0325, 0.065, 0.0975], Sew_Subsurface_M, priority=2), #Exfiltration cp.TimeDependendDistributionTransfer([SS_NT_y0, SS_NT_y1,SS_NT_y2,SS_NT_y3,SS_NT_y4,SS_NT_y5, SS_NT_y6,SS_NT_y7,SS_NT_y8,SS_NT_y9,SS_NT_y10, SS_NT_y11,SS_NT_y12,SS_NT_y13,SS_NT_y14,SS_NT_y15, SS_NT_y16,SS_NT_y17,SS_NT_y18,SS_NT_y19,SS_NT_y20], Sew_SW_M, SewageSystem_M, priority=2), #discharges of sewer to nature (when no wwtp is installed) cp.ConstTransfer(1, WWTP_M, priority=1)] SewageSystem_T.transfers= [cp.StochasticTransfer(nr.triangular,[0.016, 0.032, 0.048], STPoverflow_T, priority=2), # overflow rate cp.StochasticTransfer(nr.triangular, [0.0325, 0.065, 0.0975], Sew_Subsurface_T, priority=2), #Exfiltration cp.TimeDependendDistributionTransfer([SS_NT_y0, SS_NT_y1,SS_NT_y2,SS_NT_y3,SS_NT_y4,SS_NT_y5, SS_NT_y6,SS_NT_y7,SS_NT_y8,SS_NT_y9,SS_NT_y10, SS_NT_y11,SS_NT_y12,SS_NT_y13,SS_NT_y14,SS_NT_y15, SS_NT_y16,SS_NT_y17,SS_NT_y18,SS_NT_y19,SS_NT_y20], Sew_SW_T, SewageSystem_T, priority=2), #discharges of sewer to nature (when no wwtp is installed) cp.ConstTransfer(1, WWTP_T, priority=1)] # In[42]: # Levels of treatment in WWTP WWTP_TI_y0 = cp.TransferDistribution(nr.triangular, [0.0295, 0.059, 0.0885]) WWTP_TI_y1 = cp.TransferDistribution(nr.triangular, [0.0255, 0.051, 0.0765]) WWTP_TI_y2 = cp.TransferDistribution(nr.triangular, [0.021, 0.042, 0.063]) WWTP_TI_y3 = cp.TransferDistribution(nr.triangular, [0.0175, 0.035, 0.0525]) WWTP_TI_y4 = cp.TransferDistribution(nr.triangular, [0.014, 0.028, 0.042]) WWTP_TI_y5 = cp.TransferDistribution(nr.triangular, [0.011, 0.022, 0.033]) WWTP_TI_y6 = cp.TransferDistribution(nr.triangular, [0.012, 0.024, 0.036]) WWTP_TI_y7 = cp.TransferDistribution(nr.triangular, [0.0105, 0.021, 0.0315]) WWTP_TI_y8 = cp.TransferDistribution(nr.triangular, [0.011, 0.022, 0.033]) WWTP_TI_y9 = cp.TransferDistribution(nr.triangular, [0.008, 0.016, 0.024]) WWTP_TI_y10 = cp.TransferDistribution(nr.triangular, [0.0075, 0.015, 0.0225]) WWTP_TI_y11 = cp.TransferDistribution(nr.triangular, [0.0065, 0.013, 0.0195]) WWTP_TI_y12 = cp.TransferDistribution(nr.triangular, [0.0055, 0.011, 0.0165]) WWTP_TI_y13 = cp.TransferDistribution(nr.triangular, [0.005, 0.01, 0.015]) WWTP_TI_y14 = cp.TransferDistribution(nr.triangular, [0.004, 0.008, 0.012]) WWTP_TI_y15 = cp.TransferDistribution(nr.triangular, [0.0035, 0.007, 0.0105]) WWTP_TI_y16 = cp.TransferDistribution(nr.triangular, [0.0025, 0.005, 0.0075]) WWTP_TI_y17 = cp.TransferDistribution(nr.triangular, [0.002, 0.004, 0.006]) WWTP_TI_y18 = cp.TransferDistribution(nr.triangular, [0.0015, 0.003, 0.0045]) WWTP_TI_y19 = cp.TransferDistribution(nr.triangular, [0.0005, 0.001, 0.0015]) WWTP_TI_y20 = cp.TransferDistribution(nr.triangular, [0, 0, 0.00001]) # adapted by hand to be not equal 0 WWTP_TII_y0 = cp.TransferDistribution(nr.triangular, [0.1725, 0.345, 0.5175]) WWTP_TII_y1 = cp.TransferDistribution(nr.triangular, [0.1615, 0.323, 0.4845]) WWTP_TII_y2 = cp.TransferDistribution(nr.triangular, [0.161, 0.322, 0.483]) WWTP_TII_y3 = cp.TransferDistribution(nr.triangular, [0.1535, 0.307, 0.4605]) WWTP_TII_y4 = cp.TransferDistribution(nr.triangular, [0.145, 0.29, 0.435]) WWTP_TII_y5 = cp.TransferDistribution(nr.triangular, [0.123, 0.246, 0.369]) WWTP_TII_y6 = cp.TransferDistribution(nr.triangular, [0.1305, 0.261, 0.3915]) WWTP_TII_y7 = cp.TransferDistribution(nr.triangular, [0.124, 0.248, 0.372]) WWTP_TII_y8 = cp.TransferDistribution(nr.triangular, [0.118, 0.236, 0.354]) WWTP_TII_y9 = cp.TransferDistribution(nr.triangular, [0.1125, 0.225, 0.3375]) WWTP_TII_y10 = cp.TransferDistribution(nr.triangular, [0.108, 0.216, 0.324]) WWTP_TII_y11 = cp.TransferDistribution(nr.triangular, [0.105, 0.21, 0.315]) WWTP_TII_y12 = cp.TransferDistribution(nr.triangular, [0.1025, 0.205, 0.3075]) WWTP_TII_y13 = cp.TransferDistribution(nr.triangular, [0.0975, 0.195, 0.2925]) WWTP_TII_y14 = cp.TransferDistribution(nr.triangular, [0.0925, 0.185, 0.2775]) WWTP_TII_y15 = cp.TransferDistribution(nr.triangular, [0.0875, 0.175, 0.2625]) WWTP_TII_y16 = cp.TransferDistribution(nr.triangular, [0.083, 0.166, 0.249]) WWTP_TII_y17 = cp.TransferDistribution(nr.triangular, [0.0785, 0.157, 0.2355]) WWTP_TII_y18 = cp.TransferDistribution(nr.triangular, [0.074, 0.148, 0.222]) WWTP_TII_y19 = cp.TransferDistribution(nr.triangular, [0.0706, 0.139, 0.2074]) WWTP_TII_y20 = cp.TransferDistribution(nr.triangular, [0.0609, 0.109, 0.1571]) WWTP_N.transfers = [cp.TimeDependendDistributionTransfer([WWTP_TI_y0, WWTP_TI_y1,WWTP_TI_y2,WWTP_TI_y3,WWTP_TI_y4,WWTP_TI_y5, WWTP_TI_y6,WWTP_TI_y7,WWTP_TI_y8,WWTP_TI_y9,WWTP_TI_y10, WWTP_TI_y11,WWTP_TI_y12,WWTP_TI_y13,WWTP_TI_y14,WWTP_TI_y15, WWTP_TI_y16,WWTP_TI_y17,WWTP_TI_y18,WWTP_TI_y19,WWTP_TI_y20], TreatI_N, WWTP_N, priority=2), cp.TimeDependendDistributionTransfer([WWTP_TII_y0, WWTP_TII_y1,WWTP_TII_y2,WWTP_TII_y3,WWTP_TII_y4,WWTP_TII_y5, WWTP_TII_y6,WWTP_TII_y7,WWTP_TII_y8,WWTP_TII_y9,WWTP_TII_y10, WWTP_TII_y11,WWTP_TII_y12,WWTP_TII_y13,WWTP_TII_y14,WWTP_TII_y15, WWTP_TII_y16,WWTP_TII_y17,WWTP_TII_y18,WWTP_TII_y19,WWTP_TII_y20], TreatII_N, WWTP_N, priority=2), cp.ConstTransfer(1, TreatIII_N, priority=1)] WWTP_M.transfers = [cp.TimeDependendDistributionTransfer([WWTP_TI_y0, WWTP_TI_y1,WWTP_TI_y2,WWTP_TI_y3,WWTP_TI_y4,WWTP_TI_y5, WWTP_TI_y6,WWTP_TI_y7,WWTP_TI_y8,WWTP_TI_y9,WWTP_TI_y10, WWTP_TI_y11,WWTP_TI_y12,WWTP_TI_y13,WWTP_TI_y14,WWTP_TI_y15, WWTP_TI_y16,WWTP_TI_y17,WWTP_TI_y18,WWTP_TI_y19,WWTP_TI_y20], TreatI_M, WWTP_M, priority=2), cp.TimeDependendDistributionTransfer([WWTP_TII_y0, WWTP_TII_y1,WWTP_TII_y2,WWTP_TII_y3,WWTP_TII_y4,WWTP_TII_y5, WWTP_TII_y6,WWTP_TII_y7,WWTP_TII_y8,WWTP_TII_y9,WWTP_TII_y10, WWTP_TII_y11,WWTP_TII_y12,WWTP_TII_y13,WWTP_TII_y14,WWTP_TII_y15, WWTP_TII_y16,WWTP_TII_y17,WWTP_TII_y18,WWTP_TII_y19,WWTP_TII_y20], TreatII_M, WWTP_M, priority=2), cp.ConstTransfer(1, TreatIII_M, priority=1)] WWTP_T.transfers = [cp.TimeDependendDistributionTransfer([WWTP_TI_y0, WWTP_TI_y1,WWTP_TI_y2,WWTP_TI_y3,WWTP_TI_y4,WWTP_TI_y5, WWTP_TI_y6,WWTP_TI_y7,WWTP_TI_y8,WWTP_TI_y9,WWTP_TI_y10, WWTP_TI_y11,WWTP_TI_y12,WWTP_TI_y13,WWTP_TI_y14,WWTP_TI_y15, WWTP_TI_y16,WWTP_TI_y17,WWTP_TI_y18,WWTP_TI_y19,WWTP_TI_y20], TreatI_T, WWTP_T, priority=2), cp.TimeDependendDistributionTransfer([WWTP_TII_y0, WWTP_TII_y1,WWTP_TII_y2,WWTP_TII_y3,WWTP_TII_y4,WWTP_TII_y5, WWTP_TII_y6,WWTP_TII_y7,WWTP_TII_y8,WWTP_TII_y9,WWTP_TII_y10, WWTP_TII_y11,WWTP_TII_y12,WWTP_TII_y13,WWTP_TII_y14,WWTP_TII_y15, WWTP_TII_y16,WWTP_TII_y17,WWTP_TII_y18,WWTP_TII_y19,WWTP_TII_y20], TreatII_T, WWTP_T, priority=2), cp.ConstTransfer(1, TreatIII_T, priority=1)] # In[43]: NonRemovalData_TreatI = np.concatenate([ #inversed from removal data nr.uniform( 0, 0.1, int(s*0.87)), nr.triangular(0.1, 0.1, 0.12, int(s*0.13)), # Brant et al, (2010) nr.triangular(0.07, 0.1, 0.13, s), # Sousa et al. 2017 ]) TreatI_N.transfers = [cp.RandomChoiceTransfer(NonRemovalData_TreatI, STPeffluent_N, priority=2), cp.ConstTransfer(1,STPsludge_N , priority=1)] TreatI_T.transfers = [cp.RandomChoiceTransfer(NonRemovalData_TreatI, STPeffluent_T, priority=2), cp.ConstTransfer(1, STPsludge_N, priority=1)] TreatIOnSite_N.transfers = [cp.RandomChoiceTransfer(NonRemovalData_TreatI, OnSiteTreat_Subsurface_N, priority=2), cp.ConstTransfer(1, OnSiteSludge_N, priority=1)] TreatIOnSite_T.transfers = [cp.RandomChoiceTransfer(NonRemovalData_TreatI, OnSiteTreat_Subsurface_T, priority=2), cp.ConstTransfer(1, OnSiteSludge_T, priority=1)] # We make the assumption that matrix-embedded ENMs behave as microplastics in WWTP, since they have about the same size. # Three sources are available that studied the removal of MPs after primary treatment in Europe: # Murphy et al. 2016: 78.34% go to sludge, with a coefficient of variation of 8%. # Talvitie et al. 2015: from their data, we build a triangular distribution with min of 25.8%, mode of 50% and max of 92.5%. # Lares et al. 2018: from their data, we build a triangular distribution with min of 91.9%, mode of 99% and max of 99%. Removal_Murphy = nr.triangular(0.7834*(1-0.08), 0.7834, 0.7834*(1+0.08), 3333) Removal_Talvitie = nr.triangular(0.258, 0.5, 0.925, 3333) Removal_Lares = nr.triangular(0.919, 0.99, 0.99, 3333) RemovalData_TreatI_M = np.concatenate([Removal_Murphy, Removal_Talvitie, Removal_Lares]) TreatI_M.transfers = [cp.RandomChoiceTransfer(RemovalData_TreatI_M, STPsludge_M, priority=2), cp.ConstTransfer(1, STPeffluent_M, priority=1)] TreatIOnSite_M.transfers = [cp.RandomChoiceTransfer(RemovalData_TreatI_M, OnSiteSludge_M, priority=2), cp.ConstTransfer(1, OnSiteTreat_Subsurface_M, priority=1)] # In[44]: NonRemovalData_TreatII = np.concatenate([ #inversed from removal data nr.triangular(0.035, 0.05, 0.065, s), # Wang et al. 2012 ]) TreatII_N.transfers = [cp.RandomChoiceTransfer(NonRemovalData_TreatII, STPeffluent_N, priority=2), cp.ConstTransfer(1, STPsludge_N, priority=1)] TreatII_T.transfers = [cp.RandomChoiceTransfer(NonRemovalData_TreatII, STPeffluent_T, priority=2), cp.ConstTransfer(1, STPsludge_T, priority=1)] # We make the assumption that matrix-embedded ENMs behave as microplastics in WWTP, since they have about the same size. # Four sources are available that studied the removal of MPs after secondary treatment in Europe: # Murphy et al. 2016: 98.41% go to sludge, with a coefficient of variation of 16%. # Talvitie et al. 2015: from their data, we build a triangular distribution with min of 96.5%, mode of 97.8% and max of 98.1%. # Magnusson and Norén 2014: uniform distribution within 99.94% and 99.95%. # Lares et al. 2018: from their data, we build a triangular distribution with min of 81.1%, mode of 98.3% and max of 98.3%. Removal_Murphy = nr.triangular(0.9841*(1-0.16), 0.9841, 0.9841*(1+0.16), 2500) Removal_Talvitie = nr.triangular(0.965, 0.978, 0.981, 2500) Removal_Magnusson = nr.uniform(0.9994, 0.9995, 2500) Removal_Lares = nr.triangular(0.811, 0.983, 0.983, 2500) RemovalData_TreatII_M = np.concatenate([Removal_Murphy, Removal_Talvitie, Removal_Magnusson]) TreatII_M.transfers = [cp.RandomChoiceTransfer(RemovalData_TreatII_M, STPsludge_M, priority=2), cp.ConstTransfer(1, STPeffluent_M, priority=1)] # In[45]: NonRemovalData_TreatIII = np.concatenate([ # 1-removal data nr.triangular(0.014, 0.02, 0.026, s), # Sousa et al. 2020 nr.triangular(0.0056, 0.008, 0.0104, s), # Sousa et al. 2020 ]) TreatIII_N.transfers = [cp.RandomChoiceTransfer(NonRemovalData_TreatIII, STPeffluent_N, priority=2), cp.ConstTransfer(1, STPsludge_N, priority=1)] TreatIII_T.transfers = [cp.RandomChoiceTransfer(NonRemovalData_TreatIII, STPeffluent_T, priority=2), cp.ConstTransfer(1, STPsludge_T, priority=1)] # We make the assumption that matrix-embedded ENMs behave as microplastics in WWTP, since they have about the same size. # Two sources are available that studied the removal of MPs after secondary treatment in Europe: # Magni et al. 2019: from their data, we build a triangular distribution with min of 67.5%, mode of 84% and max of 91%. # Lares et al. 2018: from their data, we build a triangular distribution with min of 94.6%, mode of 99.3% and max of 99.8%. Removal_Magni = nr.triangular(0.675, 0.84, 0.91, 5000) Removal_Lares = nr.triangular(0.946, 0.993, 0.998, 5000) RemovalData_TreatIII_M = np.concatenate([Removal_Magni, Removal_Lares]) TreatIII_M.transfers = [cp.RandomChoiceTransfer(RemovalData_TreatIII_M, STPsludge_M, priority=2), cp.ConstTransfer(1, STPeffluent_M, priority=1)] # In[46]: # Transfers from sludge # To incineration # Eurostat, for year 2014: 73.45% incinerated in Austria. # Mean of this value for Europe for year 2014: 30.1% (this model). # It results a scaling factor of 2.416. # Before multiplication by 2.416, distributions must not exceed 0.414. SF_SG_WIP = 1 SG_WIIP_y0_data = (nr.triangular(0.1887, 0.209, 0.2293, s))*SF_SG_WIP SG_WIIP_y1_data = (nr.triangular(0.1822, 0.213, 0.2438, s))*SF_SG_WIP SG_WIIP_y2_data = (nr.triangular(0.1754, 0.217, 0.2586, s))*SF_SG_WIP SG_WIIP_y3_data = (nr.triangular(0.1792, 0.221, 0.2628, s))*SF_SG_WIP SG_WIIP_y4_data = (nr.triangular(0.1943, 0.225, 0.2557, s))*SF_SG_WIP SG_WIIP_y5_data = (nr.triangular(0.2098, 0.229, 0.2482, s))*SF_SG_WIP SG_WIIP_y6_data = (nr.triangular(0.2170, 0.256, 0.2950, s))*SF_SG_WIP SG_WIIP_y7_data = (nr.triangular(0.2228, 0.286, 0.3492, s))*SF_SG_WIP SG_WIIP_y8_data = (nr.triangular(0.2411, 0.317, 0.3929, s))*SF_SG_WIP SG_WIIP_y9_data = (nr.triangular(0.2526, 0.319, 0.3854, s))*SF_SG_WIP SG_WIIP_y10_data = (nr.triangular(0.2643, 0.321, 0.3777, s))*SF_SG_WIP SG_WIIP_y11_data = (nr.triangular(0.2231, 0.315, 0.4069, s))*SF_SG_WIP SG_WIIP_y12_data = (tf.TriangTruncDet(0.1868, 0.315, 0.4432, s, 0, 0.414))*SF_SG_WIP SG_WIIP_y13_data = (tf.TriangTruncDet(0.1780, 0.308, 0.4381, s, 0, 0.414))*SF_SG_WIP SG_WIIP_y14_data = (tf.TriangTruncDet(0.1693, 0.301, 0.4327, s, 0, 0.414))*SF_SG_WIP SG_WIIP_y15_data = (tf.TriangTruncDet(0.1614, 0.315, 0.4686, s, 0, 0.414))*SF_SG_WIP SG_WIIP_y16_data = (tf.TriangTruncDet(0.1645, 0.329, 0.4935, s, 0, 0.414))*SF_SG_WIP SG_WIIP_y17_data = (tf.TriangTruncDet(0.1715, 0.343, 0.5145, s, 0, 0.414))*SF_SG_WIP SG_WIIP_y18_data = (tf.TriangTruncDet(0.1785, 0.357, 0.5355, s, 0, 0.414))*SF_SG_WIP SG_WIIP_y19_data = (tf.TriangTruncDet(0.1901, 0.371, 0.5519, s, 0, 0.414))*SF_SG_WIP SG_WIIP_y20_data = (tf.TriangTruncDet(0.2171, 0.386, 0.5549, s, 0, 0.414))*SF_SG_WIP SG_WIIP_y0 = cp.TransferDistribution(nr.choice, [SG_WIIP_y0_data]) SG_WIIP_y1 = cp.TransferDistribution(nr.choice, [SG_WIIP_y1_data]) SG_WIIP_y2 = cp.TransferDistribution(nr.choice, [SG_WIIP_y2_data]) SG_WIIP_y3 = cp.TransferDistribution(nr.choice, [SG_WIIP_y3_data]) SG_WIIP_y4 = cp.TransferDistribution(nr.choice, [SG_WIIP_y4_data]) SG_WIIP_y5 = cp.TransferDistribution(nr.choice, [SG_WIIP_y5_data]) SG_WIIP_y6 = cp.TransferDistribution(nr.choice, [SG_WIIP_y6_data]) SG_WIIP_y7 = cp.TransferDistribution(nr.choice, [SG_WIIP_y7_data]) SG_WIIP_y8 = cp.TransferDistribution(nr.choice, [SG_WIIP_y8_data]) SG_WIIP_y9 = cp.TransferDistribution(nr.choice, [SG_WIIP_y9_data]) SG_WIIP_y10 = cp.TransferDistribution(nr.choice, [SG_WIIP_y10_data]) SG_WIIP_y11 = cp.TransferDistribution(nr.choice, [SG_WIIP_y11_data]) SG_WIIP_y12 = cp.TransferDistribution(nr.choice, [SG_WIIP_y12_data]) SG_WIIP_y13 = cp.TransferDistribution(nr.choice, [SG_WIIP_y13_data]) SG_WIIP_y14 = cp.TransferDistribution(nr.choice, [SG_WIIP_y14_data]) SG_WIIP_y15 = cp.TransferDistribution(nr.choice, [SG_WIIP_y15_data]) SG_WIIP_y16 = cp.TransferDistribution(nr.choice, [SG_WIIP_y16_data]) SG_WIIP_y17 = cp.TransferDistribution(nr.choice, [SG_WIIP_y17_data]) SG_WIIP_y18 = cp.TransferDistribution(nr.choice, [SG_WIIP_y18_data]) SG_WIIP_y19 = cp.TransferDistribution(nr.choice, [SG_WIIP_y19_data]) SG_WIIP_y20 = cp.TransferDistribution(nr.choice, [SG_WIIP_y20_data]) # Sludge to landfill # Eurostat, for year 2014: 1.99% landfilled in Austria. # Mean of this value for Europe for year 2014: 10.5% (this model). # It results a scaling factor of 0.189. SF_SG_LF = 1 SG_LF_y0_data = (nr.triangular(0.1887, 0.209, 0.2293, s))*SF_SG_LF SG_LF_y1_data = (nr.triangular(0.1780, 0.208, 0.2380, s))*SF_SG_LF SG_LF_y2_data = (nr.triangular(0.1681, 0.208, 0.2479, s))*SF_SG_LF SG_LF_y3_data = (nr.triangular(0.1678, 0.207, 0.2462, s))*SF_SG_LF SG_LF_y4_data = (nr.triangular(0.1779, 0.206, 0.2341, s))*SF_SG_LF SG_LF_y5_data = (nr.triangular(0.1878, 0.205, 0.2222, s))*SF_SG_LF SG_LF_y6_data = (nr.triangular(0.1729, 0.204, 0.2351, s))*SF_SG_LF SG_LF_y7_data = (nr.triangular(0.1581, 0.203, 0.2479, s))*SF_SG_LF SG_LF_y8_data = (nr.triangular(0.1544, 0.203, 0.2516, s))*SF_SG_LF SG_LF_y9_data = (nr.triangular(0.1457, 0.184, 0.2223, s))*SF_SG_LF SG_LF_y10_data = (nr.triangular(0.1375, 0.167, 0.1965, s))*SF_SG_LF SG_LF_y11_data = (nr.triangular(0.1069, 0.151, 0.1951, s))*SF_SG_LF SG_LF_y12_data = (nr.triangular(0.0818, 0.138, 0.1942, s))*SF_SG_LF SG_LF_y13_data = (nr.triangular(0.0699, 0.121, 0.1721, s))*SF_SG_LF SG_LF_y14_data = (nr.triangular(0.0591, 0.105, 0.1509, s))*SF_SG_LF SG_LF_y15_data = (nr.triangular(0.0567, 0.102, 0.1473, s))*SF_SG_LF SG_LF_y16_data = (nr.triangular(0.0539, 0.098, 0.1421, s))*SF_SG_LF SG_LF_y17_data = (nr.triangular(0.0516, 0.095, 0.1384, s))*SF_SG_LF SG_LF_y18_data = (nr.triangular(0.0579, 0.091, 0.1241, s))*SF_SG_LF SG_LF_y19_data = (nr.triangular(0.0642, 0.088, 0.1118, s))*SF_SG_LF SG_LF_y20_data = (nr.triangular(0.0692, 0.084, 0.0988, s))*SF_SG_LF SG_LF_y0 = cp.TransferDistribution(nr.choice, [SG_LF_y0_data]) SG_LF_y1 = cp.TransferDistribution(nr.choice, [SG_LF_y1_data]) SG_LF_y2 = cp.TransferDistribution(nr.choice, [SG_LF_y2_data]) SG_LF_y3 = cp.TransferDistribution(nr.choice, [SG_LF_y3_data]) SG_LF_y4 = cp.TransferDistribution(nr.choice, [SG_LF_y4_data]) SG_LF_y5 = cp.TransferDistribution(nr.choice, [SG_LF_y5_data]) SG_LF_y6 = cp.TransferDistribution(nr.choice, [SG_LF_y6_data]) SG_LF_y7 = cp.TransferDistribution(nr.choice, [SG_LF_y7_data]) SG_LF_y8 = cp.TransferDistribution(nr.choice, [SG_LF_y8_data]) SG_LF_y9 = cp.TransferDistribution(nr.choice, [SG_LF_y9_data]) SG_LF_y10 = cp.TransferDistribution(nr.choice, [SG_LF_y10_data]) SG_LF_y11 = cp.TransferDistribution(nr.choice, [SG_LF_y11_data]) SG_LF_y12 = cp.TransferDistribution(nr.choice, [SG_LF_y12_data]) SG_LF_y13 = cp.TransferDistribution(nr.choice, [SG_LF_y13_data]) SG_LF_y14 = cp.TransferDistribution(nr.choice, [SG_LF_y14_data]) SG_LF_y15 = cp.TransferDistribution(nr.choice, [SG_LF_y15_data]) SG_LF_y16 = cp.TransferDistribution(nr.choice, [SG_LF_y16_data]) SG_LF_y17 = cp.TransferDistribution(nr.choice, [SG_LF_y17_data]) SG_LF_y18 = cp.TransferDistribution(nr.choice, [SG_LF_y18_data]) SG_LF_y19 = cp.TransferDistribution(nr.choice, [SG_LF_y19_data]) SG_LF_y20 = cp.TransferDistribution(nr.choice, [SG_LF_y20_data]) STPsludge_N.transfers =[cp.TimeDependendDistributionTransfer([SG_WIIP_y0, SG_WIIP_y1,SG_WIIP_y2,SG_WIIP_y3,SG_WIIP_y4,SG_WIIP_y5, SG_WIIP_y6,SG_WIIP_y7,SG_WIIP_y8,SG_WIIP_y9,SG_WIIP_y10, SG_WIIP_y11,SG_WIIP_y12,SG_WIIP_y13,SG_WIIP_y14,SG_WIIP_y15, SG_WIIP_y16,SG_WIIP_y17,SG_WIIP_y18,SG_WIIP_y19,SG_WIIP_y20], WIP_N, STPsludge_N, priority=2), cp.TimeDependendDistributionTransfer([SG_LF_y0, SG_LF_y1,SG_LF_y2,SG_LF_y3,SG_LF_y4,SG_LF_y5, SG_LF_y6,SG_LF_y7,SG_LF_y8,SG_LF_y9,SG_LF_y10, SG_LF_y11,SG_LF_y12,SG_LF_y13,SG_LF_y14,SG_LF_y15, SG_LF_y16,SG_LF_y17,SG_LF_y18,SG_LF_y19,SG_LF_y20], STPSlud_Landfill_N, STPsludge_N, priority=2), cp.ConstTransfer(1, STsoil_N, priority=1)] STPsludge_M.transfers =[cp.TimeDependendDistributionTransfer([SG_WIIP_y0, SG_WIIP_y1,SG_WIIP_y2,SG_WIIP_y3,SG_WIIP_y4,SG_WIIP_y5, SG_WIIP_y6,SG_WIIP_y7,SG_WIIP_y8,SG_WIIP_y9,SG_WIIP_y10, SG_WIIP_y11,SG_WIIP_y12,SG_WIIP_y13,SG_WIIP_y14,SG_WIIP_y15, SG_WIIP_y16,SG_WIIP_y17,SG_WIIP_y18,SG_WIIP_y19,SG_WIIP_y20], WIP_M, STPsludge_M, priority=2), cp.TimeDependendDistributionTransfer([SG_LF_y0, SG_LF_y1,SG_LF_y2,SG_LF_y3,SG_LF_y4,SG_LF_y5, SG_LF_y6,SG_LF_y7,SG_LF_y8,SG_LF_y9,SG_LF_y10, SG_LF_y11,SG_LF_y12,SG_LF_y13,SG_LF_y14,SG_LF_y15, SG_LF_y16,SG_LF_y17,SG_LF_y18,SG_LF_y19,SG_LF_y20], STPSlud_Landfill_M, STPsludge_M, priority=2), cp.ConstTransfer(1, STsoil_M, priority=1)] STPsludge_T.transfers =[cp.TimeDependendDistributionTransfer([SG_WIIP_y0, SG_WIIP_y1,SG_WIIP_y2,SG_WIIP_y3,SG_WIIP_y4,SG_WIIP_y5, SG_WIIP_y6,SG_WIIP_y7,SG_WIIP_y8,SG_WIIP_y9,SG_WIIP_y10, SG_WIIP_y11,SG_WIIP_y12,SG_WIIP_y13,SG_WIIP_y14,SG_WIIP_y15, SG_WIIP_y16,SG_WIIP_y17,SG_WIIP_y18,SG_WIIP_y19,SG_WIIP_y20], WIP_T, STPsludge_T, priority=2), cp.TimeDependendDistributionTransfer([SG_LF_y0, SG_LF_y1,SG_LF_y2,SG_LF_y3,SG_LF_y4,SG_LF_y5, SG_LF_y6,SG_LF_y7,SG_LF_y8,SG_LF_y9,SG_LF_y10, SG_LF_y11,SG_LF_y12,SG_LF_y13,SG_LF_y14,SG_LF_y15, SG_LF_y16,SG_LF_y17,SG_LF_y18,SG_LF_y19,SG_LF_y20], STPSlud_Landfill_T, STPsludge_T, priority=2), cp.ConstTransfer(1, STsoil_T, priority=1)] # In[47]: ########################### COLLECTION RATE ################################# # In[48]: # Mixed Municipal Solid Waste # for MSW there is first an array built with 's' sampled values from each distribution of the different sources # In a second step from this sampled data, one value is selected and multiplied by the scaling factor for the UK. # Eurostat, for year 2014: 11.25% landfilled in Austria. # Mean of this value for Europe for year 2014: 47.87% (this model). # It results a scaling factor of 0.235. SF_MW_LF = 1 MSW_LF_Data_y0 = SF_MW_LF*(np.concatenate([nr.triangular(0.461, 0.922, 1.383, s), nr.triangular(0.411, 0.822, 1.233, s), nr.triangular(0.3978, 0.765, 1.1628, s)] )) MSW_LF_Data_y1 = SF_MW_LF*(np.concatenate([nr.triangular(0.4495, 0.899, 1.3485, s), nr.triangular(0.4035, 0.807, 1.2105, s), nr.triangular(0.38948, 0.749, 1.1385, s)] )) MSW_LF_Data_y2 = SF_MW_LF*(np.concatenate([nr.triangular(0.4708, 0.877, 1.3478, s), nr.triangular(0.3955, 0.791, 1.1865, s), nr.triangular(0.3817, 0.734, 1.11568, s)] )) MSW_LF_Data_y3 = SF_MW_LF*(np.concatenate([nr.triangular(0.5231, 0.855, 1.2056, s), nr.triangular(0.4035, 0.773, 1.1765, s), nr.triangular(0.3760, 0.723, 1.099, s)] )) MSW_LF_Data_y4 = SF_MW_LF*(np.concatenate([nr.triangular(0.5713, 0.832, 1.1149, s), nr.triangular(0.4347, 0.76, 1.1947, s), nr.triangular(0.3645, 0.701, 1.0655, s)] )) MSW_LF_Data_y5 = SF_MW_LF*(np.concatenate([nr.triangular(0.6169, 0.81, 1.0287, s), nr.triangular(0.3826, 0.733, 1.0629, s), nr.triangular(0.3484, 0.67, 0.9916, s)] )) MSW_LF_Data_y6 = SF_MW_LF*(np.concatenate([nr.triangular(0.7637, 0.913, 1.0956, s), nr.triangular(0.4021, 0.703, 1.1051, s), nr.triangular(0.3406, 0.655, 0.9694, s)] )) MSW_LF_Data_y7 = SF_MW_LF*(np.concatenate([nr.triangular(0.6412, 0.822, 1.0357, s), nr.triangular(0.3581, 0.686, 0.9947, s), nr.triangular(0.33488, 0.644, 0.95312, s)] )) MSW_LF_Data_y8 = SF_MW_LF*(np.concatenate([nr.triangular(0.5995, 0.728, 0.8565, s), nr.triangular(0.3815, 0.667, 0.9525, s), nr.triangular(0.31772, 0.611, 0.90428, s)] )) MSW_LF_Data_y9 = SF_MW_LF*(np.concatenate([nr.triangular(0.5175, 0.658, 0.7985, s), nr.triangular(0.3377, 0.647, 0.9563, s), nr.triangular(0.312, 0.6, 0.888, s)] )) MSW_LF_Data_y10 = SF_MW_LF*(np.concatenate([nr.triangular(0.5072, 0.597, 0.6868, s), nr.triangular(0.3586, 0.627, 0.8954, s), nr.triangular(0.3042, 0.585, 0.8658, s)] )) MSW_LF_Data_y11 = SF_MW_LF*(np.concatenate([nr.triangular(0.4582, 0.57, 0.6818, s), nr.triangular(0.3054, 0.585, 0.8646, s), nr.triangular(0.28652, 0.551, 0.81548, s)] )) MSW_LF_Data_y12 = SF_MW_LF*(np.concatenate([nr.triangular(0.4193, 0.553, 0.6867, s), nr.triangular(0.3077, 0.538, 0.7683, s), nr.triangular(0.27456, 0.528, 0.78144, s)] )) MSW_LF_Data_y13 = SF_MW_LF*(np.concatenate([nr.triangular(0.4355, 0.536, 0.6365, s), nr.triangular(0.2662, 0.51, 0.7538, s), nr.triangular(0.25792, 0.496, 0.73408, s)] )) MSW_LF_Data_y14 = SF_MW_LF*(np.concatenate([nr.triangular(0.4200, 0.487, 0.5540, s), nr.triangular(0.2740, 0.479, 0.6840, s), nr.triangular(0.24596, 0.473, 0.70004, s)] )) MSW_LF_Data_y15 = SF_MW_LF*(np.concatenate([nr.triangular(0.3831, 0.474, 0.5649, s), nr.triangular(0.2219, 0.425, 0.6282, s), nr.triangular(0.23296, 0.448, 0.66304, s)] )) MSW_LF_Data_y16 = SF_MW_LF*(np.concatenate([nr.triangular(0.3468, 0.46, 0.5732, s), nr.triangular(0.1875, 0.375, 0.5625, s), nr.triangular(0.206, 0.412, 0.618, s)] )) MSW_LF_Data_y17 = SF_MW_LF*(np.concatenate([nr.triangular(0.3113, 0.445, 0.5787, s), nr.triangular(0.1595, 0.319, 0.4785, s), nr.triangular(0.1895, 0.379, 0.5685, s)] )) MSW_LF_Data_y18 = SF_MW_LF*(np.concatenate([nr.triangular(0.3197, 0.429, 0.5383, s), nr.triangular(0.128, 0.256, 0.384, s), nr.triangular(0.1725, 0.345, 0.5175, s)] )) MSW_LF_Data_y19 = SF_MW_LF*(np.concatenate([nr.triangular(0.3250, 0.411, 0.4970, s), nr.triangular(0.0925, 0.185, 0.2775, s), nr.triangular(0.154, 0.308, 0.462, s)] )) MSW_LF_Data_y20 = SF_MW_LF*(np.concatenate([nr.triangular(0.3271, 0.391, 0.4549, s), nr.triangular(0.0515, 0.103, 0.1545, s), nr.triangular(0.1345, 0.269, 0.4035, s)] )) MSW_LF_y0 = cp.TransferDistribution(nr.choice, [MSW_LF_Data_y0]) MSW_LF_y1 = cp.TransferDistribution(nr.choice, [MSW_LF_Data_y1]) MSW_LF_y2 = cp.TransferDistribution(nr.choice, [MSW_LF_Data_y2]) MSW_LF_y3 = cp.TransferDistribution(nr.choice, [MSW_LF_Data_y3]) MSW_LF_y4 = cp.TransferDistribution(nr.choice, [MSW_LF_Data_y4]) MSW_LF_y5 = cp.TransferDistribution(nr.choice, [MSW_LF_Data_y5]) MSW_LF_y6 = cp.TransferDistribution(nr.choice, [MSW_LF_Data_y6]) MSW_LF_y7 = cp.TransferDistribution(nr.choice, [MSW_LF_Data_y7]) MSW_LF_y8 = cp.TransferDistribution(nr.choice, [MSW_LF_Data_y8]) MSW_LF_y9 = cp.TransferDistribution(nr.choice, [MSW_LF_Data_y9]) MSW_LF_y10 = cp.TransferDistribution(nr.choice, [MSW_LF_Data_y10]) MSW_LF_y11 = cp.TransferDistribution(nr.choice, [MSW_LF_Data_y11]) MSW_LF_y12 = cp.TransferDistribution(nr.choice, [MSW_LF_Data_y12]) MSW_LF_y13 = cp.TransferDistribution(nr.choice, [MSW_LF_Data_y13]) MSW_LF_y14 = cp.TransferDistribution(nr.choice, [MSW_LF_Data_y14]) MSW_LF_y15 = cp.TransferDistribution(nr.choice, [MSW_LF_Data_y15]) MSW_LF_y16 = cp.TransferDistribution(nr.choice, [MSW_LF_Data_y16]) MSW_LF_y17 = cp.TransferDistribution(nr.choice, [MSW_LF_Data_y17]) MSW_LF_y18 = cp.TransferDistribution(nr.choice, [MSW_LF_Data_y18]) MSW_LF_y19 = cp.TransferDistribution(nr.choice, [MSW_LF_Data_y19]) MSW_LF_y20 = cp.TransferDistribution(nr.choice, [MSW_LF_Data_y20]) # if changes are made here, don't forget to update the transfers for sorting (same same TCs) MMSW.transfers = [cp.TimeDependendDistributionTransfer([MSW_LF_y0, MSW_LF_y1,MSW_LF_y2,MSW_LF_y3,MSW_LF_y4,MSW_LF_y5, MSW_LF_y6,MSW_LF_y7,MSW_LF_y8,MSW_LF_y9,MSW_LF_y10, MSW_LF_y11,MSW_LF_y12,MSW_LF_y13,MSW_LF_y14,MSW_LF_y15, MSW_LF_y16,MSW_LF_y17,MSW_LF_y18,MSW_LF_y19,MSW_LF_y20], MMSW_Landfill_P, MMSW, priority=2), # From MMSW to Landfill cp.ConstTransfer(1, WIP_P, priority=1)] # In[49]: # Packaging Waste # Statusbericht 20115, for year 2014: 44.32% goes to MMSW in Austria. # Mean of this value for Europe in year 2014: 21.3% (this model). # It results a scaling factor of 2.081. # The overall distributions must not exceed 1, so 0.481 (=1/2.081) before multiplication by 2.081. SF_PW_MW = 1 PW_MW_y0_data = SF_PW_MW*(tf.TriangTruncDet(0.2925, 0.52, 0.7475, s, 0, 0.481)) PW_MW_y1_data = SF_PW_MW*(tf.TriangTruncDet(0.2998, 0.5, 0.7003, s, 0, 0.481)) PW_MW_y2_data = SF_PW_MW*(tf.TriangTruncDet(0.3030, 0.476, 0.6490, s, 0, 0.481)) PW_MW_y3_data = SF_PW_MW*(tf.TriangTruncDet(0.3488, 0.451, 0.5532, s, 0, 0.481)) PW_MW_y4_data = SF_PW_MW*(tf.TriangTruncDet(0.3516, 0.427, 0.5024, s, 0, 0.481)) PW_MW_y5_data = SF_PW_MW*(tf.TriangTruncDet(0.3513, 0.42, 0.4887, s, 0, 0.481)) PW_MW_y6_data = SF_PW_MW*(nr.triangular(0.3457, 0.407, 0.4683, s)) PW_MW_y7_data = SF_PW_MW*(nr.triangular(0.2887, 0.365, 0.4413, s)) PW_MW_y8_data = SF_PW_MW*(nr.triangular(0.2358, 0.322, 0.4082, s)) PW_MW_y9_data = SF_PW_MW*(nr.triangular(0.2166, 0.28, 0.3434, s)) PW_MW_y10_data = SF_PW_MW*(nr.triangular(0.1952, 0.237, 0.2788, s)) PW_MW_y11_data = SF_PW_MW*(nr.triangular(0.1869, 0.227, 0.2671, s)) PW_MW_y12_data = SF_PW_MW*(nr.triangular(0.1771, 0.215, 0.2529, s)) PW_MW_y13_data = SF_PW_MW*(nr.triangular(0.1163, 0.208, 0.2997, s)) PW_MW_y14_data = SF_PW_MW*(nr.triangular(0.1191, 0.213, 0.3069, s)) PW_MW_y15_data = SF_PW_MW*(nr.triangular(0.1191, 0.213, 0.3069, s)) PW_MW_y16_data = SF_PW_MW*(nr.triangular(0.0718, 0.141, 0.2102, s)) PW_MW_y17_data = SF_PW_MW*(nr.triangular(0.0585, 0.117, 0.1755, s)) PW_MW_y18_data = SF_PW_MW*(nr.triangular(0.0465, 0.093, 0.1395, s)) PW_MW_y19_data = SF_PW_MW*(nr.triangular(0.0345, 0.069, 0.1035, s)) PW_MW_y20_data = SF_PW_MW*(nr.triangular(0.023, 0.046, 0.069, s)) PW_MW_y0 = cp.TransferDistribution(nr.choice, [PW_MW_y0_data]) PW_MW_y1 = cp.TransferDistribution(nr.choice, [PW_MW_y1_data]) PW_MW_y2 = cp.TransferDistribution(nr.choice, [PW_MW_y2_data]) PW_MW_y3 = cp.TransferDistribution(nr.choice, [PW_MW_y3_data]) PW_MW_y4 = cp.TransferDistribution(nr.choice, [PW_MW_y4_data]) PW_MW_y5 = cp.TransferDistribution(nr.choice, [PW_MW_y5_data]) PW_MW_y6 = cp.TransferDistribution(nr.choice, [PW_MW_y6_data]) PW_MW_y7 = cp.TransferDistribution(nr.choice, [PW_MW_y7_data]) PW_MW_y8 = cp.TransferDistribution(nr.choice, [PW_MW_y8_data]) PW_MW_y9 = cp.TransferDistribution(nr.choice, [PW_MW_y9_data]) PW_MW_y10 = cp.TransferDistribution(nr.choice, [PW_MW_y10_data]) PW_MW_y11 = cp.TransferDistribution(nr.choice, [PW_MW_y11_data]) PW_MW_y12 = cp.TransferDistribution(nr.choice, [PW_MW_y12_data]) PW_MW_y13 = cp.TransferDistribution(nr.choice, [PW_MW_y13_data]) PW_MW_y14 = cp.TransferDistribution(nr.choice, [PW_MW_y14_data]) PW_MW_y15 = cp.TransferDistribution(nr.choice, [PW_MW_y15_data]) PW_MW_y16 = cp.TransferDistribution(nr.choice, [PW_MW_y16_data]) PW_MW_y17 = cp.TransferDistribution(nr.choice, [PW_MW_y17_data]) PW_MW_y18 = cp.TransferDistribution(nr.choice, [PW_MW_y18_data]) PW_MW_y19 = cp.TransferDistribution(nr.choice, [PW_MW_y19_data]) PW_MW_y20 = cp.TransferDistribution(nr.choice, [PW_MW_y20_data]) PackW.transfers = [cp.TimeDependendDistributionTransfer([PW_MW_y0, PW_MW_y1,PW_MW_y2,PW_MW_y3,PW_MW_y4,PW_MW_y5, PW_MW_y6,PW_MW_y7,PW_MW_y8,PW_MW_y9,PW_MW_y10, PW_MW_y11,PW_MW_y12,PW_MW_y13,PW_MW_y14,PW_MW_y15, PW_MW_y16,PW_MW_y17,PW_MW_y18,PW_MW_y19,PW_MW_y20], MMSW, PackW, priority=2), cp.ConstTransfer(1, Sorting_PackW, priority=1)] # In[50]: # WEEE # Eurostat, for year 2014: 93.16% goes to sorting in Austria. # Mean of this value for Europe in year 2014: 40% (this model). # It results a scaling factor of 2.329. # The overall distributions must not exceed 1, so 0.429 (=1/2.329) before multiplication by 2.329 SF_WE_SO = 1 WE_S_y0_data = SF_WE_SO*(nr.triangular(0.015, 0.03, 0.045, s)) WE_S_y1_data = SF_WE_SO*(nr.triangular(0.0335, 0.067, 0.1005, s)) WE_S_y2_data = SF_WE_SO*(nr.triangular(0.0515, 0.103, 0.1545, s)) WE_S_y3_data = SF_WE_SO*(nr.triangular(0.07, 0.14, 0.21, s)) WE_S_y4_data = SF_WE_SO*(nr.triangular(0.0907, 0.177, 0.2633, s)) WE_S_y5_data = SF_WE_SO*(tf.TriangTruncDet(0.1766, 0.314, 0.4514, s, 0, 0.429)) WE_S_y6_data = SF_WE_SO*(tf.TriangTruncDet(0.1912, 0.319, 0.4468, s, 0, 0.429)) WE_S_y7_data = SF_WE_SO*(nr.triangular(0.2231, 0.325, 0.4269, s)) WE_S_y8_data = SF_WE_SO*(nr.triangular(0.2718, 0.33, 0.3882, s)) WE_S_y9_data = SF_WE_SO*(nr.triangular(0.273, 0.35, 0.427, s)) WE_S_y10_data = SF_WE_SO*(tf.TriangTruncDet(0.3095, 0.37, 0.4305, s, 0, 0.429)) WE_S_y11_data = SF_WE_SO*(tf.TriangTruncDet(0.2566, 0.395, 0.5334, s, 0, 0.429)) WE_S_y12_data = SF_WE_SO*(tf.TriangTruncDet(0.2363, 0.42, 0.6038, s, 0, 0.429)) WE_S_y13_data = SF_WE_SO*(tf.TriangTruncDet(0.2663, 0.41, 0.5537, s, 0, 0.429)) WE_S_y14_data = SF_WE_SO*(tf.TriangTruncDet(0.3346, 0.4, 0.4654, s, 0, 0.429)) WE_S_y15_data = SF_WE_SO*(tf.TriangTruncDet(0.3582, 0.49, 0.6218, s, 0, 0.429)) WE_S_y16_data = SF_WE_SO*(tf.TriangTruncDet(0.3628, 0.58, 0.7972, s, 0, 0.429)) WE_S_y17_data = SF_WE_SO*(tf.TriangTruncDet(0.3819, 0.67, 0.9581, s, 0, 0.429)) WE_S_y18_data = SF_WE_SO*(tf.TriangTruncDet(0.3819, 0.67, 0.9581, s, 0, 0.429)) WE_S_y19_data = SF_WE_SO*(tf.TriangTruncDet(0.3819, 0.67, 0.9581, s, 0, 0.429)) WE_S_y20_data = SF_WE_SO*(tf.TriangTruncDet(0.3819, 0.67, 0.9581, s, 0, 0.429)) WE_S_y0 = cp.TransferDistribution(nr.choice, [WE_S_y0_data]) WE_S_y1 = cp.TransferDistribution(nr.choice, [WE_S_y1_data]) WE_S_y2 = cp.TransferDistribution(nr.choice, [WE_S_y2_data]) WE_S_y3 = cp.TransferDistribution(nr.choice, [WE_S_y3_data]) WE_S_y4 = cp.TransferDistribution(nr.choice, [WE_S_y4_data]) WE_S_y5 = cp.TransferDistribution(nr.choice, [WE_S_y5_data]) WE_S_y6 = cp.TransferDistribution(nr.choice, [WE_S_y6_data]) WE_S_y7 = cp.TransferDistribution(nr.choice, [WE_S_y7_data]) WE_S_y8 = cp.TransferDistribution(nr.choice, [WE_S_y8_data]) WE_S_y9 = cp.TransferDistribution(nr.choice, [WE_S_y9_data]) WE_S_y10 = cp.TransferDistribution(nr.choice, [WE_S_y10_data]) WE_S_y11 = cp.TransferDistribution(nr.choice, [WE_S_y11_data]) WE_S_y12 = cp.TransferDistribution(nr.choice, [WE_S_y12_data]) WE_S_y13 = cp.TransferDistribution(nr.choice, [WE_S_y13_data]) WE_S_y14 = cp.TransferDistribution(nr.choice, [WE_S_y14_data]) WE_S_y15 = cp.TransferDistribution(nr.choice, [WE_S_y15_data]) WE_S_y16 = cp.TransferDistribution(nr.choice, [WE_S_y16_data]) WE_S_y17 = cp.TransferDistribution(nr.choice, [WE_S_y17_data]) WE_S_y18 = cp.TransferDistribution(nr.choice, [WE_S_y18_data]) WE_S_y19 = cp.TransferDistribution(nr.choice, [WE_S_y19_data]) WE_S_y20 = cp.TransferDistribution(nr.choice, [WE_S_y20_data]) WEEE.transfers = [cp.TimeDependendDistributionTransfer([WE_S_y0, WE_S_y1,WE_S_y2,WE_S_y3,WE_S_y4,WE_S_y5, WE_S_y6,WE_S_y7,WE_S_y8,WE_S_y9,WE_S_y10, WE_S_y11,WE_S_y12,WE_S_y13,WE_S_y14,WE_S_y15, WE_S_y16,WE_S_y17,WE_S_y18,WE_S_y19,WE_S_y20], Sorting_WEEE, WEEE, priority=2), cp.ConstTransfer(1, MMSW, priority=1)] # In[51]: # Textile Waste # Eurostat, for year 2014: 28.82% goes to sorting in the UK. # Mean of this value for Europe in year 2014: 30% (this model). # It results a scaling factor of 0.961. SF_TW_SO = 1 TW_S_y0_data = (nr.triangular(0.024, 0.048, 0.072, s))*SF_TW_SO TW_S_y1_data = (nr.triangular(0.0345, 0.069, 0.1035, s))*SF_TW_SO TW_S_y2_data = (nr.triangular(0.0455, 0.091, 0.1365, s))*SF_TW_SO TW_S_y3_data = (nr.triangular(0.0574, 0.112, 0.1666, s))*SF_TW_SO TW_S_y4_data = (nr.triangular(0.1515, 0.19, 0.2285, s))*SF_TW_SO TW_S_y5_data = (nr.triangular(0.1486, 0.198, 0.2474, s))*SF_TW_SO TW_S_y6_data = (nr.triangular(0.1443, 0.205, 0.2657, s))*SF_TW_SO TW_S_y7_data = (nr.triangular(0.1400, 0.213, 0.2860, s))*SF_TW_SO TW_S_y8_data = (nr.triangular(0.1343, 0.22, 0.3057, s))*SF_TW_SO TW_S_y9_data = (nr.triangular(0.1513, 0.228, 0.3047, s))*SF_TW_SO TW_S_y10_data = (nr.triangular(0.1685, 0.235, 0.3015, s))*SF_TW_SO TW_S_y11_data = (nr.triangular(0.1872, 0.243, 0.2988, s))*SF_TW_SO TW_S_y12_data = (nr.triangular(0.2059, 0.25, 0.2941, s))*SF_TW_SO TW_S_y13_data = (nr.triangular(0.2139, 0.275, 0.3361, s))*SF_TW_SO TW_S_y14_data = (nr.triangular(0.2197, 0.3, 0.3804, s))*SF_TW_SO TW_S_y15_data = (nr.triangular(0.2556, 0.325, 0.3944, s))*SF_TW_SO TW_S_y16_data = (nr.triangular(0.2928, 0.35, 0.4072, s))*SF_TW_SO TW_S_y17_data = (nr.triangular(0.3133, 0.405, 0.4967, s))*SF_TW_SO TW_S_y18_data = (nr.triangular(0.3268, 0.46, 0.5932, s))*SF_TW_SO TW_S_y19_data = (nr.triangular(0.3850, 0.515, 0.6450, s))*SF_TW_SO TW_S_y20_data = (nr.triangular(0.4546, 0.57, 0.6854, s))*SF_TW_SO TW_S_y0 = cp.TransferDistribution(nr.choice, [TW_S_y0_data]) TW_S_y1 = cp.TransferDistribution(nr.choice, [TW_S_y1_data]) TW_S_y2 = cp.TransferDistribution(nr.choice, [TW_S_y2_data]) TW_S_y3 = cp.TransferDistribution(nr.choice, [TW_S_y3_data]) TW_S_y4 = cp.TransferDistribution(nr.choice, [TW_S_y4_data]) TW_S_y5 = cp.TransferDistribution(nr.choice, [TW_S_y5_data]) TW_S_y6 = cp.TransferDistribution(nr.choice, [TW_S_y6_data]) TW_S_y7 = cp.TransferDistribution(nr.choice, [TW_S_y7_data]) TW_S_y8 = cp.TransferDistribution(nr.choice, [TW_S_y8_data]) TW_S_y9 = cp.TransferDistribution(nr.choice, [TW_S_y9_data]) TW_S_y10 = cp.TransferDistribution(nr.choice, [TW_S_y10_data]) TW_S_y11 = cp.TransferDistribution(nr.choice, [TW_S_y11_data]) TW_S_y12 = cp.TransferDistribution(nr.choice, [TW_S_y12_data]) TW_S_y13 = cp.TransferDistribution(nr.choice, [TW_S_y13_data]) TW_S_y14 = cp.TransferDistribution(nr.choice, [TW_S_y14_data]) TW_S_y15 = cp.TransferDistribution(nr.choice, [TW_S_y15_data]) TW_S_y16 = cp.TransferDistribution(nr.choice, [TW_S_y16_data]) TW_S_y17 = cp.TransferDistribution(nr.choice, [TW_S_y17_data]) TW_S_y18 = cp.TransferDistribution(nr.choice, [TW_S_y18_data]) TW_S_y19 = cp.TransferDistribution(nr.choice, [TW_S_y19_data]) TW_S_y20 = cp.TransferDistribution(nr.choice, [TW_S_y20_data]) TextW.transfers = [cp.TimeDependendDistributionTransfer([TW_S_y0, TW_S_y1,TW_S_y2,TW_S_y3,TW_S_y4,TW_S_y5, TW_S_y6,TW_S_y7,TW_S_y8,TW_S_y9,TW_S_y10, TW_S_y11,TW_S_y12,TW_S_y13,TW_S_y14,TW_S_y15, TW_S_y16,TW_S_y17,TW_S_y18,TW_S_y19,TW_S_y20], Sorting_TextW, TextW, priority=2), cp.ConstTransfer(1, MMSW, priority=1)] # In[52]: ### SAME AS FOR EUROPE ### # CDW CDW_LF_y0 = cp.TransferDistribution(nr.triangular, [0.4028, 0.646, 0.8892]) CDW_LF_y1 = cp.TransferDistribution(nr.triangular, [0.4155, 0.617, 0.8185]) CDW_LF_y2 = cp.TransferDistribution(nr.triangular, [0.4254, 0.588, 0.7506]) CDW_LF_y3 = cp.TransferDistribution(nr.triangular, [0.4324, 0.559, 0.6856]) CDW_LF_y4 = cp.TransferDistribution(nr.triangular, [0.4447, 0.54, 0.6353]) CDW_LF_y5 = cp.TransferDistribution(nr.triangular, [0.3709, 0.507, 0.6431]) CDW_LF_y6 = cp.TransferDistribution(nr.triangular, [0.3032, 0.474, 0.6448]) CDW_LF_y7 = cp.TransferDistribution(nr.triangular, [0.2416, 0.441, 0.6404]) CDW_LF_y8 = cp.TransferDistribution(nr.triangular, [0.2268, 0.408, 0.5892]) CDW_LF_y9 = cp.TransferDistribution(nr.triangular, [0.2115, 0.375, 0.5385]) CDW_LF_y10 = cp.TransferDistribution(nr.triangular, [0.1956, 0.342, 0.4884]) CDW_LF_y11 = cp.TransferDistribution(nr.triangular, [0.1650, 0.316, 0.4670]) CDW_LF_y12 = cp.TransferDistribution(nr.triangular, [0.1665, 0.291, 0.4155]) CDW_LF_y13 = cp.TransferDistribution(nr.triangular, [0.1436, 0.275, 0.4065]) CDW_LF_y14 = cp.TransferDistribution(nr.triangular, [0.1481, 0.259, 0.3699]) CDW_LF_y15 = cp.TransferDistribution(nr.triangular, [0.1096, 0.21, 0.3104]) CDW_LF_y16 = cp.TransferDistribution(nr.triangular, [0.0905, 0.181, 0.2715]) CDW_LF_y17 = cp.TransferDistribution(nr.triangular, [0.076, 0.152, 0.228]) CDW_LF_y18 = cp.TransferDistribution(nr.triangular, [0.0615, 0.123, 0.1845]) CDW_LF_y19 = cp.TransferDistribution(nr.triangular, [0.0465, 0.093, 0.1395]) CDW_LF_y20 = cp.TransferDistribution(nr.triangular, [0.032, 0.0640, 0.0960]) CDW.transfers = [cp.TimeDependendDistributionTransfer([CDW_LF_y0, CDW_LF_y1,CDW_LF_y2,CDW_LF_y3,CDW_LF_y4,CDW_LF_y5, CDW_LF_y6,CDW_LF_y7,CDW_LF_y8,CDW_LF_y9,CDW_LF_y10, CDW_LF_y11,CDW_LF_y12,CDW_LF_y13,CDW_LF_y14,CDW_LF_y15, CDW_LF_y16,CDW_LF_y17,CDW_LF_y18,CDW_LF_y19,CDW_LF_y20], CDW_Landfill_P, CDW, priority=2), cp.ConstTransfer(1, Sorting_CDW, priority=1)] # In[53]: ############################ SORTING ######################################## # all these processes were kept static, so there are no time-dependent transfers # thus the coefficient of variance is kept to 50% # In[55]: ### SAME AS FOR EUROPE # Packaging Waste Sorting_PackW.transfers = [cp.ConstTransfer(1, WW_N, priority=1)] # WEEE Exports = np.concatenate([nr.triangular(0.055, 0.11, 0.11, int(0.134*s)), nr.uniform(0.11, 0.23, int(0.585*s)), nr.triangular(0.23, 0.23, 0.28, int(0.345*s))]) Sorting_WEEE.transfers = [cp.RandomChoiceTransfer(Exports, SortWEEE_Export_P, priority =2), cp.ConstTransfer(1, WEEE_NotExported, priority=1)] WEEE_NotExported.transfers = [cp.StochasticTransfer(nr.triangular, [0.015, 0.03, 0.045], Reprocess_WEEE_NotExW_P, priority =2), # representing light bulbs going all to reprocessing cp.ConstTransfer(1, OtherWEEE, priority =1)] OtherWEEE.transfers = [cp.StochasticTransfer(nr.triangular, [0.02, 0.04, 0.06], Reprocess_Metal_P, priority =2), #representing metals going all to reprocessing cp.StochasticTransfer(nr.triangular, [0.035, 0.07, 0.105], WEEE_Plastics, priority=2), cp.ConstTransfer(1, Sorting_Disposal, priority=1)] #representing other materials going all to WIP or landfill WEEE_Plastics.transfers = [cp.StochasticTransfer(nr.triangular, [0.035, 0.07, 0.105], Reprocess_Plas_P, priority=2), #0representing Mobile phones cp.ConstTransfer(1, WEEE_Plastics_OA, priority=1)] Resorting = np.concatenate([nr.triangular(0.15, 0.3, 0.3, int(0.273*s)), nr.uniform(0.3, 0.4, int(0.364*s)), nr.triangular(0.4, 0.4, 0.6, int(0.364*s))]) WEEE_Plastics_OA.transfers = [ cp.RandomChoiceTransfer(Resorting, WEEE_Resorting, priority=2), cp.ConstTransfer(1, WIP_P, priority=1)] WEEE_Resorting.transfers = [cp.StochasticTransfer(nr.triangular, [0.05, 0.1, 0.15], Sorting_Disposal, priority=2), cp.ConstTransfer(1, Reprocess_Plas_P, priority=1)] # Textile Waste (Clothing) TextW_Ex = np.concatenate([nr.triangular(0.15, 0.3, 0.3, int(0.273*s)), nr.uniform(0.3, 0.4, int(0.364*s)), nr.triangular(0.4, 0.4, 0.6, int(0.364*s))]) TextW_RU = np.concatenate([nr.triangular(0.015, 0.03, 0.03, int(0.073*s)), nr.uniform(0.03, 0.1, int(0.683*s)), nr.triangular(0.1, 0.1, 0.15, int(0.244*s))]) TextW_WIPLF = np.concatenate([nr.triangular(0.05, 0.1, 0.1, int(0.222*s)), nr.uniform(0.1, 0.15, int(0.444*s)), nr.triangular(0.15, 0.15, 0.225, int(0.333*s))]) Sorting_TextW.transfers = [cp.RandomChoiceTransfer(TextW_Ex, SortTextW_Export_P, priority=2), cp.RandomChoiceTransfer(TextW_RU, Textiles_Use, priority=2), cp.RandomChoiceTransfer(TextW_WIPLF, Sorting_Disposal, priority=2), # goes to incineration/landfilling in same shares as MMSW cp.ConstTransfer(1, TextW_Resorting, priority=1)] TextW_Rpr = np.concatenate([nr.triangular(0.2, 0.4, 0.4, int(0.222*s)), nr.uniform(0.4, 0.6, int(0.444*s)), nr.triangular(0.6, 0.6, 0.9, int(0.333*s))]) TextW_CR = np.concatenate([nr.triangular(0.25, 0.5, 0.5, int(0.333*s)), nr.uniform(0.5, 0.6, int(0.267*s)), nr.triangular(0.6, 0.6, 0.9, int(0.4*s))]) TextW_Resorting.transfers = [cp.RandomChoiceTransfer(TextW_Rpr, Reprocess_Text_P, priority=1), cp.RandomChoiceTransfer(TextW_CR, Sorting_Disposal, priority=1)] #Cleaning rags being disposed in same shares as MMSW # CDW ## The categories (materials) of CDW are distributed here. # For TiO2, we have: Total_mineral_CDW = nr.triangular(0, 0.0508, 0.10, s) + nr.triangular(0, 0.0286, 0.05, s) + nr.triangular(0, 0.0286, 0.05, s) Total_glass_CDW = nr.triangular(0, 0.0286, 0.05, s) Total_CDW = Total_mineral_CDW + Total_glass_CDW # Proportion of mineral and glass waste: Prop_Mine = Total_mineral_CDW/Total_CDW Prop_Glass = Total_glass_CDW/Total_CDW Sorting_CDW.transfers = [cp.StochasticTransfer(nr.choice, [Prop_Mine], Sorting_CDW_Mineral, priority=2), cp.ConstTransfer(1, Sorting_CDW_Glass, priority=1)] Glass_LF = np.concatenate([nr.uniform(0, 0.5, int(s*0.8)), nr.triangular(0.5, 0.5, 0.75, int(s*0.2))]) Sorting_CDW_Glass.transfers = [cp.RandomChoiceTransfer(Glass_LF, SortCDWGlass_Landfill_P, priority=2), cp.ConstTransfer(1, Reprocess_Glass_P, priority=1)] Reprocessing_Mineral = np.concatenate([nr.uniform( 0, 0.28, int(s*0.8)), nr.triangular(0.28, 0.28, 0.294, int(s*0.2))]) Sorting_CDW_Mineral.transfers= [cp.RandomChoiceTransfer(Reprocessing_Mineral, Reprocess_Mineral_P, priority=2), cp.ConstTransfer(1, SortCDWMiner_Landfill_P, priority=1)] # Sorting_Disposal, is a transfer to disposal with the same TCs as MMSW Sorting_Disposal.transfers = [cp.TimeDependendDistributionTransfer([MSW_LF_y0, MSW_LF_y1,MSW_LF_y2,MSW_LF_y3,MSW_LF_y4,MSW_LF_y5, MSW_LF_y6,MSW_LF_y7,MSW_LF_y8,MSW_LF_y9,MSW_LF_y10, MSW_LF_y11,MSW_LF_y12,MSW_LF_y13,MSW_LF_y14,MSW_LF_y15, MSW_LF_y16,MSW_LF_y17,MSW_LF_y18,MSW_LF_y19,MSW_LF_y20], SortDisp_Landfill_P, Sorting_Disposal, priority=2), cp.ConstTransfer(1, WIP_P, priority=1)] # In[56]: ###############################Transfers form technical compartments ########## ############################################################################### # In[57]: '''######################################################################################################################## ####---------------------------------------------Recycling System---------------------------------------------------#### ######################################################################################################################## ''' # all these processes were kept static, so there are no time-dependent transfers # thus the coefficient of variance is kept to 50% # In[58]: '''Reprocess_WEEE_NotExW_P''' # this represent the light bulbs, only the mercury containing light bulbs are considered to be recycled # after discussion with Bernd and Vero, we sepeculated that all ENMs are contained in the glass part Reprocess_WEEE_NotExW_P.transfers=[cp.ConstTransfer(1, Reprocess_Glass_P, priority=1)] # In[59]: '''Reprocess_Plas_P''' # Reprocessing Procedures Reprocess_Plas_P.transfers =[cp.StochasticTransfer(nr.triangular, [0.0001, 0.0002, 0.0003], Air_Plas, priority=2), cp.ConstTransfer(1, Washing_Plas, priority=1)] Washing_Plas.transfers=[cp.ConstTransfer(1, Melting_Plas, priority=1)] #flow to melting as product-embedded forms Melting_Plas.transfers=[cp.ConstTransfer(1, Granulation_Plas, priority=1)] #flow to melting as product-embedded forms #From granulation to landfill Granu_Plas_LF_y0 = cp.TransferDistribution(nr.triangular, [0.0471, 0.0790, 0.1109]) Granu_Plas_LF_y1 = cp.TransferDistribution(nr.triangular, [0.0471, 0.0790, 0.1109]) Granu_Plas_LF_y2 = cp.TransferDistribution(nr.triangular, [0.0471, 0.0790, 0.1109]) Granu_Plas_LF_y3 = cp.TransferDistribution(nr.triangular, [0.0471, 0.0790, 0.1109]) Granu_Plas_LF_y4 = cp.TransferDistribution(nr.triangular, [0.0471, 0.0790, 0.1109]) Granu_Plas_LF_y5 = cp.TransferDistribution(nr.triangular, [0.0592, 0.0790, 0.0988]) Granu_Plas_LF_y6 = cp.TransferDistribution(nr.triangular, [0.0592, 0.0790, 0.0988]) Granu_Plas_LF_y7 = cp.TransferDistribution(nr.triangular, [0.0592, 0.0790, 0.0988]) Granu_Plas_LF_y8 = cp.TransferDistribution(nr.triangular, [0.0592, 0.0790, 0.0988]) Granu_Plas_LF_y9 = cp.TransferDistribution(nr.triangular, [0.0592, 0.0790, 0.0988]) Granu_Plas_LF_y10 = cp.TransferDistribution(nr.triangular, [0.0610, 0.0790, 0.0970]) Granu_Plas_LF_y11 = cp.TransferDistribution(nr.triangular, [0.0628, 0.0814, 0.0999]) Granu_Plas_LF_y12 = cp.TransferDistribution(nr.triangular, [0.0646, 0.0837, 0.1028]) Granu_Plas_LF_y13 = cp.TransferDistribution(nr.triangular, [0.0625, 0.0810, 0.0994]) Granu_Plas_LF_y14 = cp.TransferDistribution(nr.triangular, [0.0604, 0.0782, 0.0961]) Granu_Plas_LF_y15 = cp.TransferDistribution(nr.triangular, [0.0619, 0.0802, 0.0985]) Granu_Plas_LF_y16 = cp.TransferDistribution(nr.triangular, [0.0634, 0.0822, 0.1009]) Granu_Plas_LF_y17 = cp.TransferDistribution(nr.triangular, [0.0616, 0.0822, 0.1028]) Granu_Plas_LF_y18 = cp.TransferDistribution(nr.triangular, [0.0616, 0.0822, 0.1028]) Granu_Plas_LF_y19 = cp.TransferDistribution(nr.triangular, [0.0616, 0.0822, 0.1028]) Granu_Plas_LF_y20 = cp.TransferDistribution(nr.triangular, [0.0616, 0.0822, 0.1028]) #From granulation to incineration Granu_Plas_WIP_y0 = cp.TransferDistribution(nr.triangular, [0.0057, 0.0096, 0.0134]) Granu_Plas_WIP_y1 = cp.TransferDistribution(nr.triangular, [0.0057, 0.0096, 0.0134]) Granu_Plas_WIP_y2 = cp.TransferDistribution(nr.triangular, [0.0057, 0.0096, 0.0134]) Granu_Plas_WIP_y3 = cp.TransferDistribution(nr.triangular, [0.0057, 0.0096, 0.0134]) Granu_Plas_WIP_y4 = cp.TransferDistribution(nr.triangular, [0.0057, 0.0096, 0.0134]) Granu_Plas_WIP_y5 = cp.TransferDistribution(nr.triangular, [0.0072, 0.0096, 0.0120]) Granu_Plas_WIP_y6 = cp.TransferDistribution(nr.triangular, [0.0072, 0.0096, 0.0120]) Granu_Plas_WIP_y7 = cp.TransferDistribution(nr.triangular, [0.0072, 0.0096, 0.0120]) Granu_Plas_WIP_y8 = cp.TransferDistribution(nr.triangular, [0.0072, 0.0096, 0.0120]) Granu_Plas_WIP_y9 = cp.TransferDistribution(nr.triangular, [0.0072, 0.0096, 0.0120]) Granu_Plas_WIP_y10 = cp.TransferDistribution(nr.triangular, [0.0074, 0.0096, 0.0118]) Granu_Plas_WIP_y11 = cp.TransferDistribution(nr.triangular, [0.0056, 0.0072, 0.0089]) Granu_Plas_WIP_y12 = cp.TransferDistribution(nr.triangular, [0.0038, 0.0049, 0.0060]) Granu_Plas_WIP_y13 = cp.TransferDistribution(nr.triangular, [0.0059, 0.0076, 0.0094]) Granu_Plas_WIP_y14 = cp.TransferDistribution(nr.triangular, [0.0080, 0.0104, 0.0127]) Granu_Plas_WIP_y15 = cp.TransferDistribution(nr.triangular, [0.0065, 0.0084, 0.0103]) Granu_Plas_WIP_y16 = cp.TransferDistribution(nr.triangular, [0.0050, 0.0064, 0.0079]) Granu_Plas_WIP_y17 = cp.TransferDistribution(nr.triangular, [0.0048, 0.0064, 0.0081]) Granu_Plas_WIP_y18 = cp.TransferDistribution(nr.triangular, [0.0048, 0.0064, 0.0081]) Granu_Plas_WIP_y19 = cp.TransferDistribution(nr.triangular, [0.0048, 0.0064, 0.0081]) Granu_Plas_WIP_y20 = cp.TransferDistribution(nr.triangular, [0.0048, 0.0064, 0.0081]) Granulation_Plas.transfers =[cp.TimeDependendDistributionTransfer([Granu_Plas_LF_y0, Granu_Plas_LF_y1,Granu_Plas_LF_y2,Granu_Plas_LF_y3,Granu_Plas_LF_y4,Granu_Plas_LF_y5, Granu_Plas_LF_y6,Granu_Plas_LF_y7,Granu_Plas_LF_y8,Granu_Plas_LF_y9,Granu_Plas_LF_y10, Granu_Plas_LF_y11,Granu_Plas_LF_y12,Granu_Plas_LF_y13,Granu_Plas_LF_y14,Granu_Plas_LF_y15, Granu_Plas_LF_y16,Granu_Plas_LF_y17,Granu_Plas_LF_y18,Granu_Plas_LF_y19,Granu_Plas_LF_y20], GranuPlas_Landfill_P, Granulation_Plas, priority=2), cp.TimeDependendDistributionTransfer([Granu_Plas_WIP_y0, Granu_Plas_WIP_y1,Granu_Plas_WIP_y2,Granu_Plas_WIP_y3,Granu_Plas_WIP_y4,Granu_Plas_WIP_y5, Granu_Plas_WIP_y6,Granu_Plas_WIP_y7,Granu_Plas_WIP_y8,Granu_Plas_WIP_y9,Granu_Plas_WIP_y10, Granu_Plas_WIP_y11,Granu_Plas_WIP_y12,Granu_Plas_WIP_y13,Granu_Plas_WIP_y14,Granu_Plas_WIP_y15, Granu_Plas_WIP_y16,Granu_Plas_WIP_y17,Granu_Plas_WIP_y18,Granu_Plas_WIP_y19,Granu_Plas_WIP_y20], RPWIP_P, Granulation_Plas, priority=2), cp.ConstTransfer(1, GranuPlas_Reuse_P, priority=1)] #Release form Air_Plas.transfers=[cp.StochasticTransfer(tf.TriangTrunc, [0.2, 1.07, 1, 0, 1], Air_NM_Plas,priority=2), cp.ConstTransfer(1, RPAir_M, priority=1)] Air_NM_Plas.transfers=[cp.StochasticTransfer(nr.triangular,[0.0959, 0.1535, 0.2211], RPAir_N, priority=2), cp.ConstTransfer(1, RPAir_M, priority=1)] # In[60]: '''Reprocess_Text_P''' # from shredding textiles to air Reprocess_Text_P.transfers=[cp.StochasticTransfer(tf.TriangTrunc, [0.02, 1.47, 1, 0, 1], Air_Text, priority=2), cp.ConstTransfer(1, Washing_Text, priority=1)] # from washing textiles to WW Washing_Text.transfers=[cp.StochasticTransfer(nr.triangular, [0.0006, 0.0009, 0.0012], WW_Text, priority=2), cp.ConstTransfer(1, Baling_Text, priority=1)] # from Baling to landfill Baling_Text_LF_y0 = cp.TransferDistribution(nr.triangular, [0.0384, 0.0644, 0.0904]) Baling_Text_LF_y1 = cp.TransferDistribution(nr.triangular, [0.0384, 0.0644, 0.0904]) Baling_Text_LF_y2 = cp.TransferDistribution(nr.triangular, [0.0384, 0.0644, 0.0904]) Baling_Text_LF_y3 = cp.TransferDistribution(nr.triangular, [0.0384, 0.0644, 0.0904]) Baling_Text_LF_y4 = cp.TransferDistribution(nr.triangular, [0.0384, 0.0644, 0.0904]) Baling_Text_LF_y5 = cp.TransferDistribution(nr.triangular, [0.0482, 0.0644, 0.0806]) Baling_Text_LF_y6 = cp.TransferDistribution(nr.triangular, [0.0482, 0.0644, 0.0806]) Baling_Text_LF_y7 = cp.TransferDistribution(nr.triangular, [0.0482, 0.0644, 0.0806]) Baling_Text_LF_y8 = cp.TransferDistribution(nr.triangular, [0.0482, 0.0644, 0.0806]) Baling_Text_LF_y9 = cp.TransferDistribution(nr.triangular, [0.0482, 0.0644, 0.0806]) Baling_Text_LF_y10 = cp.TransferDistribution(nr.triangular, [0.0497, 0.0644, 0.0791]) Baling_Text_LF_y11 = cp.TransferDistribution(nr.triangular, [0.0483, 0.0625, 0.0768]) Baling_Text_LF_y12 = cp.TransferDistribution(nr.triangular, [0.0468, 0.0607, 0.0745]) Baling_Text_LF_y13 = cp.TransferDistribution(nr.triangular, [0.0468, 0.0607, 0.0745]) Baling_Text_LF_y14 = cp.TransferDistribution(nr.triangular, [0.0468, 0.0607, 0.0745]) Baling_Text_LF_y15 = cp.TransferDistribution(nr.triangular, [0.0482, 0.0624, 0.0766]) Baling_Text_LF_y16 = cp.TransferDistribution(nr.triangular, [0.0495, 0.0642, 0.0788]) Baling_Text_LF_y17 = cp.TransferDistribution(nr.triangular, [0.0481, 0.0642, 0.0803]) Baling_Text_LF_y18 = cp.TransferDistribution(nr.triangular, [0.0481, 0.0642, 0.0803]) Baling_Text_LF_y19 = cp.TransferDistribution(nr.triangular, [0.0481, 0.0642, 0.0803]) Baling_Text_LF_y20 = cp.TransferDistribution(nr.triangular, [0.0481, 0.0642, 0.0803]) # from baling to WIP Baling_Text_WIP_y0 = cp.TransferDistribution(nr.triangular, [0.0033, 0.0056, 0.0079]) Baling_Text_WIP_y1 = cp.TransferDistribution(nr.triangular, [0.0033, 0.0056, 0.0079]) Baling_Text_WIP_y2 = cp.TransferDistribution(nr.triangular, [0.0033, 0.0056, 0.0079]) Baling_Text_WIP_y3 = cp.TransferDistribution(nr.triangular, [0.0033, 0.0056, 0.0079]) Baling_Text_WIP_y4 = cp.TransferDistribution(nr.triangular, [0.0033, 0.0056, 0.0079]) Baling_Text_WIP_y5 = cp.TransferDistribution(nr.triangular, [0.0042, 0.0056, 0.0070]) Baling_Text_WIP_y6 = cp.TransferDistribution(nr.triangular, [0.0042, 0.0056, 0.0070]) Baling_Text_WIP_y7 = cp.TransferDistribution(nr.triangular, [0.0042, 0.0056, 0.0070]) Baling_Text_WIP_y8 = cp.TransferDistribution(nr.triangular, [0.0042, 0.0056, 0.0070]) Baling_Text_WIP_y9 = cp.TransferDistribution(nr.triangular, [0.0042, 0.0056, 0.0070]) Baling_Text_WIP_y10 = cp.TransferDistribution(nr.triangular, [0.0043, 0.0056, 0.0069]) Baling_Text_WIP_y11 = cp.TransferDistribution(nr.triangular, [0.0058, 0.0075, 0.0092]) Baling_Text_WIP_y12 = cp.TransferDistribution(nr.triangular, [0.0072, 0.0093, 0.0115]) Baling_Text_WIP_y13 = cp.TransferDistribution(nr.triangular, [0.0072, 0.0093, 0.0115]) Baling_Text_WIP_y14 = cp.TransferDistribution(nr.triangular, [0.0072, 0.0093, 0.0115]) Baling_Text_WIP_y15 = cp.TransferDistribution(nr.triangular, [0.0059, 0.0076, 0.0093]) Baling_Text_WIP_y16 = cp.TransferDistribution(nr.triangular, [0.0045, 0.0058, 0.0072]) Baling_Text_WIP_y17 = cp.TransferDistribution(nr.triangular, [0.0044, 0.0058, 0.0073]) Baling_Text_WIP_y18 = cp.TransferDistribution(nr.triangular, [0.0044, 0.0058, 0.0073]) Baling_Text_WIP_y19 = cp.TransferDistribution(nr.triangular, [0.0044, 0.0058, 0.0073]) Baling_Text_WIP_y20 = cp.TransferDistribution(nr.triangular, [0.0044, 0.0058, 0.0073]) Baling_Text.transfers =[cp.TimeDependendDistributionTransfer([Baling_Text_LF_y0, Baling_Text_LF_y1,Baling_Text_LF_y2,Baling_Text_LF_y3,Baling_Text_LF_y4,Baling_Text_LF_y5, Baling_Text_LF_y6,Baling_Text_LF_y7,Baling_Text_LF_y8,Baling_Text_LF_y9,Baling_Text_LF_y10, Baling_Text_LF_y11,Baling_Text_LF_y12,Baling_Text_LF_y13,Baling_Text_LF_y14,Baling_Text_LF_y15, Baling_Text_LF_y16,Baling_Text_LF_y17,Baling_Text_LF_y18,Baling_Text_LF_y19,Baling_Text_LF_y20], BaliText_Landfill_P, Baling_Text, priority=2), cp.TimeDependendDistributionTransfer([Baling_Text_WIP_y0, Baling_Text_WIP_y1,Baling_Text_WIP_y2,Baling_Text_WIP_y3,Baling_Text_WIP_y4,Baling_Text_WIP_y5, Baling_Text_WIP_y6,Baling_Text_WIP_y7,Baling_Text_WIP_y8,Baling_Text_WIP_y9,Baling_Text_WIP_y10, Baling_Text_WIP_y11,Baling_Text_WIP_y12,Baling_Text_WIP_y13,Baling_Text_WIP_y14,Baling_Text_WIP_y15, Baling_Text_WIP_y16,Baling_Text_WIP_y17,Baling_Text_WIP_y18,Baling_Text_WIP_y19,Baling_Text_WIP_y20], RPWIP_P, Baling_Text, priority=2), cp.ConstTransfer(1, BaliText_Reuse_P, priority=1)] # Release form Air_Text.transfers=[cp.StochasticTransfer(tf.TriangTrunc, [0.2, 1.47, 1, 0, 1], Air_NM, priority=2), cp.ConstTransfer(1, RPAir_M, priority=1)] WW_Text.transfers=[cp.StochasticTransfer(nr.triangular,[0.2942, 0.46, 0.6258], RPWW_N, priority=2), cp.ConstTransfer(1, RPWW_M, priority=1)] # In[61]: '''Reprocess_Metal_P''' # from shredding metals to air Reprocess_Metal_P.transfers =[cp.StochasticTransfer(tf.TriangTrunc, [0.0002, 1.07, 1, 0, 1], Air_Metal, priority=2), cp.ConstTransfer(1, Melting_Metal, priority=1)] Melting_Metal.transfers=[cp.StochasticTransfer(nr.uniform, [0, 1], Purification_Metal_T, priority=2), cp.ConstTransfer(1, Purification_Metal_P, priority=1)] # from purification to Moulding and finally to Reuse # we assumed the TCs from purification to moulding is between 0.0001 to 0.01, #and thus consider 0.0001 as min and 0.01 as max to draw a trapezoidal distribution #3.05445180 is the calculated spread Purification_Metal_P.transfers =[cp.StochasticTransfer(tf.TrapezTrunc, [0.0001, 0.01, 3.05, 3.05, 1, 0, 1], PuriMetal_Reuse_P, priority=2), # since Tc from moulding to reuse is 1, here we skip moulding to reuse to simplify the code cp.ConstTransfer(1, Puri_slag_P, priority=1)] Purification_Metal_T.transfers =[cp.StochasticTransfer(tf.TrapezTrunc, [0.0001, 0.01, 3.05, 3.05, 1, 0, 1], PuriMetal_Reuse_T, priority=2), # since Tc from moulding to reuse is 1, here we skip moulding to reuse to simplify the code cp.ConstTransfer(1, Puri_slag_T, priority=1)] # from slags to landfill Puri_Metal_LF_y0 = cp.TransferDistribution(nr.triangular, [0.2791, 0.35, 0.4209]) Puri_Metal_LF_y1 = cp.TransferDistribution(nr.triangular, [0.2894, 0.35, 0.4106]) Puri_Metal_LF_y2 = cp.TransferDistribution(nr.triangular, [0.2791, 0.35, 0.4209]) Puri_Metal_LF_y3 = cp.TransferDistribution(nr.triangular, [0.2791, 0.35, 0.4209]) Puri_Metal_LF_y4 = cp.TransferDistribution(nr.triangular, [0.2791, 0.35, 0.4209]) Puri_Metal_LF_y5 = cp.TransferDistribution(nr.triangular, [0.2791, 0.35, 0.4209]) Puri_Metal_LF_y6 = cp.TransferDistribution(nr.triangular, [0.2791, 0.35, 0.4209]) Puri_Metal_LF_y7 = cp.TransferDistribution(nr.triangular, [0.2186, 0.35, 0.4814]) Puri_Metal_LF_y8 = cp.TransferDistribution(nr.triangular, [0.2186, 0.35, 0.4814]) Puri_Metal_LF_y9 = cp.TransferDistribution(nr.triangular, [0.2186, 0.35, 0.4814]) Puri_Metal_LF_y10 = cp.TransferDistribution(nr.triangular, [0.2186, 0.35, 0.4814]) Puri_Metal_LF_y11 = cp.TransferDistribution(nr.triangular, [0.2186, 0.35, 0.4814]) Puri_Metal_LF_y12 = cp.TransferDistribution(tf.TriangTrunc, [0.35, 1.03, 1, 0, 1]) Puri_Metal_LF_y13 = cp.TransferDistribution(tf.TriangTrunc, [0.35, 1.03, 1, 0, 1]) Puri_Metal_LF_y14 = cp.TransferDistribution(tf.TriangTrunc, [0.35, 1.03, 1, 0, 1]) Puri_Metal_LF_y15 = cp.TransferDistribution(tf.TriangTrunc, [0.35, 1.03, 1, 0, 1]) Puri_Metal_LF_y16 = cp.TransferDistribution(tf.TriangTrunc, [0.35, 1.03, 1, 0, 1]) Puri_Metal_LF_y17 = cp.TransferDistribution(tf.TriangTrunc, [0.35, 1.03, 1, 0, 1]) Puri_Metal_LF_y18 = cp.TransferDistribution(tf.TriangTrunc, [0.35, 1.03, 1, 0, 1]) Puri_Metal_LF_y19 = cp.TransferDistribution(tf.TriangTrunc, [0.35, 1.03, 1, 0, 1]) Puri_Metal_LF_y20 = cp.TransferDistribution(tf.TriangTrunc, [0.35, 1.03, 1, 0, 1]) Puri_slag_T.transfers =[cp.TimeDependendDistributionTransfer([Puri_Metal_LF_y0, Puri_Metal_LF_y1,Puri_Metal_LF_y2,Puri_Metal_LF_y3,Puri_Metal_LF_y4,Puri_Metal_LF_y5, Puri_Metal_LF_y6,Puri_Metal_LF_y7,Puri_Metal_LF_y8,Puri_Metal_LF_y9,Puri_Metal_LF_y10, Puri_Metal_LF_y11,Puri_Metal_LF_y12,Puri_Metal_LF_y13,Puri_Metal_LF_y14,Puri_Metal_LF_y15, Puri_Metal_LF_y16,Puri_Metal_LF_y17,Puri_Metal_LF_y18,Puri_Metal_LF_y19,Puri_Metal_LF_y20], PuriSlag_Landfill_T, Puri_slag_T, priority=2), cp.ConstTransfer(1, PuriSlag_Reuse_T, priority=1)] Puri_slag_P.transfers =[cp.TimeDependendDistributionTransfer([Puri_Metal_LF_y0, Puri_Metal_LF_y1,Puri_Metal_LF_y2,Puri_Metal_LF_y3,Puri_Metal_LF_y4,Puri_Metal_LF_y5, Puri_Metal_LF_y6,Puri_Metal_LF_y7,Puri_Metal_LF_y8,Puri_Metal_LF_y9,Puri_Metal_LF_y10, Puri_Metal_LF_y11,Puri_Metal_LF_y12,Puri_Metal_LF_y13,Puri_Metal_LF_y14,Puri_Metal_LF_y15, Puri_Metal_LF_y16,Puri_Metal_LF_y17,Puri_Metal_LF_y18,Puri_Metal_LF_y19,Puri_Metal_LF_y20], PuriSlag_Landfill_P, Puri_slag_P, priority=2), cp.ConstTransfer(1, PuriSlag_Reuse_P, priority=1)] # release form Air_Metal.transfers=[cp.StochasticTransfer(tf.TriangTrunc, [0.2, 1.47, 1, 0, 1], Air_NM, priority=2), cp.ConstTransfer(1, RPAir_M, priority=1)] # In[62]: '''Reprocess_Glass_P''' # from crushing galss to air Reprocess_Glass_P.transfers =[cp.StochasticTransfer(tf.TrapezTrunc, [0.07, 0.47, 0.49, 0.49, 1, 0, 1], Air_Glass, priority=2), # since Tc from moulding to reuse is 1, here we skip moulding to reuse to simplify the code cp.ConstTransfer(1, Melting_Glass, priority=1)] Melting_Glass.transfers=[cp.ConstTransfer(1, MeltGlass_Reuse_T, priority=1)] # Release form Air_Glass.transfers=[cp.StochasticTransfer(tf.TriangTrunc, [0.2, 1.47, 1, 0, 1], Air_NM,priority=2), cp.ConstTransfer(1, RPAir_M, priority=1)] # In[63]: '''#7. Reprocess_Mineral_P''' # from crushing minerals to air Reprocess_Mineral_P.transfers =[cp.StochasticTransfer(tf.TrapezTrunc, [0.05, 0.16, 0.48, 0.48, 1, 0, 1], Air_Mineral, priority=2), # since Tc from moulding to reuse is 1, here we skip moulding to reuse to simplify the code cp.ConstTransfer(1, Sorting_Mineral, priority=1)] Sorting_Mineral.transfers =[cp.StochasticTransfer(tf.TriangTrunc, [0.7944, 1.02, 1, 0, 1], SortMiner_Landfill_P, priority=2), cp.StochasticTransfer(tf.TriangTrunc, [0.05, 1.01, 1, 0, 1], Air_Mineral, priority=2), cp.ConstTransfer(1, SortMiner_Reuse_P, priority=1)] # release form Air_Mineral.transfers=[cp.StochasticTransfer(tf.TriangTrunc, [0.2, 1.47, 1, 0, 1], Air_NM, priority=2), cp.ConstTransfer(1, RPAir_M, priority=1)] # In[62]: # Allocate all release to air with mixed N and M # allocation of Air_NM for the fraction of 0.2 Air_NM.transfers=[cp.StochasticTransfer(tf.TriangTrunc, [0.1535, 1.07, 1, 0, 1], RPAir_N, priority=2), cp.ConstTransfer(1, RPAir_M,priority=1)] # In[64]: RPWIP_P.transfers=[cp.ConstTransfer(1, WIP_P, priority=1)] RPWW_N.transfers=[cp.ConstTransfer(1, WW_N, priority=1)] RPWW_M.transfers=[cp.ConstTransfer(1, WW_M, priority=1)] # In[65]: ################################## WIP ###################################### # In[66]: # Elimination while burning WIP_P.transfers =[cp.StochasticTransfer(nr.uniform, [1e-6, 1], Burning_T, priority=2), cp.ConstTransfer(1, Burning_N, priority=1)] WIP_M.transfers =[cp.StochasticTransfer(nr.uniform,[1e-6 , 1], Burning_T, priority=2), cp.ConstTransfer(1, Burning_N, priority=1)] WIP_N.transfers =[cp.StochasticTransfer(nr.uniform,[1e-6 , 1], Burning_T, priority=2), cp.ConstTransfer(1, Burning_N, priority=1)] WIP_T.transfers = [cp.ConstTransfer(1, Burning_T, priority=1)] # transfer to bottom / fly ash ToFlyash = np.concatenate([ nr.triangular(0.095, 0.19, 0.285, s), nr.triangular(0.0276, 0.029, 0.029, int(0.1759*s)), nr.uniform(0.029, 0.197, int(0.63241*s)), nr.triangular(0.197, 0.197, 0.2069, int(0.1917*s)) ]) Burning_N.transfers =[cp.RandomChoiceTransfer(ToFlyash, Flyash_N, priority=2), cp.ConstTransfer(1, Bottomash_N, priority=1)] Burning_T.transfers =[cp.RandomChoiceTransfer(ToFlyash, Flyash_T, priority=2), cp.ConstTransfer(1, Bottomash_T, priority=1)] Flyash_N.transfers =[cp.StochasticTransfer(nr.triangular, [0.000025, 0.00005, 0.000075], WIP_Air_N, priority=2), # 1- 0.99995 (filter efficiency) cp.StochasticTransfer(nr.triangular, [0, 0.0001, 0.00015], WW_N, priority=2), # flow to Wastewater during the acid washing cp.ConstTransfer(1, Filterash_N, priority=1)] Flyash_T.transfers =[cp.StochasticTransfer(nr.triangular, [0.000025, 0.00005, 0.000075], WIP_Air_T, priority=2), # 1- 0.99995 (filter efficiency) cp.StochasticTransfer(nr.triangular, [0, 0.0001, 0.00015], WW_T, priority=2), # flow to Wastewater during the acid washing cp.ConstTransfer(1, Filterash_T, priority=1)] # In[67]: # From filter ash FA_LF_y0 = cp.TransferConstant(1) FA_LF_y1 = cp.TransferConstant(1) FA_LF_y2 = cp.TransferConstant(1) FA_LF_y3 = cp.TransferConstant(1) FA_LF_y4 = cp.TransferConstant(1) FA_LF_y5 = cp.TransferConstant(1) FA_LF_y6 = cp.TransferConstant(1) FA_LF_y7 = cp.TransferConstant(1) FA_LF_y8 = cp.TransferConstant(1) FA_LF_y9 = cp.TransferConstant(1) FA_LF_y10 = cp.TransferConstant(1) FA_LF_y11 = cp.TransferConstant(1) FA_LF_y12 = cp.TransferConstant(1) FA_LF_y13 = cp.TransferConstant(1) FA_LF_y14 = cp.TransferConstant(1) FA_LF_y15 = cp.TransferConstant(1) FA_LF_y16 = cp.TransferConstant(1) FA_LF_y17 = cp.TransferConstant(1) FA_LF_y18 = cp.TransferConstant(1) FA_LF_y19 = cp.TransferConstant(1) FA_LF_y20 = cp.TransferConstant(1) Filterash_N.transfers =[cp.TimeDependendDistributionTransfer([FA_LF_y0, FA_LF_y1,FA_LF_y2,FA_LF_y3,FA_LF_y4,FA_LF_y5, FA_LF_y6,FA_LF_y7,FA_LF_y8,FA_LF_y9,FA_LF_y10, FA_LF_y11,FA_LF_y12,FA_LF_y13,FA_LF_y14,FA_LF_y15, FA_LF_y16,FA_LF_y17,FA_LF_y18,FA_LF_y19,FA_LF_y20], Filterash_Landfill_N, Filterash_N, priority=2), cp.ConstTransfer(1, Filterash_Reuse_N, priority=1)] Filterash_T.transfers =[cp.TimeDependendDistributionTransfer([FA_LF_y0, FA_LF_y1,FA_LF_y2,FA_LF_y3,FA_LF_y4,FA_LF_y5, FA_LF_y6,FA_LF_y7,FA_LF_y8,FA_LF_y9,FA_LF_y10, FA_LF_y11,FA_LF_y12,FA_LF_y13,FA_LF_y14,FA_LF_y15, FA_LF_y16,FA_LF_y17,FA_LF_y18,FA_LF_y19,FA_LF_y20], Filterash_Landfill_T, Filterash_T, priority=2), cp.ConstTransfer(1, Filterash_Reuse_T, priority=1)] # In[68]: ### SAME AS EUROPE ### # From bottom ash BA_R_y0 = cp.TransferDistribution(nr.triangular, [0.2445, 0.489, 0.7335]) BA_R_y1 = cp.TransferDistribution(nr.triangular, [0.243, 0.486, 0.729]) BA_R_y2 = cp.TransferDistribution(nr.triangular, [0.2470, 0.482, 0.7170]) BA_R_y3 = cp.TransferDistribution(nr.triangular, [0.2531, 0.45, 0.6469]) BA_R_y4 = cp.TransferDistribution(nr.triangular, [0.2531, 0.45, 0.6469]) BA_R_y5 = cp.TransferDistribution(nr.triangular, [0.2365, 0.46, 0.6835]) BA_R_y6 = cp.TransferDistribution(nr.triangular, [0.235, 0.47, 0.705]) BA_R_y7 = cp.TransferDistribution(nr.triangular, [0.24, 0.48, 0.72]) BA_R_y8 = cp.TransferDistribution(nr.triangular, [0.245, 0.49, 0.735]) BA_R_y9 = cp.TransferDistribution(nr.triangular, [0.2602, 0.5, 0.7398]) BA_R_y10 = cp.TransferDistribution(nr.triangular, [0.2917, 0.51, 0.7283]) BA_R_y11 = cp.TransferDistribution(nr.triangular, [0.2401, 0.46, 0.6799]) BA_R_y12 = cp.TransferDistribution(nr.triangular, [0.2345, 0.41, 0.5855]) BA_R_y13 = cp.TransferDistribution(nr.triangular, [0.2156, 0.413, 0.6104]) BA_R_y14 = cp.TransferDistribution(nr.triangular, [0.2380, 0.416, 0.5940]) BA_R_y15 = cp.TransferDistribution(nr.triangular, [0.2276, 0.436, 0.6444]) BA_R_y16 = cp.TransferDistribution(nr.triangular, [0.2165, 0.433, 0.6495]) BA_R_y17 = cp.TransferDistribution(nr.triangular, [0.2145, 0.429, 0.6435]) BA_R_y18 = cp.TransferDistribution(nr.triangular, [0.2125, 0.425, 0.6375]) BA_R_y19 = cp.TransferDistribution(nr.triangular, [0.211, 0.422, 0.633]) BA_R_y20 = cp.TransferDistribution(nr.triangular, [0.209, 0.418, 0.627]) Bottomash_N.transfers =[cp.TimeDependendDistributionTransfer([BA_R_y0, BA_R_y1,BA_R_y2,BA_R_y3,BA_R_y4,BA_R_y5, BA_R_y6,BA_R_y7,BA_R_y8,BA_R_y9,BA_R_y10, BA_R_y11,BA_R_y12,BA_R_y13,BA_R_y14,BA_R_y15, BA_R_y16,BA_R_y17,BA_R_y18,BA_R_y19,BA_R_y20], Bottomash_Reuse_N, Bottomash_N, priority=2), cp.ConstTransfer(1, Bottomash_Landfill_N, priority=1)] Bottomash_T.transfers =[cp.TimeDependendDistributionTransfer([BA_R_y0, BA_R_y1,BA_R_y2,BA_R_y3,BA_R_y4,BA_R_y5, BA_R_y6,BA_R_y7,BA_R_y8,BA_R_y9,BA_R_y10, BA_R_y11,BA_R_y12,BA_R_y13,BA_R_y14,BA_R_y15, BA_R_y16,BA_R_y17,BA_R_y18,BA_R_y19,BA_R_y20], Bottomash_Reuse_T, Bottomash_T, priority=2), cp.ConstTransfer(1, Bottomash_Landfill_T, priority=1)] # In[ ]: