From 73593744b71c0e3c974e76357c7b3be7e73f13f3 Mon Sep 17 00:00:00 2001 From: Peter Yefi Date: Thu, 17 Mar 2022 10:37:47 -0400 Subject: [PATCH] Modified coefficient computation for water to water hp --- .../energy_systems/water_to_water_hp.py | 19 - .../heat_pumps/w2w_parallel.txt | 1620 ++++++++-------- data/energy_systems/heat_pumps/w2w_series.txt | 1705 ++++++++--------- exports/energy_systems/heat_pump_export.py | 12 +- .../water_to_water_hp_export.py | 18 +- exports/energy_systems_factory.py | 2 +- .../water_to_water_hp_parameters.py | 10 +- 7 files changed, 1707 insertions(+), 1679 deletions(-) diff --git a/city_model_structure/energy_systems/water_to_water_hp.py b/city_model_structure/energy_systems/water_to_water_hp.py index ab81623f..e65d2247 100644 --- a/city_model_structure/energy_systems/water_to_water_hp.py +++ b/city_model_structure/energy_systems/water_to_water_hp.py @@ -20,7 +20,6 @@ class WaterToWaterHP(HeatPump): self._total_cooling_capacity = None self._power_demand = None self._flow_rate = None - self._heat_output_coff = None # b coefficients self._power_demand_coff = None # a coefficients @property @@ -112,24 +111,6 @@ class WaterToWaterHP(HeatPump): if self._flow_rate is None: self._flow_rate = value - @property - def heat_output_coff(self) -> List[float]: - """ - Get coefficients for total cooling capacity - :return: [[float]] - """ - return self._heat_output_coff - - @heat_output_coff.setter - def heat_output_coff(self, value): - """ - Set coefficients for totol cooling capacity - :param value: [[float]] - :return: - """ - if self._heat_output_coff is None: - self._heat_output_coff = value - @property def power_demand_coff(self) -> List[float]: """ diff --git a/data/energy_systems/heat_pumps/w2w_parallel.txt b/data/energy_systems/heat_pumps/w2w_parallel.txt index 1f1b6120..44bfb9ae 100644 --- a/data/energy_systems/heat_pumps/w2w_parallel.txt +++ b/data/energy_systems/heat_pumps/w2w_parallel.txt @@ -1,482 +1,429 @@ -B 1 DIV - 256.1 - 135.1 - -B 2 DIV - 86.1 - 86.2 - -B 3 DIV - 1.1 - 279.1 - -B 4 DIV - 257.1 - 119.1 - -B 5 DIV - 123.1 - 168.1 - -B 6 DIV - 179.1 - 127.1 - -B 7 DIV - 123.1 - 174.1 - -B 8 DIV - 179.1 - 167.1 - -B 9 DIV - 107.1 - 123.1 - -B 10 DIV - 171.1 - 147.1 - -B 11 DIV - 264.1 - 146.1 - -B 76 SCREEN - 2.1 -P 76 - '*' % Format - 'HP Seasonal COP' % Headline - -B 77 SCREEN - 87.1 -P 77 - '*' % Format - 'Total Cost of the Auxiliary Heater Fuel in a Year (CAD)' % Headline - -B 78 SCREEN - 85.1 -P 78 - '*' % Format - 'Total Electricity Cost in a Year (CAD)' % Headline - -B 79 SCREEN - 84.1 -P 79 - '*' % Format - 'Total Fossil Fuel consumption in a Year (m3)' % Headline - -B 80 SCREEN - 89.1 -P 80 - '*' % Format - 'Total CO2 Emissions from Electricity Grid (g)' % Headline - -B 81 SCREEN - 83.1 -P 81 - '*' % Format - 'Total CO2 Emissions from Auxiliary Heater (g)' % Headline - -B 82 SCREEN - 88.1 -P 82 - '*' % Format - 'Total Electricuty Demand of Heat Pumps in a year (kWh)' % Headline - -B 83 CUM - 140.1 - -B 84 CUM - 255.1 - -B 85 CUM - 110.1 - -B 86 CUM - 107.1 - 123.1 - -B 87 CUM - 104.1 - -B 88 CUM - 259.1 - -B 89 CUM - 120.1 - -B 90 SUM - 141.1 - 108.1 - 105.1 - 114.1 - 125.1 - 131.1 - 134.1 - 112.1 - 130.1 - 111.1 - 201.1 - -B 91 SUM - 219.1 - 268.1 - -B 92 SUM - 132.1 - 143.1 - -B 93 SUM - 148.1 - 269.2 - -B 94 SUM - 187.1 - 220.1 - -B 95 SUM - 253.1 - 212.1 - -B 96 SUM - 117.1 - 116.1 - 136.1 - 133.1 - 124.1 - 122.1 - 118.1 - 142.1 - 106.1 - 121.1 - 150.1 - -B 97 SUM - 164.1 - 218.1 - -B 98 SUM - 113.1 - 274.1 - -B 99 SUM - 268.5 - 221.1 - -B 100 SUM - 103.1 - 155.1 - -B 101 SUM - 222.1 - 215.1 - -B 102 SUM - 109.1 - 129.1 - 137.1 - -B 103 MUL - 4.1 - 199.1 - -B 104 MUL - 265.1 - 151.1 - 209.1 - -B 105 MUL - 250.1 +B 1 MUL 149.1 + 264.1 -B 106 MUL - 187.1 - 102.1 - 178.1 - -B 107 MUL - 279.1 - 96.1 - 95.1 - 98.1 - -B 108 MUL - 280.1 +B 2 MUL + 259.1 190.1 -B 109 MUL +B 3 MUL + 296.1 160.1 + 289.2 + 298.1 + 83.1 + +B 4 MUL 252.1 + 79.1 + 151.1 -B 110 MUL - 7.1 - 175.1 - -B 111 MUL - 280.1 - 187.1 - 102.1 - 188.1 - -B 112 MUL - 280.1 - 102.1 - 145.1 - -B 113 MUL - 273.1 - 138.1 - -B 114 MUL - 187.1 - 144.1 - -B 115 MUL - 100.1 - 177.1 - 97.1 - 251.1 - -B 116 MUL - 280.1 - 172.1 - -B 117 MUL - 191.1 - 245.1 - -B 118 MUL - 280.1 - 187.1 - 184.1 - -B 119 MUL - 208.1 - 154.1 - -B 120 MUL - 5.1 +B 5 MUL + 79.1 205.1 -B 121 MUL - 280.1 - 187.1 - 102.1 +B 6 MUL + 33.1 + 178.1 + +B 7 MUL + 162.1 + 260.1 + +B 8 MUL + 171.1 + 297.1 + +B 9 MUL 176.1 + 79.1 + 169.1 -B 122 MUL - 102.1 - 165.1 - -B 123 MUL - 279.1 - 90.1 - 95.1 - 98.1 - -B 124 MUL - 166.1 - 249.1 - -B 125 MUL - 189.1 - 248.1 - -B 126 MUL - 102.1 - 279.1 - -B 127 MUL - 139.1 - 156.1 - -B 128 MUL - 213.1 - 126.1 - -B 129 MUL - 3.1 - 242.1 - 254.1 - -B 130 MUL - 187.1 - 102.1 - 207.1 - -B 131 MUL - 102.1 - 206.1 - -B 132 MUL - 91.1 - 210.1 - -B 133 MUL - 187.1 - 183.1 - -B 134 MUL - 280.1 - 187.1 - 159.1 - -B 135 MUL - 181.1 - 94.1 - -B 136 MUL - 247.1 - 196.1 - -B 137 MUL - 163.1 - 214.1 - -B 138 MUL - 273.1 - 243.1 - 211.1 - -B 139 MUL - 8.1 - 195.1 - -B 140 MUL - 262.1 - 158.1 - 161.1 - -B 141 MUL - 182.1 - 246.1 - -B 142 MUL - 280.1 - 102.1 +B 10 MUL + 176.1 + 79.1 170.1 -B 143 MUL - 275.1 - 185.1 +B 11 MUL + 255.1 + 40.1 + +B 12 MUL + 176.1 + 164.1 + +B 13 MUL + 304.1 + 231.1 + 256.1 + +B 14 MUL + 245.1 + 233.1 + 299.1 + +B 15 MUL + 282.1 + 153.1 + +B 16 MUL + 79.1 + 192.1 + +B 17 MUL + 252.1 + 176.1 + 79.1 + 195.1 + +B 18 MUL + 304.1 + 13.1 + +B 19 MUL + 281.1 + 89.1 + 84.1 + 81.1 + +B 20 MUL + 202.1 + 258.1 + +B 21 MUL + 261.1 + 199.1 + +B 22 MUL + 146.1 + 263.1 + +B 23 MUL + 252.1 + 79.1 + 186.1 + +B 24 MUL + 252.1 + 189.1 + +B 25 MUL + 242.1 + 175.1 + +B 26 MUL + 283.1 + 193.1 + 148.1 + +B 27 MUL + 252.1 + 155.1 + +B 28 MUL + 289.1 + 257.1 + +B 29 MUL + 252.1 + 176.1 + 79.1 + 163.1 + +B 30 MUL + 253.1 + 88.1 + +B 31 MUL + 240.1 + 161.1 + +B 32 MUL + 286.1 + 145.1 + 144.1 + +B 33 MUL + 241.1 + 200.1 + +B 34 MUL + 194.1 + 262.1 + +B 35 MUL + 252.1 + 176.1 + 201.1 + +B 36 MUL + 172.1 + 76.1 + +B 37 MUL + 252.1 + 176.1 + 198.1 + +B 38 MUL + 295.1 + 86.1 + +B 39 MUL + 302.1 + 157.1 + +B 40 MUL + 79.1 + 281.1 + +B 41 MUL + 281.1 + 85.1 + 84.1 + 81.1 + +B 42 MUL + 176.1 + 166.1 + +B 61 NOP + 72.1 + +B 62 NOP + 30.1 + +B 63 NOP + 80.1 + +B 64 NOP + 38.1 + +B 65 NOP + 71.1 + +B 66 NOP + 80.1 + +B 67 NOP + 68.1 + +B 68 NOP + 289.2 + +B 69 NOP + 11.1 + +B 70 NOP + 279.1 + +B 71 NOP + 176.1 + +B 72 NOP + 71.1 + +B 73 NOP + 287.1 + +B 74 NOP + 38.1 + +B 75 SUM + 28.1 + 39.1 + +B 76 SUM + 176.1 + 229.1 + +B 77 SUM + 224.1 + 277.1 + +B 78 SUM + 197.1 + 280.2 + +B 79 SUM + 8.1 + 14.1 + 20.1 + +B 80 SUM + 75.1 + 227.1 + +B 81 SUM + 18.1 + 303.1 + +B 82 SUM + 67.1 + 226.1 + +B 83 SUM + 177.1 + 223.1 + +B 84 SUM + 294.1 + 254.1 + +B 85 SUM + 1.1 + 27.1 + 21.1 + 42.1 + 34.1 + 16.1 + 37.1 + 4.1 + 10.1 + 29.1 + 196.1 + +B 86 SUM + 230.1 + 234.1 + +B 87 SUM + 69.1 + 225.1 + +B 88 SUM + 69.1 + 228.1 + +B 89 SUM + 22.1 + 24.1 + 2.1 + 12.1 + 7.1 + 5.1 + 35.1 + 23.1 + 9.1 + 17.1 + 179.1 B 144 CONST P 144 - $b4 % Constant value + $FuelEF % Constant value B 145 CONST P 145 - $b8 % Constant value + 300 % Constant value B 146 CONST P 146 - $FuelLHV % Constant value + $b1 % Constant value B 147 CONST P 147 - $HPNominalCapacity % Constant value + $HPDisactivationTemperature % Constant value B 148 CONST P 148 - 10 % Constant value + $FuelPrice % Constant value B 149 CONST P 149 - $b3 % Constant value + $a1 % Constant value B 150 CONST P 150 - $a11 % Constant value + $FuelLHV % Constant value B 151 CONST P 151 - 300 % Constant value + $a8 % Constant value B 152 CONST P 152 - 40 % Constant value + $HPNominalCapacity % Constant value B 153 CONST P 153 - 2 % Constant value + $TemperatureDifference % Constant value B 154 CONST P 154 - $TemperatureDifference % Constant value + 2 % Constant value B 155 CONST P 155 - 0 % Constant value + $a2 % Constant value B 156 CONST P 156 - $Cp % Constant value + 2 % Constant value B 157 CONST P 157 - 2 % Constant value + $BuildingSuppTemp % Constant value B 158 CONST P 158 - 300 % Constant value + $AuxHeaterEfficiency % Constant value B 159 CONST P 159 - $b7 % Constant value + 5 % Constant value B 160 CONST P 160 - $LowestPossibleLoadFlow % Constant value + $Cp % Constant value B 161 CONST P 161 - $FuelEF % Constant value + $ElecGridEF % Constant value B 162 CONST P 162 - $TemperatureDifference % Constant value + $b5 % Constant value B 163 CONST P 163 - $HighestPossibleLoadFlow % Constant value + $a10 % Constant value B 164 CONST P 164 - $BuildingSuppTemp % Constant value + $b4 % Constant value B 165 CONST P 165 - $a6 % Constant value + $HPReactivationTemperature % Constant value B 166 CONST P 166 - $a5 % Constant value + $a4 % Constant value B 167 CONST P 167 - $HPNominalCapacity % Constant value + 2 % Constant value B 168 CONST P 168 - 12 % Constant value + 9 % Constant value B 169 CONST P 169 - 5 % Constant value + $b9 % Constant value B 170 CONST P 170 - $a8 % Constant value + $a9 % Constant value B 171 CONST P 171 - $MaximumHPEnergyInput % Constant value + $LowestPossibleLoadFlow % Constant value B 172 CONST P 172 - $a2 % Constant value + $Cp % Constant value B 173 CONST P 173 @@ -484,7 +431,7 @@ P 173 B 174 CONST P 174 - 12 % Constant value + $HPNominalCapacity % Constant value B 175 CONST P 175 @@ -492,487 +439,543 @@ P 175 B 176 CONST P 176 - $a10 % Constant value + $HPSupTemp % Constant value B 177 CONST P 177 - $Cp % Constant value + $BuildingSuppTemp % Constant value B 178 CONST P 178 - $a9 % Constant value + $Cp % Constant value B 179 CONST P 179 - $MaximumHPEnergyInput % Constant value + $b11 % Constant value B 180 CONST P 180 - 0 % Constant value + 40 % Constant value B 181 CONST P 181 - $Cp % Constant value + 0 % Constant value B 182 CONST P 182 - $b1 % Constant value + $Cp % Constant value B 183 CONST P 183 - $a4 % Constant value + $TemperatureDifference % Constant value B 184 CONST P 184 - $a7 % Constant value + 2 % Constant value B 185 CONST P 185 - 25 % Constant value + 0 % Constant value B 186 CONST P 186 - 2 % Constant value + $b8 % Constant value B 187 CONST P 187 - $HPSupTemp % Constant value + $MaximumHPEnergyInput % Constant value B 188 CONST P 188 - $b10 % Constant value + 12 % Constant value B 189 CONST P 189 - $b5 % Constant value + $b2 % Constant value B 190 CONST P 190 - $b2 % Constant value + $b3 % Constant value B 191 CONST P 191 - $a1 % Constant value + 2 % Constant value B 192 CONST P 192 - $MaximumHPEnergyInput % Constant value + $a6 % Constant value B 193 CONST P 193 - 9 % Constant value + 300 % Constant value B 194 CONST P 194 - 0 % Constant value + $a5 % Constant value B 195 CONST P 195 - 7.13 % Constant value + $b10 % Constant value B 196 CONST P 196 - $a3 % Constant value + $a11 % Constant value B 197 CONST P 197 - 2 % Constant value + 10 % Constant value B 198 CONST P 198 - $HPReactivationTemperature % Constant value + $a7 % Constant value B 199 CONST P 199 - 1 % Constant value + $a3 % Constant value B 200 CONST P 200 - 2 % Constant value + 7.13 % Constant value B 201 CONST P 201 - $b11 % Constant value + $b7 % Constant value B 202 CONST P 202 - $HPDisactivationTemperature % Constant value + $HighestPossibleLoadFlow % Constant value B 203 CONST P 203 - 40 % Constant value + 12 % Constant value B 204 CONST P 204 - 0.40 % Constant value + $MaximumHPEnergyInput % Constant value B 205 CONST P 205 - $ElecGridEF % Constant value + $b6 % Constant value B 206 CONST P 206 - $b6 % Constant value + 8000 % Constant value B 207 CONST P 207 - $b9 % Constant value + 0.38 % Constant value -B 208 CONST -P 208 - $Cp % Constant value +B 208 CUMC + 301.2 + 19.1 -B 209 CONST -P 209 - $FuelPrice % Constant value +B 209 CUMC + 301.3 + 31.1 + 32.1 -B 210 GT - 268.1 +B 210 CUMC + 301.3 + 25.1 + 26.1 + +B 211 CUMC + 301.2 + 25.1 + 26.1 + +B 212 CUMC + 301.4 + 19.1 + +B 213 CUMC + 301.3 + 305.1 + +B 214 CUMC + 301.3 + 19.1 + +B 215 CUMC + 301.2 + 31.1 + 32.1 + +B 216 CUMC + 301.2 + 305.1 + +B 217 SCREEN + 249.1 +P 217 + '*' % Format + 'Total Electricity Cost in a Year (CAD)' % Headline + +B 218 SCREEN + 246.1 +P 218 + '*' % Format + 'Total CO2 Emissions from Auxiliary Heater (g)' % Headline + +B 219 SCREEN + 247.1 +P 219 + '*' % Format + 'Total Fossil Fuel consumption in a Year (m3)' % Headline + +B 220 SCREEN + 251.1 +P 220 + '*' % Format + 'Total Cost of the Auxiliary Heater Fuel in a Year (CAD)' % Headline + +B 221 SCREEN + 235.1 +P 221 + '*' % Format + 'HP Seasonal COP' % Headline + +B 222 SCREEN + 248.1 +P 222 + '*' % Format + 'Total CO2 Emissions from Electricity Grid (g)' % Headline + +B 223 CHS + 289.1 + +B 224 CHS + 243.1 + +B 225 CHS + 62.1 + +B 226 CHS + 74.1 + +B 227 CHS + 183.1 + +B 228 CHS + 234.1 + +B 229 CHS + 70.1 + +B 230 CHS + 11.1 + +B 231 GE + 287.1 + 165.1 +P 231 + 0 % Error tolerance + +B 232 GE + 287.1 + 147.1 +P 232 + 0 % Error tolerance + +B 233 GE + 245.1 + 171.1 +P 233 + 0 % Error tolerance + +B 234 DIV + 265.1 + 15.1 + +B 235 DIV + 250.1 + 250.2 + +B 236 DIV + 204.1 152.1 -B 211 GT - 217.1 - 180.1 - -B 212 GT - 271.2 - 193.1 - -B 213 GT - 107.1 - 194.1 - -B 214 GT +B 237 DIV 3.1 - 163.1 + 158.1 -B 215 DELAY - 280.1 -P 215 - 10 % Initial value +B 238 DIV + 41.1 + 19.1 -B 216 DELAY - 268.5 -P 216 - 0 % Initial value +B 239 DIV + 306.1 + 36.1 -B 217 DELAY - 107.1 -P 217 - 0 % Initial value +B 240 DIV + 19.1 + 188.1 -B 218 CHS - 268.1 +B 241 DIV + 187.1 + 174.1 -B 219 CHS - 162.1 +B 242 DIV + 19.1 + 203.1 -B 220 CHS - 268.4 - -B 221 CHS - 216.1 - -B 222 CHS +B 243 DIV + 187.1 6.1 -B 223 CUMC - 271.2 - 259.1 +B 244 DIV + 237.1 + 150.1 -B 224 CUMC - 271.4 - 259.1 +B 245 DIV + 239.1 + 281.1 -B 225 CUMC - 271.2 - 120.1 - 140.1 +B 246 CUM + 32.1 -B 226 CUMC - 271.3 - 259.1 +B 247 CUM + 305.1 -B 227 CUMC - 271.2 - 110.1 - 104.1 +B 248 CUM + 31.1 -B 228 CUMC - 271.3 - 255.1 +B 249 CUM + 25.1 -B 229 CUMC - 271.3 - 110.1 - 104.1 +B 250 CUM + 41.1 + 19.1 -B 230 CUMC - 271.2 - 255.1 +B 251 CUM + 26.1 -B 231 CUMC - 271.3 - 120.1 - 140.1 +B 252 HXS + 77.1 + 33.1 + 300.2 + 207.1 +P 252 + 1 % Mode + 5000 % Overall heat transfer coefficient + 1007 % Specific heat of side 1 fluid + 1007 % Specific heat of side 2 fluid -B 232 WRITE - 228.1 - 228.2 -P 232 - 2 % Mode - $fileOut9 % File name +B 253 GT + 69.1 + 234.1 + +B 254 GT + 301.2 + 168.1 + +B 255 GT + 41.1 + 185.1 + +B 256 GT + 278.1 + 181.1 + +B 257 GT + 289.1 + 157.1 + +B 258 GT + 245.1 + 202.1 + +B 259 EXPG + 176.1 + 154.1 + +B 260 EXPG + 79.1 + 156.1 + +B 261 EXPG + 176.1 + 173.1 + +B 262 EXPG + 79.1 + 184.1 + +B 263 EXPG + 252.1 + 191.1 + +B 264 EXPG + 252.1 + 167.1 + +B 265 READ +P 265 + 1 % Number of values to be read per record + 0 % Number of records to be skipped on the first call + $HeatingDemand % File name '*' % Fortran format -B 233 WRITE - 227.1 - 227.2 - 227.3 -P 233 - 2 % Mode - $fileOut3 % File name +B 266 READ +P 266 + 1 % Number of values to be read per record + 0 % Number of records to be skipped on the first call + $WaterTemperature % File name '*' % Fortran format -B 234 WRITE - 226.1 - 226.2 -P 234 - 2 % Mode - $fileOut6 % File name - '*' % Fortran format - -B 235 WRITE - 230.1 - 230.2 -P 235 - 2 % Mode - $fileOut4 % File name - '*' % Fortran format - -B 236 WRITE - 223.1 - 223.2 -P 236 +B 267 WRITE + 208.1 + 208.2 +P 267 2 % Mode $fileOut8 % File name '*' % Fortran format -B 237 WRITE - 231.1 - 231.2 - 231.3 -P 237 - 2 % Mode - $fileOut2 % File name - '*' % Fortran format - -B 238 WRITE - 229.1 - 229.2 - 229.3 -P 238 +B 268 WRITE + 210.1 + 210.2 + 210.3 +P 268 2 % Mode $fileOut7 % File name '*' % Fortran format -B 239 WRITE - 271.1 - 271.2 - 271.3 - 224.1 - 224.2 -P 239 - 2 % Mode - $fileOut10 % File name - '*' % Fortran format - -B 240 WRITE - 225.1 - 225.2 - 225.3 -P 240 - 2 % Mode - $fileOut5 % File name - '*' % Fortran format - -B 241 WRITE - 271.1 - 271.2 - 271.3 - 271.4 - 271.5 - 107.1 - 123.1 - 9.1 - 128.1 - 100.1 - 268.1 - 268.2 - 268.3 - 268.4 - 268.5 - 258.1 - 260.1 +B 269 WRITE + 301.1 + 301.2 + 301.3 + 301.4 + 301.5 + 41.1 + 265.1 11.1 - 263.1 - 110.1 - 104.1 - 120.1 - 140.1 - 92.1 - 277.1 - 269.2 - 96.1 - 90.1 -P 241 + 234.1 + 30.1 + 87.1 + 65.1 + 64.1 + 73.1 + 289.2 + 289.1 + 244.1 + 285.1 + 68.1 + 75.1 + 80.1 + 74.1 + 66.1 + 82.1 + 63.1 + 62.1 + 287.4 + 288.2 + 288.1 + 287.1 + 287.2 + 287.3 + 287.4 + 287.5 + 19.1 + 238.1 + 280.2 + 25.1 + 26.1 + 31.1 + 32.1 +P 269 2 % Mode $fileOut1 % File name '*' % Fortran format -B 242 GE - 3.1 - 160.1 -P 242 - 0 % Error tolerance +B 270 WRITE + 213.1 + 213.2 +P 270 + 2 % Mode + $fileOut9 % File name + '*' % Fortran format -B 243 GE - 268.1 - 198.1 -P 243 - 0 % Error tolerance +B 271 WRITE + 301.1 + 301.2 + 301.3 + 212.1 + 212.2 +P 271 + 2 % Mode + $fileOut10 % File name + '*' % Fortran format -B 244 GE - 268.1 - 202.1 -P 244 - 0 % Error tolerance +B 272 WRITE + 209.1 + 209.2 + 209.3 +P 272 + 2 % Mode + $fileOut2 % File name + '*' % Fortran format -B 245 EXPG - 280.1 - 157.1 +B 273 WRITE + 211.1 + 211.2 + 211.3 +P 273 + 2 % Mode + $fileOut3 % File name + '*' % Fortran format -B 246 EXPG - 280.1 - 153.1 +B 274 WRITE + 214.1 + 214.2 +P 274 + 2 % Mode + $fileOut6 % File name + '*' % Fortran format -B 247 EXPG - 187.1 - 197.1 +B 275 WRITE + 216.1 + 216.2 +P 275 + 2 % Mode + $fileOut4 % File name + '*' % Fortran format -B 248 EXPG - 102.1 - 186.1 +B 276 WRITE + 215.1 + 215.2 + 215.3 +P 276 + 2 % Mode + $fileOut5 % File name + '*' % Fortran format -B 249 EXPG - 102.1 - 173.1 +B 277 DELAY + 252.1 +P 277 + 10 % Initial value -B 250 EXPG - 187.1 - 200.1 +B 278 DELAY + 41.1 +P 278 + 0 % Initial value -B 251 LT - 268.1 - 203.1 +B 279 DELAY + 288.1 +P 279 + 25 % Initial value -B 252 LT - 3.1 - 160.1 - -B 253 LT - 271.2 - 169.1 - -B 254 LE - 3.1 - 163.1 -P 254 - 0 % Error tolerance - -B 255 GAIN - 261.1 -P 255 - 300 % Gain factor g - -B 256 GAIN - 171.1 -P 256 - 1000 % Gain factor g - -B 257 GAIN - 277.1 -P 257 - 1000 % Gain factor g - -B 258 ATT - 268.5 -P 258 - 3600000 % Attenuation factor a - -B 259 ATT - 123.1 -P 259 - 12 % Attenuation factor a - -B 260 ATT - 99.1 -P 260 - 3600000 % Attenuation factor a - -B 261 ATT - 11.1 -P 261 - $FuelDensity % Attenuation factor a - -B 262 ATT - 11.1 -P 262 - $FuelDensity % Attenuation factor a - -B 263 ATT - 115.1 -P 263 - 1000 % Attenuation factor a - -B 264 ATT - 115.1 -P 264 - $AuxHeaterEfficiency % Attenuation factor a - -B 265 ATT - 11.1 -P 265 - $FuelDensity % Attenuation factor a - -B 268 TANKST - 187.1 - 128.1 - 92.1 - 100.1 - 93.1 - 272.1 -P 268 - $TESCapacity % Tank volume - 4 % Number of temperature nodes - $TESDiameter % Tank diameter - $Cp % Specfic heat of fluid - $Rhow % Fluid density - 0 % Overall heat-loss coefficient - 1 % Effective heat conductivity - 30 % Initial tank temperature - -B 269 GENGT2 - 270.1 - 270.3 - 270.4 - 270.5 - 270.7 - 270.8 - 271.1 - 271.2 - 271.3 - 271.4 -P 269 +B 280 GENGT2 + 308.1 + 308.3 + 308.4 + 308.5 + 308.7 + 308.8 + 301.1 + 301.2 + 301.3 + 301.4 +P 280 45.5 % Latitude -73.62 % Longitude -5 % UTC Time zone @@ -984,13 +987,105 @@ P 269 2 % Maximum allowed mean temperature deviation 100 % Maximum number of iterations -B 270 MTM2 - 271.2 -P 270 - 'Montreal' % Location +B 281 INT + 236.1 -B 271 CLOCK -P 271 +B 282 ATT + 182.1 +P 282 + 1000 % Attenuation factor a + +B 283 ATT + 244.1 +P 283 + $FuelDensity % Attenuation factor a + +B 284 ATT + 244.1 +P 284 + $FuelDensity % Attenuation factor a + +B 285 ATT + 3.1 +P 285 + 1000 % Attenuation factor a + +B 286 ATT + 244.1 +P 286 + $FuelDensity % Attenuation factor a + +B 287 TANKST + 61.1 + 62.1 + 66.1 + 74.1 + 78.1 + 307.1 +P 287 + $TESCapacity % Tank volume + 4 % Number of temperature nodes + $TESDiameter % Tank diameter + $Cp % Specfic heat of fluid + $Rhow % Fluid density + 0 % Overall heat-loss coefficient + 1 % Effective heat conductivity + 30 % Initial tank temperature + +B 288 MIXER % Point 11 + 287.4 + 30.1 + 63.1 + 82.1 + +B 289 MIXER % Point 5 + 65.1 + 87.1 + 73.1 + 64.1 + +B 294 LT + 301.2 + 159.1 + +B 295 LT + 11.1 + 234.1 + +B 296 LT + 11.1 + 234.1 + +B 297 LT + 245.1 + 171.1 + +B 298 LT + 289.1 + 180.1 + +B 299 LE + 245.1 + 202.1 +P 299 + 0 % Error tolerance + +B 300 UBHLOSS + 206.1 + 266.1 +P 300 + 400 % bp1 + 0.5 % bp2 + 0.05 % bp3 + 3 % bp4 + 4200 % bp5 + 1000 % bp6 + 5 % bp7 + 20 % bp8 + 300 % bp9 + +B 301 CLOCK +P 301 $StartYear % Start year $StartMonth % Start month $StartDay % Start day @@ -1006,62 +1101,35 @@ P 271 5 % Increment 'm' % Unit -B 272 SOY - 271.1 - 271.2 - 271.3 - 271.4 - 271.5 - 271.6 +B 302 INV + 257.1 -B 273 INV - 244.1 +B 303 INV + 231.1 -B 274 INV - 243.1 +B 304 INV + 232.1 -B 275 INV - 210.1 +B 305 GAIN + 284.1 +P 305 + 300 % Gain factor g -B 276 READ -P 276 - 1 % Number of values to be read per record - 0 % Number of records to be skipped on the first call - $HeatingDemand % File name - '*' % Fortran format - -B 277 READ -P 277 - 1 % Number of values to be read per record - 0 % Number of records to be skipped on the first call - $WaterTemperature % File name - '*' % Fortran format - -B 278 UBHLOSS - 192.1 - 276.1 -P 278 - 300 % L - 0.5 % D - 0.05 % U - 3 % Tg - 4190 % cp - 1000 % rho - 5 % delta_T - 20 % dx - 300 % df - -B 279 INT - 10.1 - -B 280 HXS - 101.1 - 139.1 - 278.2 +B 306 GAIN 204.1 -P 280 - 1 % Mode - 5000 % Overall heat transfer coefficient - 1007 % Specific heat of side 1 fluid - 1007 % Specific heat of side 2 fluid +P 306 + 1000 % Gain factor g + +B 307 SOY + 301.1 + 301.2 + 301.3 + 301.4 + 301.5 + 301.6 + +B 308 MTM2 + 301.2 +P 308 + 'Montreal' % Location diff --git a/data/energy_systems/heat_pumps/w2w_series.txt b/data/energy_systems/heat_pumps/w2w_series.txt index 1f1b6120..1528a0ee 100644 --- a/data/energy_systems/heat_pumps/w2w_series.txt +++ b/data/energy_systems/heat_pumps/w2w_series.txt @@ -1,510 +1,470 @@ +B 65 INT + 198.1 -B 1 DIV - 256.1 +B 66 CUM + 242.1 + +B 67 CUM + 95.1 + +B 68 CUM + 84.1 + +B 69 CUM + 247.1 + +B 70 CUM + 96.1 + +B 71 CUM + 91.1 + +B 72 CUM + 81.1 + 112.1 + +B 73 MUL + 181.1 135.1 -B 2 DIV - 86.1 - 86.2 - -B 3 DIV - 1.1 - 279.1 - -B 4 DIV - 257.1 - 119.1 - -B 5 DIV - 123.1 - 168.1 - -B 6 DIV - 179.1 - 127.1 - -B 7 DIV - 123.1 - 174.1 - -B 8 DIV - 179.1 - 167.1 - -B 9 DIV - 107.1 - 123.1 - -B 10 DIV - 171.1 - 147.1 - -B 11 DIV - 264.1 - 146.1 - -B 76 SCREEN - 2.1 -P 76 - '*' % Format - 'HP Seasonal COP' % Headline - -B 77 SCREEN - 87.1 -P 77 - '*' % Format - 'Total Cost of the Auxiliary Heater Fuel in a Year (CAD)' % Headline - -B 78 SCREEN - 85.1 -P 78 - '*' % Format - 'Total Electricity Cost in a Year (CAD)' % Headline - -B 79 SCREEN - 84.1 -P 79 - '*' % Format - 'Total Fossil Fuel consumption in a Year (m3)' % Headline - -B 80 SCREEN - 89.1 -P 80 - '*' % Format - 'Total CO2 Emissions from Electricity Grid (g)' % Headline - -B 81 SCREEN - 83.1 -P 81 - '*' % Format - 'Total CO2 Emissions from Auxiliary Heater (g)' % Headline - -B 82 SCREEN - 88.1 -P 82 - '*' % Format - 'Total Electricuty Demand of Heat Pumps in a year (kWh)' % Headline - -B 83 CUM - 140.1 - -B 84 CUM - 255.1 - -B 85 CUM - 110.1 - -B 86 CUM - 107.1 - 123.1 - -B 87 CUM - 104.1 - -B 88 CUM - 259.1 - -B 89 CUM - 120.1 - -B 90 SUM - 141.1 - 108.1 - 105.1 - 114.1 - 125.1 - 131.1 - 134.1 - 112.1 - 130.1 - 111.1 - 201.1 - -B 91 SUM - 219.1 - 268.1 - -B 92 SUM - 132.1 - 143.1 - -B 93 SUM - 148.1 - 269.2 - -B 94 SUM - 187.1 - 220.1 - -B 95 SUM - 253.1 - 212.1 - -B 96 SUM - 117.1 - 116.1 - 136.1 - 133.1 - 124.1 - 122.1 - 118.1 - 142.1 - 106.1 - 121.1 - 150.1 - -B 97 SUM - 164.1 - 218.1 - -B 98 SUM - 113.1 - 274.1 - -B 99 SUM - 268.5 - 221.1 - -B 100 SUM - 103.1 - 155.1 - -B 101 SUM - 222.1 - 215.1 - -B 102 SUM - 109.1 - 129.1 - 137.1 - -B 103 MUL - 4.1 - 199.1 - -B 104 MUL - 265.1 - 151.1 - 209.1 - -B 105 MUL - 250.1 - 149.1 - -B 106 MUL - 187.1 - 102.1 - 178.1 - -B 107 MUL - 279.1 - 96.1 - 95.1 - 98.1 - -B 108 MUL - 280.1 - 190.1 - -B 109 MUL - 160.1 - 252.1 - -B 110 MUL - 7.1 - 175.1 - -B 111 MUL - 280.1 - 187.1 - 102.1 - 188.1 - -B 112 MUL - 280.1 - 102.1 - 145.1 - -B 113 MUL - 273.1 - 138.1 - -B 114 MUL - 187.1 - 144.1 - -B 115 MUL - 100.1 - 177.1 - 97.1 - 251.1 - -B 116 MUL - 280.1 - 172.1 - -B 117 MUL +B 74 MUL 191.1 - 245.1 + 179.1 + 194.1 + 226.1 -B 118 MUL - 280.1 - 187.1 - 184.1 +B 75 MUL + 140.1 + 265.1 -B 119 MUL - 208.1 - 154.1 +B 76 MUL + 181.1 + 196.1 + 128.1 -B 120 MUL - 5.1 - 205.1 +B 77 MUL + 196.1 + 169.1 -B 121 MUL - 280.1 - 187.1 - 102.1 - 176.1 - -B 122 MUL - 102.1 - 165.1 - -B 123 MUL - 279.1 - 90.1 - 95.1 - 98.1 - -B 124 MUL - 166.1 - 249.1 - -B 125 MUL - 189.1 - 248.1 - -B 126 MUL - 102.1 - 279.1 - -B 127 MUL - 139.1 - 156.1 - -B 128 MUL - 213.1 - 126.1 - -B 129 MUL - 3.1 - 242.1 - 254.1 - -B 130 MUL - 187.1 - 102.1 - 207.1 - -B 131 MUL - 102.1 - 206.1 - -B 132 MUL - 91.1 - 210.1 - -B 133 MUL - 187.1 - 183.1 - -B 134 MUL - 280.1 - 187.1 +B 78 MUL + 263.1 159.1 -B 135 MUL +B 79 MUL + 246.1 181.1 - 94.1 + 146.1 -B 136 MUL - 247.1 - 196.1 +B 80 MUL + 268.1 + 148.1 -B 137 MUL - 163.1 - 214.1 +B 81 MUL + 65.1 + 187.1 + 189.1 + 192.1 -B 138 MUL - 273.1 - 243.1 - 211.1 +B 82 MUL + 145.1 + 264.1 -B 139 MUL - 8.1 +B 83 MUL + 155.1 195.1 -B 140 MUL - 262.1 - 158.1 - 161.1 +B 84 MUL + 202.1 + 157.1 -B 141 MUL - 182.1 +B 85 MUL + 219.1 + 89.1 + +B 86 MUL + 196.1 + 65.1 + +B 87 MUL 246.1 + 196.1 + 158.1 -B 142 MUL - 280.1 - 102.1 +B 88 MUL + 246.1 + 152.1 + +B 89 MUL + 219.1 + 257.1 + 227.1 + +B 90 MUL + 172.1 + 267.1 + +B 91 MUL + 207.1 + 174.1 + +B 92 MUL + 199.1 + 256.1 + 269.1 + +B 93 MUL + 150.1 + 230.1 + +B 94 MUL + 200.1 + 183.1 + +B 95 MUL + 249.1 + 143.1 + 126.1 + +B 96 MUL + 251.1 + 177.1 + 171.1 + +B 97 MUL + 147.1 + 225.1 + +B 98 MUL + 203.1 + 137.1 + +B 99 MUL + 246.1 + 181.1 + 196.1 + 167.1 + +B 100 MUL + 149.1 + 164.1 + +B 101 MUL + 181.1 + 129.1 + +B 102 MUL 170.1 + 266.1 -B 143 MUL - 275.1 +B 103 MUL + 246.1 + 196.1 + 178.1 + +B 104 MUL + 221.1 + 127.1 + +B 105 MUL + 98.1 + 142.1 + +B 106 MUL + 196.1 + 176.1 + +B 107 MUL + 246.1 + 181.1 + 196.1 + 153.1 + +B 108 MUL + 229.1 + 86.1 + +B 109 MUL + 197.1 + 228.1 + +B 110 MUL + 181.1 + 196.1 + 154.1 + +B 111 MUL + 246.1 + 181.1 + 121.1 + +B 112 MUL + 65.1 185.1 + 189.1 + 192.1 + +B 113 MUL + 246.1 + 119.1 + +B 114 CHS + 278.4 + +B 115 CHS + 261.1 + +B 116 CHS + 144.1 + +B 117 CHS + 278.1 + +B 118 CHS + 206.1 + +B 119 CONST +P 119 + $b2 % Constant value + +B 120 CONST +P 120 + 2 % Constant value + +B 121 CONST +P 121 + $a7 % Constant value + +B 122 CONST +P 122 + 0.40 % Constant value + +B 123 CONST +P 123 + $HPNominalCapacity % Constant value + +B 124 CONST +P 124 + $a11 % Constant value + +B 125 CONST +P 125 + 2 % Constant value + +B 126 CONST +P 126 + $FuelEF % Constant value + +B 127 CONST +P 127 + 25 % Constant value + +B 128 CONST +P 128 + $b9 % Constant value + +B 129 CONST +P 129 + $b4 % Constant value + +B 130 CONST +P 130 + 10 % Constant value + +B 131 CONST +P 131 + $BuildingSuppTemp % Constant value + +B 132 CONST +P 132 + 12 % Constant value + +B 133 CONST +P 133 + $MaximumHPEnergyInput % Constant value + +B 134 CONST +P 134 + 2 % Constant value + +B 135 CONST +P 135 + $a4 % Constant value + +B 136 CONST +P 136 + 0 % Constant value + +B 137 CONST +P 137 + 7.13 % Constant value + +B 138 CONST +P 138 + $FuelLHV % Constant value + +B 139 CONST +P 139 + 2 % Constant value + +B 140 CONST +P 140 + $b1 % Constant value + +B 141 CONST +P 141 + 2 % Constant value + +B 142 CONST +P 142 + $Cp % Constant value + +B 143 CONST +P 143 + 300 % Constant value B 144 CONST P 144 - $b4 % Constant value + $TemperatureDifference % Constant value B 145 CONST P 145 - $b8 % Constant value + $b5 % Constant value B 146 CONST P 146 - $FuelLHV % Constant value + $b7 % Constant value B 147 CONST P 147 - $HPNominalCapacity % Constant value + $LowestPossibleLoadFlow % Constant value B 148 CONST P 148 - 10 % Constant value + $a3 % Constant value B 149 CONST P 149 - $b3 % Constant value + $Cp % Constant value B 150 CONST P 150 - $a11 % Constant value + $HighestPossibleLoadFlow % Constant value B 151 CONST P 151 - 300 % Constant value + $MaximumHPEnergyInput % Constant value B 152 CONST P 152 - 40 % Constant value + $a2 % Constant value B 153 CONST P 153 - 2 % Constant value + $a10 % Constant value B 154 CONST P 154 - $TemperatureDifference % Constant value + $a9 % Constant value B 155 CONST P 155 - 0 % Constant value + $Cp % Constant value B 156 CONST P 156 - $Cp % Constant value + $HPDisactivationTemperature % Constant value B 157 CONST P 157 - 2 % Constant value + $ElecGridEF % Constant value B 158 CONST P 158 - 300 % Constant value + $b8 % Constant value B 159 CONST P 159 - $b7 % Constant value + $b3 % Constant value B 160 CONST P 160 - $LowestPossibleLoadFlow % Constant value + 40 % Constant value B 161 CONST P 161 - $FuelEF % Constant value + 5 % Constant value B 162 CONST P 162 - $TemperatureDifference % Constant value + 2 % Constant value B 163 CONST P 163 - $HighestPossibleLoadFlow % Constant value + $MaximumHPEnergyInput % Constant value B 164 CONST P 164 - $BuildingSuppTemp % Constant value + $TemperatureDifference % Constant value B 165 CONST P 165 - $a6 % Constant value + 40 % Constant value B 166 CONST P 166 - $a5 % Constant value + $HPNominalCapacity % Constant value B 167 CONST P 167 - $HPNominalCapacity % Constant value + $b10 % Constant value B 168 CONST P 168 - 12 % Constant value + 0 % Constant value B 169 CONST P 169 - 5 % Constant value + $a6 % Constant value B 170 CONST P 170 - $a8 % Constant value + $a1 % Constant value B 171 CONST P 171 - $MaximumHPEnergyInput % Constant value + $FuelPrice % Constant value B 172 CONST P 172 - $a2 % Constant value + $a5 % Constant value B 173 CONST P 173 - 2 % Constant value + 12 % Constant value B 174 CONST P 174 - 12 % Constant value + $ElectricityPrice % Constant value B 175 CONST P 175 - $ElectricityPrice % Constant value + $HPReactivationTemperature % Constant value B 176 CONST P 176 - $a10 % Constant value + $b6 % Constant value B 177 CONST P 177 - $Cp % Constant value + 300 % Constant value B 178 CONST P 178 - $a9 % Constant value + $a8 % Constant value B 179 CONST P 179 - $MaximumHPEnergyInput % Constant value + $Cp % Constant value B 180 CONST P 180 @@ -512,467 +472,545 @@ P 180 B 181 CONST P 181 - $Cp % Constant value + $HPSupTemp % Constant value B 182 CONST P 182 - $b1 % Constant value + 9 % Constant value B 183 CONST P 183 - $a4 % Constant value + 1 % Constant value B 184 CONST P 184 - $a7 % Constant value - -B 185 CONST -P 185 - 25 % Constant value - -B 186 CONST -P 186 - 2 % Constant value - -B 187 CONST -P 187 - $HPSupTemp % Constant value - -B 188 CONST -P 188 - $b10 % Constant value - -B 189 CONST -P 189 - $b5 % Constant value - -B 190 CONST -P 190 - $b2 % Constant value - -B 191 CONST -P 191 - $a1 % Constant value - -B 192 CONST -P 192 - $MaximumHPEnergyInput % Constant value - -B 193 CONST -P 193 - 9 % Constant value - -B 194 CONST -P 194 - 0 % Constant value - -B 195 CONST -P 195 - 7.13 % Constant value - -B 196 CONST -P 196 - $a3 % Constant value - -B 197 CONST -P 197 - 2 % Constant value - -B 198 CONST -P 198 - $HPReactivationTemperature % Constant value - -B 199 CONST -P 199 - 1 % Constant value - -B 200 CONST -P 200 - 2 % Constant value - -B 201 CONST -P 201 $b11 % Constant value -B 202 CONST -P 202 - $HPDisactivationTemperature % Constant value +B 185 SUM + 75.1 + 113.1 + 78.1 + 101.1 + 82.1 + 106.1 + 79.1 + 87.1 + 76.1 + 99.1 + 184.1 -B 203 CONST -P 203 - 40 % Constant value +B 186 SUM + 278.5 + 115.1 -B 204 CONST -P 204 - 0.40 % Constant value - -B 205 CONST -P 205 - $ElecGridEF % Constant value - -B 206 CONST -P 206 - $b6 % Constant value - -B 207 CONST -P 207 - $b9 % Constant value - -B 208 CONST -P 208 - $Cp % Constant value - -B 209 CONST -P 209 - $FuelPrice % Constant value - -B 210 GT - 268.1 - 152.1 - -B 211 GT - 217.1 - 180.1 - -B 212 GT - 271.2 - 193.1 - -B 213 GT - 107.1 - 194.1 - -B 214 GT - 3.1 - 163.1 - -B 215 DELAY - 280.1 -P 215 - 10 % Initial value - -B 216 DELAY - 268.5 -P 216 - 0 % Initial value - -B 217 DELAY - 107.1 -P 217 - 0 % Initial value - -B 218 CHS - 268.1 - -B 219 CHS - 162.1 - -B 220 CHS - 268.4 - -B 221 CHS - 216.1 - -B 222 CHS - 6.1 - -B 223 CUMC - 271.2 - 259.1 - -B 224 CUMC - 271.4 - 259.1 - -B 225 CUMC - 271.2 - 120.1 - 140.1 - -B 226 CUMC - 271.3 - 259.1 - -B 227 CUMC - 271.2 +B 187 SUM + 102.1 + 88.1 + 80.1 + 73.1 + 90.1 + 77.1 + 111.1 + 103.1 110.1 - 104.1 + 107.1 + 124.1 -B 228 CUMC - 271.3 - 255.1 +B 188 SUM + 118.1 + 262.1 -B 229 CUMC - 271.3 - 110.1 - 104.1 - -B 230 CUMC - 271.2 - 255.1 - -B 231 CUMC - 271.3 - 120.1 - 140.1 - -B 232 WRITE - 228.1 - 228.2 -P 232 - 2 % Mode - $fileOut9 % File name - '*' % Fortran format - -B 233 WRITE - 227.1 - 227.2 - 227.3 -P 233 - 2 % Mode - $fileOut3 % File name - '*' % Fortran format - -B 234 WRITE - 226.1 - 226.2 -P 234 - 2 % Mode - $fileOut6 % File name - '*' % Fortran format - -B 235 WRITE - 230.1 - 230.2 -P 235 - 2 % Mode - $fileOut4 % File name - '*' % Fortran format - -B 236 WRITE - 223.1 - 223.2 -P 236 - 2 % Mode - $fileOut8 % File name - '*' % Fortran format - -B 237 WRITE - 231.1 - 231.2 - 231.3 -P 237 - 2 % Mode - $fileOut2 % File name - '*' % Fortran format - -B 238 WRITE - 229.1 - 229.2 - 229.3 -P 238 - 2 % Mode - $fileOut7 % File name - '*' % Fortran format - -B 239 WRITE - 271.1 - 271.2 - 271.3 +B 189 SUM 224.1 - 224.2 -P 239 - 2 % Mode - $fileOut10 % File name - '*' % Fortran format + 231.1 -B 240 WRITE - 225.1 - 225.2 - 225.3 -P 240 +B 190 SUM + 130.1 + 277.2 + +B 191 SUM + 94.1 + 168.1 + +B 192 SUM + 85.1 + 220.1 + +B 193 SUM + 109.1 + 104.1 + +B 194 SUM + 131.1 + 117.1 + +B 195 SUM + 181.1 + 114.1 + +B 196 SUM + 97.1 + 92.1 + 93.1 + +B 197 SUM + 116.1 + 278.1 + +B 198 DIV + 151.1 + 123.1 + +B 199 DIV + 201.1 + 65.1 + +B 200 DIV + 243.1 + 100.1 + +B 201 DIV + 241.1 + 83.1 + +B 202 DIV + 112.1 + 173.1 + +B 203 DIV + 133.1 + 166.1 + +B 204 DIV + 81.1 + 112.1 + +B 205 DIV + 253.1 + 138.1 + +B 206 DIV + 133.1 + 105.1 + +B 207 DIV + 112.1 + 132.1 + +B 208 DIV + 72.1 + 72.2 + +B 209 WRITE + 236.1 + 236.2 + 236.3 +P 209 2 % Mode $fileOut5 % File name '*' % Fortran format -B 241 WRITE - 271.1 - 271.2 - 271.3 - 271.4 - 271.5 - 107.1 - 123.1 - 9.1 - 128.1 - 100.1 - 268.1 - 268.2 - 268.3 - 268.4 - 268.5 - 258.1 - 260.1 - 11.1 - 263.1 - 110.1 - 104.1 - 120.1 - 140.1 - 92.1 - 277.1 - 269.2 +B 210 WRITE + 280.1 + 280.2 + 280.3 + 280.4 + 280.5 + 81.1 + 112.1 + 204.1 + 108.1 + 191.1 + 278.1 + 278.2 + 278.3 + 278.4 + 278.5 + 248.1 + 254.1 + 205.1 + 252.1 + 91.1 96.1 - 90.1 -P 241 + 84.1 + 95.1 + 193.1 + 222.1 + 277.2 + 187.1 + 185.1 +P 210 2 % Mode $fileOut1 % File name '*' % Fortran format -B 242 GE - 3.1 +B 211 WRITE + 233.1 + 233.2 +P 211 + 2 % Mode + $fileOut6 % File name + '*' % Fortran format + +B 212 WRITE + 280.1 + 280.2 + 280.3 + 239.1 + 239.2 +P 212 + 2 % Mode + $fileOut10 % File name + '*' % Fortran format + +B 213 WRITE + 235.1 + 235.2 +P 213 + 2 % Mode + $fileOut4 % File name + '*' % Fortran format + +B 214 WRITE + 232.1 + 232.2 +P 214 + 2 % Mode + $fileOut8 % File name + '*' % Fortran format + +B 215 WRITE + 237.1 + 237.2 + 237.3 +P 215 + 2 % Mode + $fileOut7 % File name + '*' % Fortran format + +B 216 WRITE + 234.1 + 234.2 + 234.3 +P 216 + 2 % Mode + $fileOut2 % File name + '*' % Fortran format + +B 217 WRITE + 240.1 + 240.2 + 240.3 +P 217 + 2 % Mode + $fileOut3 % File name + '*' % Fortran format + +B 218 WRITE + 238.1 + 238.2 +P 218 + 2 % Mode + $fileOut9 % File name + '*' % Fortran format + +B 219 INV + 258.1 + +B 220 INV + 257.1 + +B 221 INV + 228.1 + +B 222 READ +P 222 + 1 % Number of values to be read per record + 0 % Number of records to be skipped on the first call + $HeatingDemand % File name + '*' % Fortran format + +B 223 READ +P 223 + 1 % Number of values to be read per record + 0 % Number of records to be skipped on the first call + $WaterTemperature % File name + '*' % Fortran format + +B 224 LT + 280.2 + 161.1 + +B 225 LT + 199.1 + 147.1 + +B 226 LT + 278.1 + 165.1 + +B 227 GT + 260.1 + 180.1 + +B 228 GT + 278.1 160.1 + +B 229 GT + 81.1 + 136.1 + +B 230 GT + 199.1 + 150.1 + +B 231 GT + 280.2 + 182.1 + +B 232 CUMC + 280.2 + 247.1 + +B 233 CUMC + 280.3 + 247.1 + +B 234 CUMC + 280.3 + 84.1 + 95.1 + +B 235 CUMC + 280.2 + 242.1 + +B 236 CUMC + 280.2 + 84.1 + 95.1 + +B 237 CUMC + 280.3 + 91.1 + 96.1 + +B 238 CUMC + 280.3 + 242.1 + +B 239 CUMC + 280.4 + 247.1 + +B 240 CUMC + 280.2 + 91.1 + 96.1 + +B 241 GAIN + 151.1 +P 241 + 1000 % Gain factor g + +B 242 GAIN + 250.1 P 242 - 0 % Error tolerance - -B 243 GE - 268.1 - 198.1 -P 243 - 0 % Error tolerance - -B 244 GE - 268.1 - 202.1 -P 244 - 0 % Error tolerance - -B 245 EXPG - 280.1 - 157.1 - -B 246 EXPG - 280.1 - 153.1 - -B 247 EXPG - 187.1 - 197.1 - -B 248 EXPG - 102.1 - 186.1 - -B 249 EXPG - 102.1 - 173.1 - -B 250 EXPG - 187.1 - 200.1 - -B 251 LT - 268.1 - 203.1 - -B 252 LT - 3.1 - 160.1 - -B 253 LT - 271.2 - 169.1 - -B 254 LE - 3.1 - 163.1 -P 254 - 0 % Error tolerance - -B 255 GAIN - 261.1 -P 255 300 % Gain factor g -B 256 GAIN - 171.1 -P 256 +B 243 GAIN + 222.1 +P 243 1000 % Gain factor g -B 257 GAIN - 277.1 -P 257 - 1000 % Gain factor g +B 246 HXS + 188.1 + 98.1 + 259.2 + 122.1 +P 246 + 1 % Mode + 5000 % Overall heat transfer coefficient + 1007 % Specific heat of side 1 fluid + 1007 % Specific heat of side 2 fluid -B 258 ATT - 268.5 -P 258 - 3600000 % Attenuation factor a - -B 259 ATT - 123.1 -P 259 +B 247 ATT + 112.1 +P 247 12 % Attenuation factor a -B 260 ATT - 99.1 -P 260 +B 248 ATT + 278.5 +P 248 3600000 % Attenuation factor a -B 261 ATT - 11.1 -P 261 +B 249 ATT + 205.1 +P 249 $FuelDensity % Attenuation factor a -B 262 ATT - 11.1 -P 262 +B 250 ATT + 205.1 +P 250 $FuelDensity % Attenuation factor a -B 263 ATT - 115.1 -P 263 +B 251 ATT + 205.1 +P 251 + $FuelDensity % Attenuation factor a + +B 252 ATT + 74.1 +P 252 1000 % Attenuation factor a -B 264 ATT - 115.1 -P 264 +B 253 ATT + 74.1 +P 253 $AuxHeaterEfficiency % Attenuation factor a -B 265 ATT - 11.1 -P 265 - $FuelDensity % Attenuation factor a +B 254 ATT + 186.1 +P 254 + 3600000 % Attenuation factor a -B 268 TANKST - 187.1 - 128.1 - 92.1 - 100.1 - 93.1 - 272.1 -P 268 - $TESCapacity % Tank volume - 4 % Number of temperature nodes - $TESDiameter % Tank diameter - $Cp % Specfic heat of fluid - $Rhow % Fluid density - 0 % Overall heat-loss coefficient - 1 % Effective heat conductivity - 30 % Initial tank temperature +B 255 MTM2 + 280.2 +P 255 + 'Montreal' % Location -B 269 GENGT2 - 270.1 - 270.3 - 270.4 - 270.5 - 270.7 - 270.8 - 271.1 - 271.2 - 271.3 - 271.4 +B 256 GE + 199.1 + 147.1 +P 256 + 0 % Error tolerance + +B 257 GE + 278.1 + 175.1 +P 257 + 0 % Error tolerance + +B 258 GE + 278.1 + 156.1 +P 258 + 0 % Error tolerance + +B 259 UBHLOSS + 163.1 + 223.1 +P 259 + 300 % bp1 + 0.5 % bp2 + 0.05 % bp3 + 3 % bp4 + 4190 % bp5 + 1000 % bp6 + 5 % bp7 + 20 % bp8 + 300 % bp9 + +B 260 DELAY + 81.1 +P 260 + 0 % Initial value + +B 261 DELAY + 278.5 +P 261 + 0 % Initial value + +B 262 DELAY + 246.1 +P 262 + 10 % Initial value + +B 263 EXPG + 181.1 + 120.1 + +B 264 EXPG + 196.1 + 141.1 + +B 265 EXPG + 246.1 + 125.1 + +B 266 EXPG + 246.1 + 139.1 + +B 267 EXPG + 196.1 + 162.1 + +B 268 EXPG + 181.1 + 134.1 + +B 269 LE + 199.1 + 150.1 P 269 + 0 % Error tolerance + +B 270 SCREEN + 71.1 +P 270 + '*' % Format + 'Total Electricity Cost in a Year (CAD)' % Headline + +B 271 SCREEN + 67.1 +P 271 + '*' % Format + 'Total CO2 Emissions from Auxiliary Heater (g)' % Headline + +B 272 SCREEN + 208.1 +P 272 + '*' % Format + 'HP Seasonal COP' % Headline + +B 273 SCREEN + 66.1 +P 273 + '*' % Format + 'Total Fossil Fuel consumption in a Year (m3)' % Headline + +B 274 SCREEN + 69.1 +P 274 + '*' % Format + 'Total Electricuty Demand of Heat Pumps in a year (kWh)' % Headline + +B 275 SCREEN + 68.1 +P 275 + '*' % Format + 'Total CO2 Emissions from Electricity Grid (g)' % Headline + +B 276 SCREEN + 70.1 +P 276 + '*' % Format + 'Total Cost of the Auxiliary Heater Fuel in a Year (CAD)' % Headline + +B 277 GENGT2 + 255.1 + 255.3 + 255.4 + 255.5 + 255.7 + 255.8 + 280.1 + 280.2 + 280.3 + 280.4 +P 277 45.5 % Latitude -73.62 % Longitude -5 % UTC Time zone @@ -984,13 +1022,33 @@ P 269 2 % Maximum allowed mean temperature deviation 100 % Maximum number of iterations -B 270 MTM2 - 271.2 -P 270 - 'Montreal' % Location +B 278 TANKST + 181.1 + 108.1 + 193.1 + 191.1 + 190.1 + 279.1 +P 278 + $TESCapacity % Tank volume + 4 % Number of temperature nodes + $TESDiameter % Tank diameter + $Cp % Specfic heat of fluid + $Rhow % Fluid density + 0 % Overall heat-loss coefficient + 1 % Effective heat conductivity + 30 % Initial tank temperature -B 271 CLOCK -P 271 +B 279 SOY + 280.1 + 280.2 + 280.3 + 280.4 + 280.5 + 280.6 + +B 280 CLOCK +P 280 $StartYear % Start year $StartMonth % Start month $StartDay % Start day @@ -1006,62 +1064,3 @@ P 271 5 % Increment 'm' % Unit -B 272 SOY - 271.1 - 271.2 - 271.3 - 271.4 - 271.5 - 271.6 - -B 273 INV - 244.1 - -B 274 INV - 243.1 - -B 275 INV - 210.1 - -B 276 READ -P 276 - 1 % Number of values to be read per record - 0 % Number of records to be skipped on the first call - $HeatingDemand % File name - '*' % Fortran format - -B 277 READ -P 277 - 1 % Number of values to be read per record - 0 % Number of records to be skipped on the first call - $WaterTemperature % File name - '*' % Fortran format - -B 278 UBHLOSS - 192.1 - 276.1 -P 278 - 300 % L - 0.5 % D - 0.05 % U - 3 % Tg - 4190 % cp - 1000 % rho - 5 % delta_T - 20 % dx - 300 % df - -B 279 INT - 10.1 - -B 280 HXS - 101.1 - 139.1 - 278.2 - 204.1 -P 280 - 1 % Mode - 5000 % Overall heat transfer coefficient - 1007 % Specific heat of side 1 fluid - 1007 % Specific heat of side 2 fluid - diff --git a/exports/energy_systems/heat_pump_export.py b/exports/energy_systems/heat_pump_export.py index bdcfcdae..b6731579 100644 --- a/exports/energy_systems/heat_pump_export.py +++ b/exports/energy_systems/heat_pump_export.py @@ -71,17 +71,6 @@ class HeatPumpExport: Write headers to the various csv file generated by insel :return: """ - file_one_headers = ['Year', ' Month', ' Day', 'Hour', 'Minute', 'HP Heat Output (kW)', - 'HP Electricity Consumption (kW)', 'HP COP', 'TES Charging Rate (kg/s)', - 'TES Discharging Rate (kg/s)', 'TES Node 1 Temperature', 'TES Node 2 Temperature', - 'TES Node 3 Temperature', 'TES Node 4 Temperature', 'TES Energy Content (J)', - 'TES Energy Content (kWh)', 'TES Energy Content Variation (kWh)', - 'Auxiliary Heater Fuel Flow Rate (kg/s)', 'Auxiliary Heater Energy Input (kW)', - 'HP Operational Cost (CAD)', 'Auxiliary Heater Operational Cost (CAD)', - 'Operational CO2 Emissions of HP (g)', - 'Operational CO2 Emissions of Auxiliary Heater (g)', - 'Return Temperature', 'Demand (kW)', 'Test Column'] - header_data = { self._input_data['fileOut1']: ['Year', ' Month', ' Day', 'Hour', 'Minute', 'HP Heat Output (kW)', 'HP Electricity Consumption (kW)', 'HP COP', 'TES Charging Rate (kg/s)', @@ -200,6 +189,7 @@ class HeatPumpExport: and water to water source heat pump :return: """ + self._input_data["a1"] = a_coeff[0] self._input_data["a2"] = a_coeff[1] self._input_data["a3"] = a_coeff[2] diff --git a/exports/energy_systems/water_to_water_hp_export.py b/exports/energy_systems/water_to_water_hp_export.py index 7f7c4750..07787218 100644 --- a/exports/energy_systems/water_to_water_hp_export.py +++ b/exports/energy_systems/water_to_water_hp_export.py @@ -26,22 +26,19 @@ class WaterToWaterHPExport(HeatPumpExport): water_temp = (base_path / 'heat_pumps/wt_hourly3.txt') super().__init__(base_path, city, output_path, template_path, water_temp) - def _extract_model_coff(self, hp_model: str, data_type='heat') -> Union[Tuple[List, List], None]: + def _extract_model_coff(self, hp_model: str) -> Union[List, None]: """ Extracts heat pump coefficient data for a specific - model. e.g 012, 140 + model. e.g ClimateMaster 156 kW, etc :param hp_model: the model type - :param data_type: indicates whether we're extracting cooling - or heating perfarmcn coefficients :return: """ for energy_system in self._city.energy_systems: if energy_system.water_to_water_hp.model == hp_model: - return energy_system.water_to_water_hp.power_demand_coff, \ - energy_system.water_to_water_hp.heat_output_coff + return energy_system.water_to_water_hp.power_demand_coff return None - def execute_insel(self, user_input, hp_model, data_type): + def execute_insel(self, user_input, hp_model): """ Runs insel and produces output files Runs insel and write the necessary files @@ -49,10 +46,7 @@ class WaterToWaterHPExport(HeatPumpExport): values necessary to run insel :param hp_model: a string that indicates the heat pump model to be used e.g. 012, 015 - :param data_type: a string that indicates whether - insel should run for heat or cooling performance - :return: :return: """ - pow_demand_coeff, heat_output_coeff = self._extract_model_coff(hp_model, data_type) - super(WaterToWaterHPExport, self)._run_insel(user_input, pow_demand_coeff, heat_output_coeff, 'w2w.insel') + pow_demand_coeff = self._extract_model_coff(hp_model) + super(WaterToWaterHPExport, self)._run_insel(user_input, pow_demand_coeff, 'w2w.insel') diff --git a/exports/energy_systems_factory.py b/exports/energy_systems_factory.py index e24e0848..f3b52e11 100644 --- a/exports/energy_systems_factory.py +++ b/exports/energy_systems_factory.py @@ -47,7 +47,7 @@ class EnergySystemsExportFactory: .execute_insel(self._user_input, self._hp_model, self._data_type) elif source == 'water': WaterToWaterHPExport(self._base_path, self._city, self._output_path, self._sim_type)\ - .execute_insel(self._user_input, self._hp_model, self._data_type) + .execute_insel(self._user_input, self._hp_model) def export(self, source='air'): """ diff --git a/imports/energy_systems/water_to_water_hp_parameters.py b/imports/energy_systems/water_to_water_hp_parameters.py index 9348832a..346c5e3f 100644 --- a/imports/energy_systems/water_to_water_hp_parameters.py +++ b/imports/energy_systems/water_to_water_hp_parameters.py @@ -125,24 +125,20 @@ class WaterToWaterHPParameters: heat_pump.flow_rate = data['fr'] heat_pump.entering_water_temp = data['ewt'] heat_pump.leaving_water_temp = data['lwt'] - heat_pump.power_demand_coff = self._compute_coefficients(data, data_type='power') - heat_pump.heat_output_coff = self._compute_coefficients(data) - + heat_pump.power_demand_coff = self._compute_coefficients(data) energy_system = EnergySystem(heat_pump.model, 0, [], None) energy_system.water_to_water_hp = heat_pump self._city.add_city_object(energy_system) return self._city - def _compute_coefficients(self, heat_pump_data: Dict, data_type="heat_output") -> List[float]: + def _compute_coefficients(self, heat_pump_data: Dict) -> List[float]: """ Compute heat output and electrical demand coefficients from heating performance data :param heat_pump_data: a dictionary of heat pump data. - :param data_type: string to indicate whether coefficient are power demands - or heat output coefficients. Default is heat output :return: Tuple[Dict, Dict] """ - demand = heat_pump_data['tc'] if data_type == "heat_output" else heat_pump_data['pd'] + demand = [i / j for i, j in zip(heat_pump_data['tc'], heat_pump_data['pd'])] # Compute heat output coefficients popt, _ = curve_fit(self._objective_function, [heat_pump_data['ewt'], heat_pump_data['lwt'], heat_pump_data['fr']],