Dryer modelling

[MMArro]

Hello!
I have recently begun to attempt to use tespy for the modelling of heat recovery solutions. However, as an initial step, I wanted to create a simple model of a dryer.
For this, I based myself off of the SimpleHeatExchanger component and attempted to create a new, custom component. I have attached the code for this custom component:

This custom component would allow input of an amount of water to be "added" to the inlet, which will correspond to the water evaporated from the material we are drying. However, any attempt to manipulate the composition of the fluid has ended in failure. I have been able to manipulate the enthalpy calculations to obtain an acceptable shell from which to work around (though it is still a work in progres). I am however, stumped as to how the composition could be altered. I would be very grateful if anyone could help.
Respectfully,
Manuel.

from tespy.components.component import Component
from tespy.tools import ComponentCharacteristics as dc_cc
from tespy.tools import ComponentProperties as dc_cp
from tespy.tools.document_models import generate_latex_eq

      class dryer(Component):
          
        This is a custom component.
    
        You can add your documentation here. From this part, it should be clear
        for the user, which parameters are available, which mandatory equations
        are applied and which optional equations can be applied using the
        component parameters.
        
    
        def component(self):
            return 'dryer'
        
        
        
        def get_parameters(self):
            return {
                'Q': dc_cp(
                    min_val=-1e10, max_val = 1e10, num_eq=1,
                    deriv =self.Q_deriv, func = self.Q_func,
                    func_params={'Q':'Q'}, latex=self.Q_func_doc),
                'm': dc_cp(
                    min_val=-1e10, max_val = 1e10, num_eq=1,
                    deriv =self.m_deriv, func = self.m_func,
                    func_params={'m':'m'}, latex = self.m_func_doc),
                'Losses': dc_cp(
                    min_val = 0, max_val = 1e10, num_eq=1,
                    deriv = self.loss_func_deriv, func = self.loss_func,
                    func_params={'loss':'loss'}, latex = self.loss_func_doc),
                'evap': dc_cp(
                    min_val = 0, max_val = 1e3, d=1e-4,
                    ),
                    
                'Qevap': dc_cp(
                    min_val = 0, max_val = 1e10, num_eq=1,
                    deriv = self.qevap_func_deriv, func = self.qevap_func,
                    func_params={'qevap':'qevap'}, 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,
                     func_params={'qheatin':'qheatin'}, 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 inlets(self):
            return{'in1'}
        
        
        
        def outlets(self):
            return{'out1'}
        
        
        
    
        
        
        
        def Q_func(self):
            r"""
            Equation for dryer heat load.
    
            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.Q.val
            )
        def Q_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 Q_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 m_func(self):
            r"""
            Equation for dryer mass expansion.
    
            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].m_val_SI - self.evap.val
                
            )
        def m_deriv(self,increment_filter,k):
            r"""
            Partial derivatives for dryer max expansion.
    
            Parameters
            ----------
            increment_filter : ndarray
                Matrix for filtering non-changing variables.
    
            k : int
                Position of derivatives in Jacobian matrix (k-th equation).
            """
            return()
        
        def m_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 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.loss.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 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)
            """
            return (
                (self.mevap.val*2273*1e3)-self.qevap.val
            )
        def qevap_func_deriv(self,increment_filter,k):
            r"""
            Partial derivatives for water evaporation.
    
            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_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"""
            calculating heat required for material .
    
            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"""
            Partial derivatives for material 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 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 calc_parameters(self):
            self.Qevap.val = 2273* self.mevap.val
            self.Qheatin.val = self.Mmat.val * self.cp.val * self.Dt.val 
            
            self.Q.val=self.inl[0].m.val_SI * (
                self.outl[0].h.val_SI - self.inl[0].h.val_SI) 
            self.outl[0].fluid.val['H2O']= self.inl[0].fluid.val['H2O'] + self.evap.val
            
            self.Losses.val= self.Q.val - self.Qevap.val - self.Qheatin.val 

[MMArro]

Upon reinspection, I have much refactoring to do. I will create a new issue when I have a clearer idea in mind for what exactly is going wrong.