$SB1-X5.226
Function NTU_epsilon(C_c,C_f,epsilon,tipo,n)
If (n=0) Then CALL ERROR ('El valor de n no puede ser cero, toma el valor del número de pasos por la carcasa, cuando el valor de la variable tipo=4, y toma el valor 1 en cualquier otro caso')
If (tipo<>1) and (tipo<>2) and (tipo<>3) and (tipo<>4) and (tipo<>5) and (tipo<>6) Then CALL ERROR ('El valor de la variable tipo no es correcto, solo puede tomar valores enteros desde 1 hasta 6, ambos incluidos')
Cmin=min(C_c,C_f)
Cmax=max(C_c,C_f)
Cr=(Cmin/Cmax)
If (Cr=0) Then
NTU_epsilon=-ln(1-epsilon)
EndIf
If (tipo=1) Then
NTU_epsilon=-((ln(1-(epsilon*(1+Cr))))/(1+Cr))
EndIf
If (tipo=2) and (Cr<1) Then
NTU_epsilon=(1/(Cr-1))*ln((epsilon-1)/((epsilon*Cr)-1))
EndIf
If (tipo=2) and (Cr=1) Then
NTU_epsilon=epsilon/(1-epsilon)
EndIf
If (tipo=3) and (Cr>0) Then
E=((2/epsilon)-(1+Cr))/((1+Cr^2)^(1/2))
NTU_epsilon=-(((1+Cr^2)^(-1/2))*ln((E-1)/(E+1)))
EndIf
If (tipo=4) and (Cr<1) Then
F=(((epsilon*Cr)-1)/(epsilon-1))^(1/n)
epsilon_1=(F-1)/(F-Cr)
E=((2/epsilon_1)-(1+Cr))/((1+Cr^2)^(1/2))
NTU=-(((1+Cr^2)^(-1/2))*ln((E-1)/(E+1)))
NTU_epsilon=n*NTU
EndIf
If (tipo=4) and (Cr=1) Then
Cr_2=Cr-0.001
F=(((epsilon*Cr_2)-1)/(epsilon-1))^(1/n)
epsilon_1=(F-1)/(F-Cr_2)
E=((2/epsilon_1)-(1+Cr_2))/((1+Cr_2^2)^(1/2))
NTU=-(((1+Cr_2^2)^(-1/2))*ln((E-1)/(E+1)))
NTU_epsilon=n*NTU
EndIf
If (tipo=5) and (Cr>0) Then
NTU_epsilon=-(ln(1+((1/Cr)*ln(1-(epsilon*Cr)))))
EndIf
If (tipo=6) and (Cr>0) Then
NTU_epsilon=-((1/Cr)*ln((Cr*ln(1-epsilon))+1))
EndIf
End
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
