{$DS.}Function EFICIENCIA(C_c,C_f,NTU,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) and (tipo<>7) Then CALL ERROR ('el valor de la variable tipo no es correcto, solo puede tomar valores enteros desde 1 hasta 7, ambos incluidos')
Cmin=min(C_c,C_f)
Cmax=max(C_c,C_f)
Cr=(Cmin/Cmax)
If (Cr=0) Then
EFICIENCIA=1-exp(-NTU)
EndIf
If (tipo=1) Then
EFICIENCIA=(1-exp(-NTU*(1+Cr)))/(1+Cr)
EndIf
If (tipo=2) and (Cr<1) Then
EFICIENCIA=(1-exp(-NTU*(1-Cr)))/(1-(Cr*exp(-NTU*(1-Cr))))
EndIf
If (tipo=2) and (Cr=1) Then
EFICIENCIA=NTU/(1+NTU)
EndIf
If (tipo=3) Then
A=1+exp(-NTU*((1+Cr^2)^(1/2)))
B=1-exp(-NTU*((1+Cr^2)^(1/2)))
EFICIENCIA=2*(1/(1+Cr+(((1+Cr^2)^(1/2))*(A/B))))
EndIf
If (tipo=4) and (Cr<1) Then
A=1+exp(-(NTU/n)*((1+Cr^2)^(1/2)))
B=1-exp(-(NTU/n)*((1+Cr^2)^(1/2)))
C=2*(1/(1+Cr+(((1+Cr^2)^(1/2))*(A/B))))
X=(((1-(C*Cr))/(1-C))^n)-1
Y=(((1-(C*Cr))/(1-C))^n)-Cr
EFICIENCIA=X*(Y^(-1))
EndIf
If (tipo=4) and (Cr=1) Then
Cr_2=Cr-0.0001
A=1+exp(-(NTU/n)*((1+Cr_2^2)^(1/2)))
B=1-exp(-(NTU/n)*((1+Cr_2^2)^(1/2)))
C=2*(1/(1+Cr_2+(((1+Cr_2^2)^(1/2))*(A/B))))
X=(((1-(C*Cr_2))/(1-C))^n)-1
Y=(((1-(C*Cr_2))/(1-C))^n)-Cr_2
EFICIENCIA=X*(Y^(-1))
EndIf
If (tipo=5) and (Cr>0) Then
A1=((Cr)^(-1))*NTU^(0.22)
A2=exp(-Cr*(NTU^0.78))-1
EFICIENCIA=1-exp(A1*A2)
EndIf
If (tipo=6) and (Cr>0) Then
EFICIENCIA=(1/Cr)*(1-exp(-Cr*(1-exp(-NTU))))
EndIf
If (tipo=7) and (Cr>0) Then
EFICIENCIA=1-exp(-(1/Cr)*(1-exp(-Cr*NTU)))
EndIf
End
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
