add exergoeconomic analysis

bash.exe.stackdump

@@ -0,0 +1,16 @@
+Stack trace:
+Frame        Function    Args
+000FFFFCD30  00210062B0E (0021028A770, 00210275E51, 00000000059, 000FFFFB710)
+000FFFFCD30  0021004846A (000FFFFCD30, 00200000000, 00000000000, 00000000001)
+000FFFFCD30  002100484A2 (00000000000, 00000000000, 00000000059, 99B853482A57)
+000FFFFCD30  0021006E496 (00210045323, 00210358A70, 00000000000, 0000000000D)
+000FFFFCD30  0021006E4A9 (00210045170, 0021023D7E0, 002100448F2, 000FFFFC910)
+000FFFFCD30  00210070DE4 (00000000013, 00000000001, 000FFFFC910, 00210278640)
+000FFFFCD30  0021005AB65 (000FFFFCA60, 00000000000, 00000000000, 008FFFFFFFF)
+000FFFFCD30  0021005B335 (00000000002, 00000000000, 000FFFFCD30, 00000030201)
+000FFFFCD30  0021005B847 (002100DF73E, 00000000000, 00000000000, 00000000000)
+000FFFFCD30  0021005BB86 (00000000000, 000FFFFCD30, FFFFFFFFFFFFFFC6, 00000000000)
+000FFFFCD30  00210048C0C (00000000000, 00000000000, 00000000000, 00000000000)
+000FFFFFFF0  00210047716 (00000000000, 00000000000, 00000000000, 00000000000)
+000FFFFFFF0  002100477C4 (00000000000, 00000000000, 00000000000, 00000000000)
+End of stack trace

src/tespy/components/basics/cycle_closer.py

@@ -174,3 +174,9 @@ def calc_parameters(self):
         d2 = self.outl[0].fluid.val
         diff = [d1[key] - d2[key] for key in d1.keys()]
         self.fluid_deviation.val = np.linalg.norm(diff)
+
+
+    """+F+F+F+F++++START++++F+F+F+F+"""
+
+    def exergoeconomic_balance(self, T0):
+        return

src/tespy/components/basics/source.py

@@ -122,3 +122,31 @@ def exergy_balance(self, T0):
         }
         self.E_D = np.nan
         self.epsilon = self._calc_epsilon()
+
+
+    """+F+F+F+F++++START++++F+F+F+F+    von Jubran"""
+    def set_source_costs(self, c_tot=None):
+        # c_tot must be in SI unit ($/J)
+
+        if c_tot is not None:
+            self.outl[0].c_tot = c_tot
+            # calculate outlet
+            self.outl[0].Ex_tot = self.outl[0].Ex_physical + self.outl[0].Ex_chemical
+            self.outl[0].C_tot = self.outl[0].c_tot * self.outl[0].Ex_tot
+
+            # approximate costs per exergy unit fot T, M, PH and CH
+            # c_tot = c_therm = c_mech = c_physical = c_chemical
+            self.outl[0].c_therm = c_tot
+            self.outl[0].c_mech = c_tot
+            self.outl[0].c_physical = c_tot
+            self.outl[0].c_chemical = c_tot
+
+            self.outl[0].C_therm = self.outl[0].c_therm * self.outl[0].Ex_therm
+            self.outl[0].C_mech = self.outl[0].c_mech * self.outl[0].Ex_mech
+            self.outl[0].C_physical = self.outl[0].c_physical * self.outl[0].Ex_physical
+            self.outl[0].C_chemical = self.outl[0].c_chemical * self.outl[0].Ex_chemical
+
+    def set_source_costs_standard(self):
+        # determine source costs depending on material, temperature, ...
+        return
+    """+F+F+F+F++++END++++F+F+F+F+    von Jubran"""

src/tespy/components/combustion/base.py

@@ -1399,4 +1399,46 @@ def exergy_balance(self, T0):
 
         self.E_D = self.E_F - self.E_P
         self.epsilon = self._calc_epsilon()
-        self.E_bus = np.nan
+        self.E_bus = {"chemical": np.nan, "physical": np.nan, "massless": np.nan}
+
+
+    """+F+F+F+F++++START++++F+F+F+F+"""
+
+    def exergoeconomic_balance(self, T0):
+        self.C_P = self.outl[0].C_therm - (
+                self.inl[0].C_therm + self.inl[1].C_therm
+        )
+        self.C_F = (
+                self.inl[0].C_chemical + self.inl[1].C_chemical -
+                self.outl[0].C_chemical + self.inl[0].C_mech +
+                self.inl[1].C_mech - self.outl[0].C_mech
+        )
+        self.c_F = self.C_F / self.E_F
+        self.c_P = self.C_P / self.E_P
+        self.C_D = self.c_F * self.E_D
+        self.r = (self.c_P - self.c_F) / self.c_F
+        self.f = self.Z_costs / (self.Z_costs + self.C_D)
+
+    def aux_eqs(self, exergy_cost_matrix, exergy_cost_vector, counter, T0):
+        if self.outl[0].Ex_mech != 0 and self.inl[0].Ex_mech  != 0 and self.inl[1].Ex_mech  != 0:
+            exergy_cost_matrix[counter+0, self.outl[0].Ex_C_col["mech"]] = -1 / self.outl[0].Ex_mech
+            exergy_cost_matrix[counter+0, self.inl[0].Ex_C_col["mech"]] = 1 / self.inl[0].Ex_mech * self.inl[0].m.val / (self.inl[0].m.val + self.inl[1].m.val)
+            exergy_cost_matrix[counter+0, self.inl[1].Ex_C_col["mech"]] = 1 / self.inl[1].Ex_mech * self.inl[1].m.val / (self.inl[0].m.val + self.inl[1].m.val)
+        else:
+            exergy_cost_matrix[counter+0, self.outl[0].Ex_C_col["mech"]] = 1
+
+        if self.outl[0].Ex_chemical != 0 and self.inl[0].Ex_chemical  != 0 and self.inl[1].Ex_chemical  != 0:
+            exergy_cost_matrix[counter+1, self.outl[0].Ex_C_col["chemical"]] = -1 / self.outl[0].Ex_chemical
+            exergy_cost_matrix[counter+1, self.inl[0].Ex_C_col["chemical"]] = 1 / self.inl[0].Ex_chemical * self.inl[0].m.val / (self.inl[0].m.val + self.inl[1].m.val)
+            exergy_cost_matrix[counter+1, self.inl[1].Ex_C_col["chemical"]] = 1 / self.inl[1].Ex_chemical * self.inl[1].m.val / (self.inl[0].m.val + self.inl[1].m.val)
+        elif self.in1[0].Ex_chemical == 0:
+            exergy_cost_matrix[counter+1, self.inl[0].Ex_C_col["chemical"]] = 1
+        elif self.inl[1].Ex_chemical == 0:
+            exergy_cost_matrix[counter+1, self.inl[1].Ex_C_col["chemical"]] = 1
+
+        for i in range(2):
+            exergy_cost_vector[counter+i]=0
+
+        return [exergy_cost_matrix, exergy_cost_vector, counter+2]
+
+    """+F+F+F+F++++END++++F+F+F+F+"""