Sub Tsat()
Dim A0, A1, A2, A3, A4, A5, A6, A7 As Double ' Double p. cifras significativas
Dim P, Tsat As Double ' Double p. cifras significativas
'Parametros para el agua
A0 = -20.69171
A1 = 13.610819
A2 = 7.50043
A3 = -0.37146631
A4 = 0.26070017
A5 = -0.0240092
A6 = 0.0013378831
A7 = -8.9845187E-33
P = Cells(2, "B")
X = Log(A2 * P)
Tsat = A0 + A1 * X + A3 * X ^ 2 + A4 * X ^ 3 + A5 * X ^ 4 + A6 * X ^ 5 + A7 * X ^ 30
Cells(3, "B") = Tsat
End Sub
Sub Psat()
Dim T, Psat As Double
Dim Tabs, B0, B1, B2, B3, B4 As Double
'Parametros para el agua
B0 = 98.0904
B1 = 21.675595
B2 = -0.017091984
B3 = -6249.2
B4 = 0.000010643944
T = Cells(5, "B")
Tabs = T + 273.15 'K
Psat = B0 * Exp(B1 + B2 * Tabs + B3 / Tabs + B4 * (Tabs ^ 2))
Cells(6, "B") = Psat
End Sub
Sub EntalpiaLiquidoSaturado()
Dim Psat1, Hlsat As Double
Dim Pkgf, C0, C1, C2, C3, C4, C5, C6, C7, C8, C9 As Double
'Parametros para el agua
C0 = -38196.597
C1 = 8156.7769
C2 = 0.062913098
C3 = 5.9407557E-13
C4 = 1.1973266E-59
C5 = 135722.09
C6 = -170161.83
C7 = 66738.661
C8 = -2246.2754
C9 = 402.37285
Psat1 = Cells(8, "b") 'Kpa
Pkgf = Psat1 / 98.0904
Hlsat = C0 + (C1 * Pkgf ^ 0.1) + (C2 * Pkgf ^ 1.5) + (C3 * Pkgf ^ 6) + (C4 * Pkgf ^ 26) + (C5 / Pkgf ^
0.1) + (C6 / Pkgf ^ 0.15) + (C7 / Pkgf ^ 0.2) + (C8 / Pkgf ^ 0.4) + (C9 / Pkgf ^ 0.5)
Cells(9, "B") = Hlsat
End Sub
Sub EntalpiaVaporSaturado()
Dim Psat2, Hvsat As Double
Dim Pkgf1, D0, D1, D2, D3, D4, D5, D6, D7, D8, D9 As Double
'Parametros para el agua
D0 = 16394.097
D1 = -2326.9083
D2 = -0.12038685
D3 = -1.11101E-12
D4 = -1.2893239E-59
D5 = -46574.48
D6 = 54360.668
D7 = -19508.785
D8 = 368.64023
D9 = -38.23317
Psat2 = Cells(11, "B")
Pkgf1 = Psat2 / 98.0904
Hvsat = D0 + D1 * Pkgf1 ^ 0.1 + D2 * Pkgf1 ^ 1.5 + D3 * Pkgf1 ^ 6 + D4 * Pkgf1 ^ 26 + (D5 / Pkgf1 ^
0.1) + (D6 / Pkgf1 ^ 0.15) + (D7 / Pkgf1 ^ 0.2) + (D8 / Pkgf1 ^ 0.4) + (D9 / Pkgf1 ^ 0.5)
Cells(12, "B") = Hvsat
End Sub
Sub EntalpiaVaporSobrecalentado()
Dim Tvsc, Tsat2, P2, Hvsat2, Hvsc, Fact, DT, VSC(10, 6) As Double
Dim E(1 To 6) As Variant
Dim i, j As Integer
Tvsc = Cells(14, "B") '°C
Tsat2 = Cells(15, "B") '°C
P2 = Cells(16, "B") 'kPa
Hvsat2 = Cells(17, "B") 'kJ/kg
For i = 1 To 10
For j = 1 To 6
VSC(i, j) = Sheets("Base de datos").Cells(1 + i, 1 + j)
Next
Next
Pkg1 = P2 / 98.0904 'Unidades inglesas
For j = 1 To 6
E(j) = VSC(1, j) + VSC(2, j) * Pkg1 + VSC(3, j) * (P ^ 2) + VSC(4, j) * (Pkg1 ^ 3) + VSC(5, j) * (Pkg1 ^ 4) +
VSC(6, j) * (Pkg1 ^ 5) + VSC(7, j) * (Pkg1 ^ 6) + VSC(8, j) / (Pkg1) + VSC(9, j) / (Pkg1 ^ 2) + VSC(10, j) *
(Pkg1 ^ 0.5)
Next
DT = Tvsc - Tsat2 '°C
Fact = E(1) + E(2) * DT + E(3) / DT + E(4) * DT ^ 2 + E(5) / DT ^ 2 + E(6) * DT ^ 0.5
Hvsc = Hvsat2 + Fact * DT 'kJ/kg
Cells(18, "B") = Hvsc ' kJ/kg
End Sub
Sub EntropiaSat()
Dim S(1 To 2) As Variant
Dim P3, Psia, BD(7, 2) As Double
Dim i, j As Integer
P3 = Cells(2, "G") ' kPa
For i = 1 To 7
For j = 1 To 2
BD(i, j) = Sheets("Base de datos").Cells(1 + i, 11 + j)
Next
Next
Psia = P3 * 0.14504 ' Psia
For j = 1 To 2
S(j) = BD(1, j) * Psia + BD(2, j) / Psia + BD(3, j) * (Psia ^ (1 / 2)) + BD(4, j) * Log(Psia) + BD(5, j) * (Psia
^ 2) + BD(6, j) * (Psia ^ 3) + BD(7, j)
Next
Cells(3, "G") = S(1)
Cells(4, "G") = S(2)
End Sub
Sub EntalpíaVap()
Dim Tb, Tc, T1 As Double
Dim DHvap1, DHvap2 As Double
Tb = Cells(9, "G") + 273.15 'Conversión a Kelvin
Tc = Cells(10, "G") + 273.15
T2 = Cells(11, "G") + 273.15
DHvap1 = 0.109 * Tb
DHvap2 = DHvap1 * ((Tc - T2) / (Tc - Tb)) ^ 0.38
Cells(12, "G") = DHvap1
Cells(13, "G") = DHvap2
End Sub
Sub EntropíaVaporización()
Dim DHvap3, Tb1, DSvap As Double
Tb1 = Cells(15, "G")
DHvap3 = Cells(16, "G")
DSvap3 = DHvap3 / Tb1
Cells(17, "G") = DSvap3
End Sub
Sub EntropiaVaporSobrecalentado()
Dim DSvap4, Hvsc, Hvsat, Tvsc, Tsat, Tprom, Ssc As Double
DSvap4 = Cells(2, "L")
Hvsc = Cells(3, "L")
Hvsat = Cells(4, "L")
Tvsc = Cells(5, "L")
Tsat = Cells(6, "L")
Tprom = (Tvsc * Tsat) ^ 1 / 2
Ssc = DSvap4 + (Hvsc - Hvsat) / (Tprom)
Cells(7, "L") = Ssc
End Sub
Sub ViscosidadLiquidos()
Dim T, u As Double
T = Cells(9, "L")
u = Exp(-3.63148 + ((542.05) / (T + 129)))
Cells(10, "L") = u
End Sub
Sub DensidadLiquidos()
Dim Tf, D, Drel As Double
Tf = Cells(12, "L") * (1.8) + 32 ' Celsis a Kelvin
Drel = 0.997375 + 0.0001201 * Tf - 0.000001601 * Tf ^ 2 + 0.000000001601 * Tf ^ 3
D = 62.47 * Drel ' lb/ft3
Cells(13, "L") = D
End Sub
Sub TensiónSuperficialLíquidos()
Dim Tf1, Tsup As Double
Tf1 = Cells(15, "L") * (1.8) + 32 ' Celsius a Kelvin
Tsup = 79.5118 - 0.09605 * Tf1 ' dina/cm
Cells(16, "L") = Tf1
End Sub
Sub ConductividadTérmica()
Dim Tf2, K As Double
Tf2 = Cells(18, "L") * (1.8) + 32 ' Celsius a Kelvin
K = 0.31431 + 0.00047673 * Tf2
Cells(19, "L") = K
End Sub
