Construction of a "Drying" custom component

[MMArro]

Hello,

I have been attempting to construct a custom component to represent an industrial dryer. The initial idea is to make it as simple as possible. For this purpose, I have attempted to repurpose the heat exchanger component already existing in Tespy, intending to add a component of mass transfer which would correspond to the evaporation and "addition" of water mass to the drying air. Under this system, the "material" side of the component is made up of a "fictional" water current, which would lose mass equal to the amount evaporated. The drying air would then have it's composition and mass recalculated with this added water mass ( I am currently not preoccupied with the enthalpy for the material side). The attached image is a schematic representation of what I intend for the custom component to do:

The attached code has been working reasonably well for this purpose, and has been able to calculate the heat load transferred in the drying process (to which I can subtract the minimum theoretical heat load required to obtain an estimation for the losses occuring in the dryer). However, I have run into issues when I attempt to attach this to a combustion system. As is typical for an industrial continuous dryer, the drying air is heated through combustion of a fuel (such as natural gas, or in the example provided, pure Methane). The problem comes when I create a system such as the one present in the following figure:

This results in either a jacobian error, or an overcharacterization of the system. What's confounding me is that when I simply perform the calculations for the combustion separately, and feed the resulting exhaust stream as parameters using the set_attr method for the dryer inlet, it functions and is able to converge. The problem only comes when the entire network is connected, and no combination of parameters has been able to allow it to converge.

Is there anything I am doing wrong with the component definition, or the equations present therein that could generate this problem? It would be highly beneficial for the dryer to be represented as a single subnetwork.

Kind regards and a preemptive thanks,
Manuel

class dryer(Component):
    r"""
    Class for Drying Component

    The component dryer is loosely based on the existing Heat Exchanger component;
    The values for the second side are dummy values and represent the water content trapped 
    within the material we are drying.


    **Mandatory Equations**

    - :py:meth:`tespy.components.component.Component.fluid_func`
    - :py:meth:`tespy.components.component.Component.mass_flow_func`
    - :py:meth:`tespy.components.heat_exchangers.base.HeatExchanger.energy_balance_func`

    **Optional Equations**

    - :py:meth:`tespy.components.heat_exchangers.base.HeatExchanger.energy_balance_hot_func`
    - :py:meth:`tespy.components.customs.dryer.dryer.qevap_func
    - Functions are still missing WIP 

    Inlets/Outlets

    - in1, in2 (index 1: air side, index 2: "fictional" material side)
    - out1, out2 (index 1: air side, index 2: "fictional" material side)

    Image
// 
    Parameters
    ----------
    label : str
        The label of the component.

    design : list
        List containing design parameters (stated as String).

    offdesign : list
        List containing offdesign parameters (stated as String).

    design_path : str
        Path to the components design case.

    local_offdesign : boolean
        Treat this component in offdesign mode in a design calculation.

    local_design : boolean
        Treat this component in design mode in an offdesign calculation.

    char_warnings : boolean
        Ignore warnings on default characteristics usage for this component.

    printout : boolean
        Include this component in the network's results printout.

    Q : float, dict
        Heat transfer, :math:`Q/\text{W}`.

    Losses  : float, dict, 
        Heat transfer that wasn't useful work

    Qevap : float, dict, :code:`"var"`
        Heat load of evaporation of water
        
    mevap : float, dict, :code:`"var"`
        Mass of water evaporated from material
        
    Qheatin: float, dict 
        Heat Transfer referring to material heating
    
    Mmat
        Mass of  material being dried
    
    Dt
        temperature difference between material at inlet and outlet
    
    cp
        Specific heat capacity of material being dried.
    

    
    Note
    ----
    

    Example
    -------
    A ceramic current of 500 kg/hour is being dried by an air current, which enters at 450 ºC and exits at 75ºC.
    We define the quantity of material, the amount of water evaporated and the air composition at the inlet and outlet.
    
    
    >>>from tespy.networks import Network
    >>>from tespy.connections import Connection, Ref, Bus
    >>>from tespy.components.customs import dryer

    >>>from tespy.components import Sink,Source
    
    >>>nw = Network(p_unit="bar", T_unit="C",h_unit='J / kg',)

    >>>dry=dryer.dryer('dryer',)

    >>>s1=Source('s1')
    >>>s2=Source('s2')
    >>>s3=Sink('s3')
    >>>s4=Sink('s4')
    
    
    >>>c1=Connection(s1,"out1", dry,'in1',label='1')
    >>>c3=Connection(s2,"out1", dry, 'in2', label ="3")

    >>>c2=Connection(dry,"out1", s3,'in1',label='2')
    >>>c4=Connection(dry,'out2', s4, 'in1',label="4")

    >>>nw.add_conns(c1,c2,c3,c4)

    >>>xH2O=0.02
    >>>xH2Oe=0.04
    
    dry.set_attr(
        Mmat=500/3600, cp = 1000, Dt = 40, mevap = (500/3600)*(0.04), 
        )

    c1.set_attr(
        T=450, m=0.02885,p=1,
        fluid={'O2': 0.192585 , 'N2': 0.748354, "Ar": 0.012781,  'H2O': 0.020655,'CO2':0.02562
    },
        )
    c2.set_attr(
        p=1, m0=0.9, T=75,
        fluid={'O2': 0.192585 * (1 - xH2Oe), 'N2': 0.748354 * (1 - xH2Oe), "Ar": 0.012781 * (1 - xH2Oe), 'H2O': xH2Oe,'CO2': 0.02562 * (1-xH2Oe)}
        ),
        
    c3.set_attr( 
        T=20,p=1,
        fluid={'H2O': 1},   
        )

    c4.set_attr(
        p=1,m=0.089,
        fluid={'H2O': 1}
        )
    
    
    nw2.solve(mode='design')
    nw2.print_results()

    """

    @staticmethod
    def component():
        return 'dryer'

    def get_parameters(self):
        return {
            'Q': dc_cp(
                max_val=0, func=self.energy_balance_hot_func, num_eq=1,
                deriv=self.energy_balance_hot_deriv,
                latex=self.energy_balance_hot_func_doc),
       
            'Losses': dc_cp(
                min_val = 0, max_val = 1e10, num_eq=1,
                deriv = self.loss_func_deriv, func = self.loss_func,
                latex = self.loss_func_doc),
        
                
            'Qevap': dc_cp(
                min_val = 0, max_val = 1e10, num_eq=1,
                deriv = self.qevap_func_deriv, func = self.qevap_func,
                 latex = self.qevap_func_doc),
            'mevap': dc_cp(
                min_val = 0, max_val = 1e10, d=1e-4,
               ),
             'Qheatin': dc_cp(
                 min_val = 0 , max_val = 1e10, num_eq=1 ,
                 deriv = self.qheatin_deriv, func = self.qheatin_func,
                  latex=self.qheatin_func_doc),
             'Mmat': dc_cp(
                 min_val = 0, max_val = 1e4, d = 1e-4),
             'Dt': dc_cp(
                 min_val = 0 , max_val = 5e3, d = 1e-4),
             'cp': dc_cp(
                 min_val = 0 , max_val = 1e5, d= 1e-4),
            
           
        }

    def get_mandatory_constraints(self):
        self.variable_fluids = set(
            [fluid for c in self.inl + self.outl for fluid in c.fluid.is_var]
        )
        num_fluid_eq = len(self.variable_fluids)
        if num_fluid_eq == 0:
            num_fluid_eq = 1
            self.variable_fluids = [list(self.inl[0].fluid.is_set)[0]]
        return {
            'energy_balance_constraints': {
                'func': self.energy_balance_func,
                'deriv': self.energy_balance_deriv,
                'constant_deriv': False, 'latex': self.energy_balance_func_doc,
                'num_eq': 1},
            'fluid_constraints': {
                'func': self.fluid_func, 'deriv': self.fluid_deriv,
                'constant_deriv': False, 'latex': self.fluid_func_doc,
                'num_eq': num_fluid_eq},
            'mass constraints': {
                'func':self.mass_balance_func, 'deriv': self.mass_balance_deriv,
                'constant_deriv': True, 'latex': self.mass_balance_doc,
                'num_eq': 1,}
        }

    @staticmethod
    def inlets():
        return ['in1', 'in2']

    @staticmethod
    def outlets():
        return ['out1', 'out2']
    
    @staticmethod
    def is_branch_source():
        return True

    def start_branch(self):
        branches = {}
        for outconn in self.outl:
            branch = {
                "connections": [outconn],
                "components": [self, outconn.target],
                "subbranches": {}
            }
            outconn.target.propagate_to_target(branch)

            branches[outconn.label] = branch
        return branches

    def propagate_to_target(self, branch):
        return

    def propagate_wrapper_to_target(self, branch):
        branch["components"] += [self]
        for outconn in self.outl:
            branch["connections"] += [outconn]
            outconn.target.propagate_wrapper_to_target(branch)

        
    def energy_balance_func(self):
        r"""
        Equation for dryer energy balance.

        Returns
        -------
        residual : float
            Residual value of equation.

            .. math::

                0 = \dot{m}_{in,1} \cdot \left(h_{out,1} - h_{in,1} \right) +
                \dot{m}_{in,2} \cdot \left(h_{out,2} - h_{in,2} \right)
        """
        return (
            self.inl[0].m.val_SI
            * (self.outl[0].h.val_SI - self.inl[0].h.val_SI)
            + self.inl[1].m.val_SI
            * (self.outl[1].h.val_SI - self.inl[1].h.val_SI)
        )

    def energy_balance_func_doc(self, label):
        r"""
        Equation for dryer energy balance.

        Parameters
        ----------
        label : str
            Label for equation.

        Returns
        -------
        latex : str
            LaTeX code of equations applied.
        """
        latex = (
            r'0 = \dot{m}_\mathrm{in,1} \cdot \left(h_\mathrm{out,1} -'
            r' h_\mathrm{in,1} \right) +\dot{m}_\mathrm{in,2} \cdot '
            r'\left(h_\mathrm{out,2} - h_\mathrm{in,2} \right)')
        return generate_latex_eq(self, latex, label)
    
    def energy_balance_deriv(self, increment_filter, k):
        r"""
        Partial derivatives of energy balance function.

        Parameters
        ----------
        increment_filter : ndarray
            Matrix for filtering non-changing variables.

        k : int
            Position of derivatives in Jacobian matrix (k-th equation).
        """
        for _c_num, i in enumerate(self.inl):
            o = self.outl[_c_num]
            if self.is_variable(i.m, increment_filter):
                self.jacobian[k, i.m.J_col] = o.h.val_SI - i.h.val_SI
            if self.is_variable(i.h, increment_filter):
                self.jacobian[k, i.h.J_col] = -i.m.val_SI
            if self.is_variable(o.h, increment_filter):
                self.jacobian[k, o.h.J_col] = i.m.val_SI
    def loss_func(self):
        r"""
        calculating thermal losses in dryer .

        Returns
        -------
        residual : float
            Residual value of equation.

            .. math::

                0 = \dot{m}_{in,1} \cdot \left(h_{out,1} - h_{in,1} \right) +
                \dot{m}_{in,2} \cdot \left(h_{out,2} - h_{in,2} \right)
        """
        return (
            self.inl[0].m.val_SI
            * (self.outl[0].h.val_SI - self.inl[0].h.val_SI)
            - self.Losses.val
        )
    def loss_func_deriv(self,increment_filter,k):
        r"""
        Partial derivatives for dryer heat load.

        Parameters
        ----------
        increment_filter : ndarray 
            Matrix for filtering non-changing variables.

        k : int
            Position of derivatives in Jacobian matrix (k-th equation).
        """
        i = self.inl[0]
        o = self.outl[0]
        if self.is_variable(i.m):
            self.jacobian[k, i.m.J_col] = o.h.val_SI - i.h.val_SI
        if self.is_variable(i.h):
            self.jacobian[k, i.h.J_col] = -i.m.val_SI
        if self.is_variable(o.h):
            self.jacobian[k, o.h.J_col] = i.m.val_SI
    def loss_func_doc(self, label):
        r"""
        Equation for hot side heat exchanger energy balance.

        Parameters
        ----------
        label : str
            Label for equation.

        Returns
        -------
        latex : str
            LaTeX code of equations applied.
        """ 
        latex = (
            r'0 =\dot{m}_{in,1} \cdot \left(h_{out,1}-'
            r'h_{in,1}\right)-\dot{Q}')
        return generate_latex_eq(self, latex, label)
    
    def mass_balance_func(self):
        r"""
        Equation for heat exchanger mass balance.

        Returns
        -------
        residual : float
            Residual value of equation.

            .. math::

                0 = \dot{m}_{in,1} \cdot \left(h_{out,1} - h_{in,1} \right) +
                \dot{m}_{in,2} \cdot \left(h_{out,2} - h_{in,2} \right)
        """
        return (
            self.inl[0].m.val_SI + self.inl[1].m.val_SI
            - (self.outl[0].m.val_SI + self.outl[1].m.val_SI)
                
            
            
        )

    
    def mass_balance_deriv(self, k):
        r"""
        Partial derivatives of energy balance function.

        Parameters
        ----------
        increment_filter : ndarray
            Matrix for filtering non-changing variables.

        k : int
            Position of derivatives in Jacobian matrix (k-th equation).
        """
        for _c_num, i in enumerate(self.inl):
            o = self.outl[_c_num]
            if self.is_variable(i.m):
                self.jacobian[k, i.m.J_col] =  1 
            if self.is_variable(o.m):
                self.jacobian[k, o.m.J_col] = -1
    
    def mass_balance_doc(self, label):
        r"""
        Equation for heat exchanger energy balance.

        Parameters
        ----------
        label : str
            Label for equation.

        Returns
        -------
        latex : str
            LaTeX code of equations applied.
        """
        latex = (
            r'0 = \dot{m}_\mathrm{in,1} \cdot \left(h_\mathrm{out,1} -'
            r' h_\mathrm{in,1} \right) +\dot{m}_\mathrm{in,2} \cdot '
            r'\left(h_\mathrm{out,2} - h_\mathrm{in,2} \right)')
        return generate_latex_eq(self, latex, label)

  
                
    def fluid_func(self):
        r'''
        Equation for determination of mass and mass fractions of fluid

        Returns
        -------
        Residuals of the fluid function

        '''
        i1=self.inl[0]
        i2=self.inl[1]
        o1=self.outl[0]
        o2=self.outl[1]
        residual = []
        for fluid in self.variable_fluids:
            res = i1.fluid.val[fluid] * i1.m.val_SI - o1.fluid.val[fluid] * o1.m.val_SI +  i2.fluid.val[fluid] * i2.m.val_SI - o2.fluid.val[fluid] * o2.m.val_SI
            residual += [res]
            
            
            
        
        return(residual)
    def fluid_deriv(self,increment_filter,k):
        r'''
        Derivatives for the determination of mass and mass fractions of fluid

        Parameters
        ----------
        increment_filter : ndarray
            Matrix for filtering non-changing variables.

        k : int
            Position of derivatives in Jacobian matrix (k-th equation).
        """
        Returns
        -------
        The jacobian for the residual functions

        '''
        i1=self.inl[0]
        i2=self.inl[1]
        o1=self.outl[0]
        o2=self.outl[1]
        for fluid in self.variable_fluids:
            
            if self.is_variable(i1.m):
                self.jacobian[k, i1.m.J_col] = i1.fluid.val[fluid]
            if fluid in i1.fluid.is_var:
                self.jacobian[k, i1.fluid.J_col[fluid]] =  i1.m.val_SI
            if self.is_variable(i2.m):
                self.jacobian[k ,i2.m.J_col] = i2.fluid.val[fluid]
            if fluid in i2.fluid.is_var:
                self.jacobian[k, i2.fluid.J_col[fluid]] =  i2.m.val_SI
            if self.is_variable(o1.m):
                self.jacobian[k,o1.m.J_col] = -o1.fluid.val[fluid]
            if fluid in o1.fluid.is_var:
                self.jacobian[k, o1.fluid.J_col[fluid]] = - o1.m.val_SI
            if self.is_variable(o2.m):
                self.jacobian[k,o2.m.J_col] = -o2.fluid.val[fluid]
            if fluid in o2.fluid.is_var:
                self.jacobian[k, o2.fluid.J_col[fluid]] = - o2.m.val_SI
        return()
    def fluid_func_doc(self,label):
        
        return()
    
    def energy_balance_hot_func(self):
        r"""
        Equation for heat transferred in the dryer (Based on specific enthalpies of entering and exiting drying air)

        Returns
        -------
        residual : float
            Residual value of equation.

            .. math::

                0 =\dot{m}_{in,1} \cdot \left(h_{out,1}-h_{in,1}\right)-\dot{Q}
        """
        return self.inl[0].m.val_SI * (
            self.outl[0].h.val_SI - self.inl[0].h.val_SI
        )  - self.Q.val

    def energy_balance_hot_func_doc(self, label):
        r"""
        Equation for heat transferred in the dryer (Based on specific enthalpies of entering and exiting drying air)

        Parameters
        ----------
        label : str
            Label for equation.

        Returns
        -------
        latex : str
            LaTeX code of equations applied.
        """
        latex = (
            r'0 =\dot{m}_{in,1} \cdot \left(h_{out,1}-'
            r'h_{in,1}\right)-\dot{Q}')
        return generate_latex_eq(self, latex, label)

    def energy_balance_hot_deriv(self, increment_filter, k):
        r"""
        Partial derivatives for dryer heat transfer based on entering and exiting enthalpies

        Parameters
        ----------
        increment_filter : ndarray
            Matrix for filtering non-changing variables.

        k : int
            Position of derivatives in Jacobian matrix (k-th equation).
        """
        i = self.inl[0]
        o = self.outl[0]
        if self.is_variable(i.m):
            self.jacobian[k, i.m.J_col] = o.h.val_SI - i.h.val_SI
        if self.is_variable(i.h):
            self.jacobian[k, i.h.J_col] = -i.m.val_SI
        if self.is_variable(o.h):
            self.jacobian[k, o.h.J_col] = i.m.val_SI

    
    def qevap_func(self):
        r"""
        calculating heat required to evaporate water.
    
        Returns
        -------
        residual : float
            Residual value of equation.
    
            .. math::
    
                0 = \dot{m}_{in,1} \cdot \left(h_{out,1} - h_{in,1} \right) +
                \dot{m}_{in,2} \cdot \left(h_{out,2} - h_{in,2} \right)
        """
        self.inl[0].m.val_SI * (self.outl[0].h.val_SI - self.inl[0].h.val_SI)
        -(self.mevap.val*2273*1e3)-self.qevap.val
        
    def qevap_func_deriv(self,increment_filter,k):
        r"""
        Partial derivatives for water evaporation.
        falta definir corretamente as derivadas para cada uma das cosias que podem ser variaveis
    
        Parameters
        ----------
        increment_filter : ndarray 
            Matrix for filtering non-changing variables.
    
        k : int
            Position of derivatives in Jacobian matrix (k-th equation).
        """
        if self.is_variable(self.inl[0].m):
            self.jacobian[k,self.inl[0].m.J_col] = self.outl[0].h.val_SI - self.inl[0].h.val_SI
        if self.is_variable(self.mevap):
            self.jacobian[k,self.mevap.J_col] = 2273*1e3
            
        
    def qevap_func_doc(self, label):
        r"""
        Equation for heat of water evaporation.
    
        Parameters
        ----------
        label : str
            Label for equation.
    
        Returns
        -------
        latex : str
            LaTeX code of equations applied.
        """
        latex = (
            r'0 =\dot{m}_{in,1} \cdot \left(h_{out,1}-'
            r'h_{in,1}\right)-\dot{Q}')
        return generate_latex_eq(self, latex, label)
    def qheatin_func(self):
        r"""
        Calculates Heat used for material heating. 
    
        Returns
        -------
        residual : float
            Residual value of equation.
    
            .. math::
    
                0 = \dot{m}_{in,1} \cdot \left(h_{out,1} - h_{in,1} \right) +
                \dot{m}_{in,2} \cdot \left(h_{out,2} - h_{in,2} \right)
        """
        return (
            self.Mmat.val*self.cp.val*self.Dt.val
        )
    def qheatin_deriv(self,increment_filter,k):
        r"""
        Derivatives for heat in material heating
        Falta uma estrutura do tipo, if 1 deles variavel = multiplicação dos outros 2 
        Parameters
        ----------
        increment_filter : ndarray 
            Matrix for filtering non-changing variables.
    
        k : int
            Position of derivatives in Jacobian matrix (k-th equation).
        """
        if self.is_variable(self.Mmat.val):
            self.jacobian[k,self.Mmat.J_col] = self.cp.val*self.Dt.val
        if self.is_variable(self.cp.val):
            self.jacobian[k,self.cp.J_col] = self.Mmat.val*self.Dt.val
        if self.is_variable(self.Dt.val):
            self.jacobian[k,self.Dt.J_col] = self.cp.val*self.Mmat.val
              
    def qheatin_func_doc(self, label):
        r"""
        Equation for material heat load.
    
        Parameters
        ----------
        label : str
            Label for equation.
    
        Returns
        -------
        latex : str
            LaTeX code of equations applied.
        """
        latex = (
            r'0 =\dot{m}_{in,1} \cdot \left(h_{out,1}-'
            r'h_{in,1}\right)-\dot{Q}')
        return generate_latex_eq(self, latex, label)
    def bus_func(self, bus):
        r"""
        Calculate the value of the bus function.

        Parameters
        ----------
        bus : tespy.connections.bus.Bus
            TESPy bus object.

        Returns
        -------
        val : float
            Value of energy transfer :math:`\dot{E}`. This value is passed to
            :py:meth:`tespy.components.component.Component.calc_bus_value`
            for value manipulation according to the specified characteristic
            line of the bus.

            .. math::

                \dot{E} = \dot{m}_{in,1} \cdot \left(
                h_{out,1} - h_{in,1} \right)
        """
        return self.inl[0].m.val_SI * (
            self.outl[0].h.val_SI - self.inl[0].h.val_SI
        )

    def bus_func_doc(self, bus):
        r"""
        Return LaTeX string of the bus function.

        Parameters
        ----------
        bus : tespy.connections.bus.Bus
            TESPy bus object.

        Returns
        -------
        latex : str
            LaTeX string of bus function.
        """
        return (
            r'\dot{m}_\mathrm{in,1} \cdot \left(h_\mathrm{out,1} - '
            r'h_\mathrm{in,1} \right)')

    def bus_deriv(self, bus):
        r"""
        Calculate partial derivatives of the bus function.

        Parameters
        ----------
        bus : tespy.connections.bus.Bus
            TESPy bus object.

        Returns
        -------
        deriv : ndarray
            Matrix of partial derivatives.
        """
        f = self.calc_bus_value
        if self.inl[0].m.is_var:
            if self.inl[0].m.J_col not in bus.jacobian:
                bus.jacobian[self.inl[0].m.J_col] = 0
            bus.jacobian[self.inl[0].m.J_col] -= self.numeric_deriv(f, 'm', self.inl[0], bus=bus)

        if self.inl[0].h.is_var:
            if self.inl[0].h.J_col not in bus.jacobian:
                bus.jacobian[self.inl[0].h.J_col] = 0
            bus.jacobian[self.inl[0].h.J_col] -= self.numeric_deriv(f, 'h', self.inl[0], bus=bus)

        if self.outl[0].h.is_var:
            if self.outl[0].h.J_col not in bus.jacobian:
                bus.jacobian[self.outl[0].h.J_col] = 0
            bus.jacobian[self.outl[0].h.J_col] -= self.numeric_deriv(f, 'h', self.outl[0], bus=bus)

    def initialise_source(self, c, key):
        r"""
        Return a starting value for pressure and enthalpy at outlet.

        Parameters
        ----------
        c : tespy.connections.connection.Connection
            Connection to perform initialisation on.

        key : str
            Fluid property to retrieve.

        Returns
        -------
        val : float
            Starting value for pressure/enthalpy in SI units.

            .. math::

                val = \begin{cases}
                4 \cdot 10^5 & \text{key = 'p'}\\
                h\left(p, 200 \text{K} \right) & \text{key = 'h' at outlet 1}\\
                h\left(p, 250 \text{K} \right) & \text{key = 'h' at outlet 2}
                \end{cases}
        """
        if key == 'p':
            return 50e5
        elif key == 'h':
            if c.source_id == 'out1':
                T = 200 + 273.15
            else:
                T = 250 + 273.15
            return h_mix_pT(c.p.val_SI, T, c.fluid_data, c.mixing_rule)

    def initialise_target(self, c, key):
        r"""
        Return a starting value for pressure and enthalpy at inlet.

        Parameters
        ----------
        c : tespy.connections.connection.Connection
            Connection to perform initialisation on.

        key : str
            Fluid property to retrieve.

        Returns
        -------
        val : float
            Starting value for pressure/enthalpy in SI units.

            .. math::

                val = \begin{cases}
                4 \cdot 10^5 & \text{key = 'p'}\\
                h\left(p, 300 \text{K} \right) & \text{key = 'h' at inlet 1}\\
                h\left(p, 220 \text{K} \right) & \text{key = 'h' at outlet 2}
                \end{cases}
        """
        if key == 'p':
            return 50e5
        elif key == 'h':
            if c.target_id == 'in1':
                T = 300 + 273.15
            else:
                T = 220 + 273.15
            return h_mix_pT(c.p.val_SI, T, c.fluid_data, c.mixing_rule)

    def calc_parameters(self):
        r"""Postprocessing parameter calculation."""
        # component parameters
        self.Q.val = self.inl[0].m.val_SI * (
            self.outl[0].h.val_SI - self.inl[0].h.val_SI)
        
        self.Qevap.val = self.mevap.val * 2273 * 1e3 
        
        self.Qheatin.val = self.Mmat.val * self.Dt.val * self.cp.val

        self.Losses.val = self.Q.val +self.Qheatin.val + self.Qevap.val

[fwitte]

Hi @MMArro,

thanks for reaching out. That is a lot to process so have a little patience with me. I would actually like to invite you to the next online user meeting, so we can discuss this in person?

Best

Francesco

[fwitte]

Great, looking forward to meeting you!

[risinggard]

Hi!
I’m looking into the feasibility of using TESPy for dryer modeling. Do you think you could share your code that fixes the initial issue, @MMArro ?

[MMArro]

Hi! I’m looking into the feasibility of using TESPy for dryer modeling. Do you think you could share your code that fixes the initial issue, @MMArro ?

Hello, sure, no problem. Please do take into account that this does not accurately model physical systems for which you have real data. I would advise including a sink/source to include any potential mass losses/gains occuring from bad insulation, as these are often not airtight systems.

class dryer(Component):
r"""
Class for Drying Component

The component dryer is loosely based on the existing Heat Exchanger component;
The values for the second side are dummy values and represent the water content trapped 
within the material we are drying.


**Mandatory Equations**

- :py:meth:`tespy.components.component.Component.fluid_func`
- :py:meth:`tespy.components.component.Component.mass_flow_func`
- :py:meth:`tespy.components.heat_exchangers.base.HeatExchanger.energy_balance_func`

**Optional Equations**

- :py:meth:`tespy.components.heat_exchangers.base.HeatExchanger.energy_balance_hot_func`
- :py:meth:`tespy.components.customs.dryer.dryer.qevap_func
- Functions are still missing WIP 

Inlets/Outlets

- in1, in2 (index 1: air side, index 2: "fictional" material side)
- out1, out2 (index 1: air side, index 2: "fictional" material side)

Image

//
Parameters
----------
label : str
The label of the component.

design : list
    List containing design parameters (stated as String).

offdesign : list
    List containing offdesign parameters (stated as String).

design_path : str
    Path to the components design case.

local_offdesign : boolean
    Treat this component in offdesign mode in a design calculation.

local_design : boolean
    Treat this component in design mode in an offdesign calculation.

char_warnings : boolean
    Ignore warnings on default characteristics usage for this component.

printout : boolean
    Include this component in the network's results printout.

Q : float, dict
    Heat transfer, :math:`Q/\text{W}`.

Losses  : float, dict, 
    Heat transfer that wasn't useful work

Qevap : float, dict, :code:`"var"`
    Heat load of evaporation of water
    
mevap : float, dict, :code:`"var"`
    Mass of water evaporated from material
    
Qheatin: float, dict 
    Heat Transfer referring to material heating

Mmat
    Mass of  material being dried

Dt
    temperature difference between material at inlet and outlet

cp
    Specific heat capacity of material being dried.



Note
----


Example
-------
A ceramic current of 500 kg/hour is being dried by an air current, which enters at 450 ºC and exits at 75ºC.
We define the quantity of material, the amount of water evaporated and the air composition at the inlet and outlet.


>>>from tespy.networks import Network
>>>from tespy.connections import Connection, Ref, Bus
>>>from tespy.components.customs import dryer

>>>from tespy.components import Sink,Source

>>>nw = Network(p_unit="bar", T_unit="C",h_unit='J / kg',)

>>>dry=dryer.dryer('dryer',)

>>>s1=Source('s1')
>>>s2=Source('s2')
>>>s3=Sink('s3')
>>>s4=Sink('s4')


>>>c1=Connection(s1,"out1", dry,'in1',label='1')
>>>c3=Connection(s2,"out1", dry, 'in2', label ="3")

>>>c2=Connection(dry,"out1", s3,'in1',label='2')
>>>c4=Connection(dry,'out2', s4, 'in1',label="4")

>>>nw.add_conns(c1,c2,c3,c4)

>>>xH2O=0.02
>>>xH2Oe=0.04

dry.set_attr(
    Mmat=500/3600, cp = 1000, Dt = 40, mevap = (500/3600)*(0.04), 
    )

c1.set_attr(
    T=450, m=0.02885,p=1,
    fluid={'O2': 0.192585 , 'N2': 0.748354, "Ar": 0.012781,  'H2O': 0.020655,'CO2':0.02562
},
    )
c2.set_attr(
    p=1, m0=0.9, T=75,
    fluid={'O2': 0.192585 * (1 - xH2Oe), 'N2': 0.748354 * (1 - xH2Oe), "Ar": 0.012781 * (1 - xH2Oe), 'H2O': xH2Oe,'CO2': 0.02562 * (1-xH2Oe)}
    ),
    
c3.set_attr( 
    T=20,p=1,
    fluid={'H2O': 1},   
    )

c4.set_attr(
    p=1,m=0.089,
    fluid={'H2O': 1}
    )


nw2.solve(mode='design')
nw2.print_results()

"""

@staticmethod
def component():
    return 'dryer'

def get_parameters(self):
    return {
        'Q': dc_cp(
            max_val=0, func=self.energy_balance_hot_func, num_eq=1,
            deriv=self.energy_balance_hot_deriv,
            latex=self.energy_balance_hot_func_doc),
   
        'Losses': dc_cp(
            min_val = 0, max_val = 1e10, num_eq=1,
            deriv = self.loss_func_deriv, func = self.loss_func,
            latex = self.loss_func_doc),
    
            
        'Qevap': dc_cp(
            min_val = 0, max_val = 1e10, num_eq=1,
            deriv = self.qevap_func_deriv, func = self.qevap_func,
             latex = self.qevap_func_doc),
        'mevap': dc_cp(
            min_val = 0, max_val = 1e10, d=1e-4,
           ),
         'Qheatin': dc_cp(
             min_val = 0 , max_val = 1e10, num_eq=1 ,
             deriv = self.qheatin_deriv, func = self.qheatin_func,
              latex=self.qheatin_func_doc),
         'Mmat': dc_cp(
             min_val = 0, max_val = 1e4, d = 1e-4),
         'Dt': dc_cp(
             min_val = 0 , max_val = 5e3, d = 1e-4),
         'cp': dc_cp(
             min_val = 0 , max_val = 1e5, d= 1e-4),
        
       
    }

def get_mandatory_constraints(self):
    self.variable_fluids = set(
        [fluid for c in self.inl + self.outl for fluid in c.fluid.is_var]
    )
    num_fluid_eq = len(self.variable_fluids)
    if num_fluid_eq == 0:
        num_fluid_eq = 1
        self.variable_fluids = [list(self.inl[0].fluid.is_set)[0]]
    return {
        'energy_balance_constraints': {
            'func': self.energy_balance_func,
            'deriv': self.energy_balance_deriv,
            'constant_deriv': False, 'latex': self.energy_balance_func_doc,
            'num_eq': 1},
        'fluid_constraints': {
            'func': self.fluid_func, 'deriv': self.fluid_deriv,
            'constant_deriv': False, 'latex': self.fluid_func_doc,
            'num_eq': num_fluid_eq},
        'mass constraints': {
            'func':self.mass_balance_func, 'deriv': self.mass_balance_deriv,
            'constant_deriv': True, 'latex': self.mass_balance_doc,
            'num_eq': 1,}
    }

@staticmethod
def inlets():
    return ['in1', 'in2']

@staticmethod
def outlets():
    return ['out1', 'out2']

@staticmethod
def is_branch_source():
    return True

def start_branch(self):
    branches = {}
    for outconn in self.outl:
        branch = {
            "connections": [outconn],
            "components": [self, outconn.target],
            "subbranches": {}
        }
        outconn.target.propagate_to_target(branch)

        branches[outconn.label] = branch
    return branches

def propagate_to_target(self, branch):
    return

def propagate_wrapper_to_target(self, branch):
    branch["components"] += [self]
    for outconn in self.outl:
        branch["connections"] += [outconn]
        outconn.target.propagate_wrapper_to_target(branch)

    
def energy_balance_func(self):
    r"""
    Equation for dryer energy balance.

    Returns
    -------
    residual : float
        Residual value of equation.

        .. math::

            0 = \dot{m}_{in,1} \cdot \left(h_{out,1} - h_{in,1} \right) +
            \dot{m}_{in,2} \cdot \left(h_{out,2} - h_{in,2} \right)
    """
    return (
        self.inl[0].m.val_SI
        * (self.outl[0].h.val_SI - self.inl[0].h.val_SI)
        + self.inl[1].m.val_SI
        * (self.outl[1].h.val_SI - self.inl[1].h.val_SI)
    )

def energy_balance_func_doc(self, label):
    r"""
    Equation for dryer energy balance.

    Parameters
    ----------
    label : str
        Label for equation.

    Returns
    -------
    latex : str
        LaTeX code of equations applied.
    """
    latex = (
        r'0 = \dot{m}_\mathrm{in,1} \cdot \left(h_\mathrm{out,1} -'
        r' h_\mathrm{in,1} \right) +\dot{m}_\mathrm{in,2} \cdot '
        r'\left(h_\mathrm{out,2} - h_\mathrm{in,2} \right)')
    return generate_latex_eq(self, latex, label)

def energy_balance_deriv(self, increment_filter, k):
    r"""
    Partial derivatives of energy balance function.

    Parameters
    ----------
    increment_filter : ndarray
        Matrix for filtering non-changing variables.

    k : int
        Position of derivatives in Jacobian matrix (k-th equation).
    """
    for _c_num, i in enumerate(self.inl):
        o = self.outl[_c_num]
        if self.is_variable(i.m, increment_filter):
            self.jacobian[k, i.m.J_col] = o.h.val_SI - i.h.val_SI
        if self.is_variable(i.h, increment_filter):
            self.jacobian[k, i.h.J_col] = -i.m.val_SI
        if self.is_variable(o.h, increment_filter):
            self.jacobian[k, o.h.J_col] = i.m.val_SI
def loss_func(self):
    r"""
    calculating thermal losses in dryer .

    Returns
    -------
    residual : float
        Residual value of equation.

        .. math::

            0 = \dot{m}_{in,1} \cdot \left(h_{out,1} - h_{in,1} \right) +
            \dot{m}_{in,2} \cdot \left(h_{out,2} - h_{in,2} \right)
    """
    return (
        self.inl[0].m.val_SI
        * (self.outl[0].h.val_SI - self.inl[0].h.val_SI)
        - self.Losses.val
    )
def loss_func_deriv(self,increment_filter,k):
    r"""
    Partial derivatives for dryer heat load.

    Parameters
    ----------
    increment_filter : ndarray 
        Matrix for filtering non-changing variables.

    k : int
        Position of derivatives in Jacobian matrix (k-th equation).
    """
    i = self.inl[0]
    o = self.outl[0]
    if self.is_variable(i.m):
        self.jacobian[k, i.m.J_col] = o.h.val_SI - i.h.val_SI
    if self.is_variable(i.h):
        self.jacobian[k, i.h.J_col] = -i.m.val_SI
    if self.is_variable(o.h):
        self.jacobian[k, o.h.J_col] = i.m.val_SI
def loss_func_doc(self, label):
    r"""
    Equation for hot side heat exchanger energy balance.

    Parameters
    ----------
    label : str
        Label for equation.

    Returns
    -------
    latex : str
        LaTeX code of equations applied.
    """ 
    latex = (
        r'0 =\dot{m}_{in,1} \cdot \left(h_{out,1}-'
        r'h_{in,1}\right)-\dot{Q}')
    return generate_latex_eq(self, latex, label)

def mass_balance_func(self):
    r"""
    Equation for heat exchanger mass balance.

    Returns
    -------
    residual : float
        Residual value of equation.

        .. math::

            0 = \dot{m}_{in,1} \cdot \left(h_{out,1} - h_{in,1} \right) +
            \dot{m}_{in,2} \cdot \left(h_{out,2} - h_{in,2} \right)
    """
    return (
        self.inl[0].m.val_SI + self.inl[1].m.val_SI
        - (self.outl[0].m.val_SI + self.outl[1].m.val_SI)
            
        
        
    )


def mass_balance_deriv(self, k):
    r"""
    Partial derivatives of energy balance function.

    Parameters
    ----------
    increment_filter : ndarray
        Matrix for filtering non-changing variables.

    k : int
        Position of derivatives in Jacobian matrix (k-th equation).
    """
    for _c_num, i in enumerate(self.inl):
        o = self.outl[_c_num]
        if self.is_variable(i.m):
            self.jacobian[k, i.m.J_col] =  1 
        if self.is_variable(o.m):
            self.jacobian[k, o.m.J_col] = -1

def mass_balance_doc(self, label):
    r"""
    Equation for heat exchanger energy balance.

    Parameters
    ----------
    label : str
        Label for equation.

    Returns
    -------
    latex : str
        LaTeX code of equations applied.
    """
    latex = (
        r'0 = \dot{m}_\mathrm{in,1} \cdot \left(h_\mathrm{out,1} -'
        r' h_\mathrm{in,1} \right) +\dot{m}_\mathrm{in,2} \cdot '
        r'\left(h_\mathrm{out,2} - h_\mathrm{in,2} \right)')
    return generate_latex_eq(self, latex, label)


            
def fluid_func(self):
    r'''
    Equation for determination of mass and mass fractions of fluid

    Returns
    -------
    Residuals of the fluid function

    '''
    i1=self.inl[0]
    i2=self.inl[1]
    o1=self.outl[0]
    o2=self.outl[1]
    residual = []
    for fluid in self.variable_fluids:
        res = i1.fluid.val[fluid] * i1.m.val_SI - o1.fluid.val[fluid] * o1.m.val_SI +  i2.fluid.val[fluid] * i2.m.val_SI - o2.fluid.val[fluid] * o2.m.val_SI
        residual += [res]
        
        
        
    
    return(residual)
def fluid_deriv(self,increment_filter,k):
    r'''
    Derivatives for the determination of mass and mass fractions of fluid

    Parameters
    ----------
    increment_filter : ndarray
        Matrix for filtering non-changing variables.

    k : int
        Position of derivatives in Jacobian matrix (k-th equation).
    """
    Returns
    -------
    The jacobian for the residual functions

    '''
    i1=self.inl[0]
    i2=self.inl[1]
    o1=self.outl[0]
    o2=self.outl[1]
    for fluid in self.variable_fluids:
        
        if self.is_variable(i1.m):
            self.jacobian[k, i1.m.J_col] = i1.fluid.val[fluid]
        if fluid in i1.fluid.is_var:
            self.jacobian[k, i1.fluid.J_col[fluid]] =  i1.m.val_SI
        if self.is_variable(i2.m):
            self.jacobian[k ,i2.m.J_col] = i2.fluid.val[fluid]
        if fluid in i2.fluid.is_var:
            self.jacobian[k, i2.fluid.J_col[fluid]] =  i2.m.val_SI
        if self.is_variable(o1.m):
            self.jacobian[k,o1.m.J_col] = -o1.fluid.val[fluid]
        if fluid in o1.fluid.is_var:
            self.jacobian[k, o1.fluid.J_col[fluid]] = - o1.m.val_SI
        if self.is_variable(o2.m):
            self.jacobian[k,o2.m.J_col] = -o2.fluid.val[fluid]
        if fluid in o2.fluid.is_var:
            self.jacobian[k, o2.fluid.J_col[fluid]] = - o2.m.val_SI
        k+=1
    return()
def fluid_func_doc(self,label):
    
    return()

def energy_balance_hot_func(self):
    r"""
    Equation for heat transferred in the dryer (Based on specific enthalpies of entering and exiting drying air)

    Returns
    -------
    residual : float
        Residual value of equation.

        .. math::

            0 =\dot{m}_{in,1} \cdot \left(h_{out,1}-h_{in,1}\right)-\dot{Q}
    """
    return self.inl[0].m.val_SI * (
        self.outl[0].h.val_SI - self.inl[0].h.val_SI
    )  - self.Q.val

def energy_balance_hot_func_doc(self, label):
    r"""
    Equation for heat transferred in the dryer (Based on specific enthalpies of entering and exiting drying air)

    Parameters
    ----------
    label : str
        Label for equation.

    Returns
    -------
    latex : str
        LaTeX code of equations applied.
    """
    latex = (
        r'0 =\dot{m}_{in,1} \cdot \left(h_{out,1}-'
        r'h_{in,1}\right)-\dot{Q}')
    return generate_latex_eq(self, latex, label)

def energy_balance_hot_deriv(self, increment_filter, k):
    r"""
    Partial derivatives for dryer heat transfer based on entering and exiting enthalpies

    Parameters
    ----------
    increment_filter : ndarray
        Matrix for filtering non-changing variables.

    k : int
        Position of derivatives in Jacobian matrix (k-th equation).
    """
    i = self.inl[0]
    o = self.outl[0]
    if self.is_variable(i.m):
        self.jacobian[k, i.m.J_col] = o.h.val_SI - i.h.val_SI
    if self.is_variable(i.h):
        self.jacobian[k, i.h.J_col] = -i.m.val_SI
    if self.is_variable(o.h):
        self.jacobian[k, o.h.J_col] = i.m.val_SI


def qevap_func(self):
    r"""
    calculating heat required to evaporate water.

    Returns
    -------
    residual : float
        Residual value of equation.

        .. math::

            0 = \dot{m}_{in,1} \cdot \left(h_{out,1} - h_{in,1} \right) +
            \dot{m}_{in,2} \cdot \left(h_{out,2} - h_{in,2} \right)
    """
    self.inl[0].m.val_SI * (self.outl[0].h.val_SI - self.inl[0].h.val_SI)
    -(self.mevap.val*2273*1e3)-self.qevap.val
    
def qevap_func_deriv(self,increment_filter,k):
    r"""
    Partial derivatives for water evaporation.
    falta definir corretamente as derivadas para cada uma das cosias que podem ser variaveis

    Parameters
    ----------
    increment_filter : ndarray 
        Matrix for filtering non-changing variables.

    k : int
        Position of derivatives in Jacobian matrix (k-th equation).
    """
    if self.is_variable(self.inl[0].m):
        self.jacobian[k,self.inl[0].m.J_col] = self.outl[0].h.val_SI - self.inl[0].h.val_SI
    if self.is_variable(self.mevap):
        self.jacobian[k,self.mevap.J_col] = 2273*1e3
        
    
def qevap_func_doc(self, label):
    r"""
    Equation for heat of water evaporation.

    Parameters
    ----------
    label : str
        Label for equation.

    Returns
    -------
    latex : str
        LaTeX code of equations applied.
    """
    latex = (
        r'0 =\dot{m}_{in,1} \cdot \left(h_{out,1}-'
        r'h_{in,1}\right)-\dot{Q}')
    return generate_latex_eq(self, latex, label)
def qheatin_func(self):
    r"""
    Calculates Heat used for material heating. 

    Returns
    -------
    residual : float
        Residual value of equation.

        .. math::

            0 = \dot{m}_{in,1} \cdot \left(h_{out,1} - h_{in,1} \right) +
            \dot{m}_{in,2} \cdot \left(h_{out,2} - h_{in,2} \right)
    """
    return (
        self.Mmat.val*self.cp.val*self.Dt.val
    )
def qheatin_deriv(self,increment_filter,k):
    r"""
    Derivatives for heat in material heating
    Falta uma estrutura do tipo, if 1 deles variavel = multiplicação dos outros 2 
    Parameters
    ----------
    increment_filter : ndarray 
        Matrix for filtering non-changing variables.

    k : int
        Position of derivatives in Jacobian matrix (k-th equation).
    """
    if self.is_variable(self.Mmat.val):
        self.jacobian[k,self.Mmat.J_col] = self.cp.val*self.Dt.val
    if self.is_variable(self.cp.val):
        self.jacobian[k,self.cp.J_col] = self.Mmat.val*self.Dt.val
    if self.is_variable(self.Dt.val):
        self.jacobian[k,self.Dt.J_col] = self.cp.val*self.Mmat.val
          
def qheatin_func_doc(self, label):
    r"""
    Equation for material heat load.

    Parameters
    ----------
    label : str
        Label for equation.

    Returns
    -------
    latex : str
        LaTeX code of equations applied.
    """
    latex = (
        r'0 =\dot{m}_{in,1} \cdot \left(h_{out,1}-'
        r'h_{in,1}\right)-\dot{Q}')
    return generate_latex_eq(self, latex, label)
def bus_func(self, bus):
    r"""
    Calculate the value of the bus function.

    Parameters
    ----------
    bus : tespy.connections.bus.Bus
        TESPy bus object.

    Returns
    -------
    val : float
        Value of energy transfer :math:`\dot{E}`. This value is passed to
        :py:meth:`tespy.components.component.Component.calc_bus_value`
        for value manipulation according to the specified characteristic
        line of the bus.

        .. math::

            \dot{E} = \dot{m}_{in,1} \cdot \left(
            h_{out,1} - h_{in,1} \right)
    """
    return self.inl[0].m.val_SI * (
        self.outl[0].h.val_SI - self.inl[0].h.val_SI
    )

def bus_func_doc(self, bus):
    r"""
    Return LaTeX string of the bus function.

    Parameters
    ----------
    bus : tespy.connections.bus.Bus
        TESPy bus object.

    Returns
    -------
    latex : str
        LaTeX string of bus function.
    """
    return (
        r'\dot{m}_\mathrm{in,1} \cdot \left(h_\mathrm{out,1} - '
        r'h_\mathrm{in,1} \right)')

def bus_deriv(self, bus):
    r"""
    Calculate partial derivatives of the bus function.

    Parameters
    ----------
    bus : tespy.connections.bus.Bus
        TESPy bus object.

    Returns
    -------
    deriv : ndarray
        Matrix of partial derivatives.
    """
    f = self.calc_bus_value
    if self.inl[0].m.is_var:
        if self.inl[0].m.J_col not in bus.jacobian:
            bus.jacobian[self.inl[0].m.J_col] = 0
        bus.jacobian[self.inl[0].m.J_col] -= self.numeric_deriv(f, 'm', self.inl[0], bus=bus)

    if self.inl[0].h.is_var:
        if self.inl[0].h.J_col not in bus.jacobian:
            bus.jacobian[self.inl[0].h.J_col] = 0
        bus.jacobian[self.inl[0].h.J_col] -= self.numeric_deriv(f, 'h', self.inl[0], bus=bus)

    if self.outl[0].h.is_var:
        if self.outl[0].h.J_col not in bus.jacobian:
            bus.jacobian[self.outl[0].h.J_col] = 0
        bus.jacobian[self.outl[0].h.J_col] -= self.numeric_deriv(f, 'h', self.outl[0], bus=bus)

def initialise_source(self, c, key):
    r"""
    Return a starting value for pressure and enthalpy at outlet.

    Parameters
    ----------
    c : tespy.connections.connection.Connection
        Connection to perform initialisation on.

    key : str
        Fluid property to retrieve.

    Returns
    -------
    val : float
        Starting value for pressure/enthalpy in SI units.

        .. math::

            val = \begin{cases}
            4 \cdot 10^5 & \text{key = 'p'}\\
            h\left(p, 200 \text{K} \right) & \text{key = 'h' at outlet 1}\\
            h\left(p, 250 \text{K} \right) & \text{key = 'h' at outlet 2}
            \end{cases}
    """
    if key == 'p':
        return 50e5
    elif key == 'h':
        if c.source_id == 'out1':
            T = 200 + 273.15
        else:
            T = 250 + 273.15
        return h_mix_pT(c.p.val_SI, T, c.fluid_data, c.mixing_rule)

def initialise_target(self, c, key):
    r"""
    Return a starting value for pressure and enthalpy at inlet.

    Parameters
    ----------
    c : tespy.connections.connection.Connection
        Connection to perform initialisation on.

    key : str
        Fluid property to retrieve.

    Returns
    -------
    val : float
        Starting value for pressure/enthalpy in SI units.

        .. math::

            val = \begin{cases}
            4 \cdot 10^5 & \text{key = 'p'}\\
            h\left(p, 300 \text{K} \right) & \text{key = 'h' at inlet 1}\\
            h\left(p, 220 \text{K} \right) & \text{key = 'h' at outlet 2}
            \end{cases}
    """
    if key == 'p':
        return 50e5
    elif key == 'h':
        if c.target_id == 'in1':
            T = 300 + 273.15
        else:
            T = 220 + 273.15
        return h_mix_pT(c.p.val_SI, T, c.fluid_data, c.mixing_rule)

def calc_parameters(self):
    r"""Postprocessing parameter calculation."""
    # component parameters
    self.Q.val = self.inl[0].m.val_SI * (
        self.outl[0].h.val_SI - self.inl[0].h.val_SI)
    
    self.Qevap.val = self.mevap.val * 2273 * 1e3 
    
    self.Qheatin.val = self.Mmat.val * self.Dt.val * self.cp.val

    self.Losses.val = self.Q.val +self.Qheatin.val + self.Qevap.val

[risinggard]

Thanks! I’ll keep working on this to see if I can adapt it to suit my needs :)

[MMArro]

Thanks! I’ll keep working on this to see if I can adapt it to suit my needs :)

If you need any help or wanna discuss this further at any point, feel free to ask or shoot me an email,

Thanks for showing interest!