Feature/add stratified tes

docs/parameters.rst

@@ -254,6 +254,60 @@ Some parameters can be calculated automatically by *pvcompare* and do not need t
 
         b. **input power** and **output power**: kW
 
+* storage_02.csv:
+    1. **age_installed**: year, 0 (for all components of the stratified thermal storage)
+    2. **development_costs**: currency, 0 (for all components of the stratified thermal storage)
+    3. **specific_costs**: currency/unit
+        a. **storage capacity**: 410, See `Danish energy agency's technology data of small-scale hot water tanks [dea_swt] <https://ens.dk/sites/ens.dk/files/Analyser/technology_data_catalogue_for_energy_storage.pdf>`_ on p.66 - However investment costs of stratified TES could be higher.
+        b. **input power** and **output power**: 0
+    4. **c_rate**: factor of total capacity (kWh)
+        a. **storage capacity**: NA (does not apply)
+        b. **input power** and **output power**: 1 (this just means that the whole capacity of the stratified thermal storage would be used during charging and discharging cycles)
+    5. **efficiency**: factor
+        a. **storage capacity**: 1, or "NA" if calculated
+        b. **input power** and **output power**: 1
+    6. **installedCap**: unit 0, or "NA" if calculated
+        a. **storage capacity**: 0, or "NA" if calculated
+        b. **input power** and **output power**: 0
+    7. **lifetime**: year, 30 (applies for all the parameters of the stratified thermal energy storage)
+    8. **specific_costs_om**: currency/unit/year
+        a. **storage capacity**: 16.67, [dea_swt] p.66 - however fix om costs of stratified TES could differ
+        b. **input power** and **output power**: 0
+    9. **dispatch_price**: currency/kWh
+        a. **storage capacity**: NA (does not apply)
+        b. **input power** and **output power**: 0
+    10. **soc_initial**: None or factor
+
+        a. **storage capacity**: None
+
+        b. **input power** and **output power**: NA
+    11. **soc_max**: factor
+
+        a. **storage capacity**: 0.925 (7.5% unused volume see `European Commission study large-scale heating and cooling in EU [EUC_heat] <https://op.europa.eu/en/publication-detail/-/publication/312f0f62-dfbd-11e7-9749-01aa75ed71a1/language-en>`_ p.168 - This applies for large scale TES but could be validated for a small scale storage too.)
+
+        b. **input power** and **output power**: NA
+    12. **soc_min**: factor
+
+        a. **storage capacity**: 0.075 (7.5% unused volume see [EUC_heat] p.168 - This applies for large scale TES but could be validated for a small scale storage too.)
+
+        b. **input power** and **output power**: NA
+    13. **unit**: str
+
+        a. **storage capacity**: kWh
+
+        b. **input power** and **output power**: kW
+    14. **fixed_thermal_losses_relative**: factor
+
+        a. **storage capacity**: "{'file_name': 'None', 'header': 'no_unit', 'unit': ''}", is calculated in pvcompare
+
+        b. **input power** and **output power**: NA (does not apply)
+    15. **fixed_thermal_losses_absolute**: kWh
+
+        a. **storage capacity**: "{'file_name': 'None', 'header': 'no_unit', 'unit': ''}", is calculated in pvcompare
+
+        b. **input power** and **output power**: NA (does not apply)
+
+
 
 2. pvcompare-specific parameters
 ================================

examples/example_user_inputs/mvs_inputs_sector_coupling/csv_elements/energyStorage.csv

@@ -1,8 +1,8 @@
-,unit
-label,str
-optimizeCap,bool
-outflow_direction,str
-inflow_direction,str
-storage_filename,str
-energyVector,str
-type_oemof,str
+,unit,storage_02
+label,str,TES
+optimizeCap,bool,True
+outflow_direction,str,Heat bus
+inflow_direction,str,Heat bus
+storage_filename,str,storage_02.csv
+energyVector,str,Heat
+type_oemof,str,storage

examples/example_user_inputs/mvs_inputs_sector_coupling/csv_elements/storage_02.csv

@@ -0,0 +1,16 @@
+,unit,storage capacity,input power,output power
+unit,str,kWh,kW,kW
+installedCap,kWh,NaN,0,0
+age_installed,year,0,0,0
+lifetime,year,30,30,30
+development_costs,currency,0.0,0.0,0.0
+specific_costs,currency/kWh,410,0,0
+specific_costs_om,currency/kWh/year,16.67,0,0
+dispatch_price,currency/kWh,NaN,0,0
+c_rate,factor of total capacity (kWh),NaN,1,1
+soc_initial,None or factor,None,NaN,NaN
+soc_max,factor,0.925,NaN,NaN
+soc_min,factor,0.075,NaN,NaN
+efficiency,factor,NaN,1,1
+fixed_thermal_losses_relative,factor,"{'file_name': 'None', 'header': 'no_unit', 'unit': ''}",,
+fixed_thermal_losses_absolute,kWh,"{'file_name': 'None', 'header': 'kWh', 'unit': ''}",,

examples/example_user_inputs/pvcompare_inputs/stratified_thermal_storage.csv

@@ -0,0 +1,9 @@
+var_name,var_value,var_unit
+height,,m
+diameter,0.99,m
+temp_h,60,degC
+temp_c,35,degC
+s_iso,100,mm
+lamb_iso,0.03,W/(m*K)
+alpha_inside,10,W/(m2*K)
+alpha_outside,2.31,W/(m2*K)

pvcompare/check_inputs.py

@@ -611,3 +611,58 @@ def load_parameter_from_mvs_file(
         f"mvs input file {mvs_filename}."
     )
     return pvcompare_parameter
+
+
+def add_parameters_to_storage_xx_file(
+    nominal_storage_capacity, loss_rate, storage_csv, user_inputs_mvs_directory=None
+):
+
+    """
+    Enters new parameters into storage_xx.csv
+
+    Parameters
+    ---------
+    nominal_storage_capacity : numeric
+        Maximum amount of stored thermal energy [MWh]
+
+    loss_rate : numeric (sequence or scalar)
+        The relative loss of the storage capacity between two consecutive
+        timesteps [-]
+
+    storage_csv: str
+        Name of the storage specific file
+
+    mvs_input_directory : str
+        directory to "mvs_inputs/"
+
+    Returns
+    -------
+    None
+    """
+
+    if user_inputs_mvs_directory == None:
+        user_inputs_mvs_directory = constants.DEFAULT_USER_INPUTS_MVS_DIRECTORY
+
+    # Read storage_xx.csv from input value
+    storage_xx_path = os.path.join(
+        user_inputs_mvs_directory, "csv_elements", storage_csv
+    )
+    storage_xx = pd.read_csv(storage_xx_path, header=0, index_col=0,)
+
+    parameters = {"installedCap": nominal_storage_capacity, "efficiency": 1 - loss_rate}
+
+    for name, param in parameters.items():
+        # Check if efficiency and nominal storage capacity already exist and if not
+        # replace with calculated value from loss_rate
+        try:
+            int(float(storage_xx.at[name, "storage capacity"]))
+            logging.info(
+                f"The {name} of the storage already exists in {storage_csv}. If you want to calculate the storage's {name} with pvcompare please delete the value in {storage_csv}."
+            )
+        except ValueError:
+            # insert parameter values
+            storage_xx.loc[[name], ["storage capacity"]] = param
+            logging.info(f"The {name} of the storage has been added to {storage_csv}.")
+
+    # Save values in storage_xx.csv
+    storage_xx.to_csv(storage_xx_path)

pvcompare/data/user_inputs_collection/mvs_inputs/csv_elements/storage_02.csv

@@ -0,0 +1,16 @@
+,unit,storage capacity,input power,output power
+unit,str,kWh,kW,kW
+installedCap,kWh,NaN,0,0
+age_installed,year,0,0,0
+lifetime,year,30,30,30
+development_costs,currency,0.0,0.0,0.0
+specific_costs,currency/kWh,410,0,0
+specific_costs_om,currency/kWh/year,16.67,0,0
+dispatch_price,currency/kWh,NaN,0,0
+c_rate,factor of total capacity (kWh),NaN,1,1
+efficiency,factor,1,1,1
+soc_initial,None or factor,None,NaN,NaN
+soc_max,factor,0.925,NaN,NaN
+soc_min,factor,0.075,NaN,NaN
+fixed_thermal_losses_relative,factor,"{'file_name': 'None', 'header': 'no_unit', 'unit': ''}",,
+fixed_thermal_losses_absolute,kWh,"{'file_name': 'None', 'header': 'no_unit', 'unit': ''}",,

pvcompare/data/user_inputs_collection/pvcompare_inputs/stratified_thermal_storage.csv

@@ -0,0 +1,9 @@
+var_name,var_value,var_unit
+height,,m
+diameter,0.99,m
+temp_h,95,degC
+temp_c,60,degC
+s_iso,100,mm
+lamb_iso,0.03,W/(m*K)
+alpha_inside,10,W/(m2*K)
+alpha_outside,2.31,W/(m2*K)

pvcompare/main.py

@@ -12,6 +12,7 @@
 from pvcompare import pv_feedin
 from pvcompare import constants
 from pvcompare import heat_pump_and_chiller
+from pvcompare import stratified_thermal_storage
 from pvcompare import check_inputs
 
 
@@ -179,6 +180,14 @@ def apply_pvcompare(
         weather=weather,
     )
 
+    stratified_thermal_storage.add_strat_tes(
+        weather=weather,
+        lat=latitude,
+        lon=longitude,
+        user_inputs_pvcompare_directory=user_inputs_pvcompare_directory,
+        user_inputs_mvs_directory=user_inputs_mvs_directory,
+    )
+
 
 def apply_mvs(
     scenario_name,