Dim ß As Double, DB As Double, DR As Double, Ep As Double, P1 As Double, P2 As Double
Dim E As Double, Rho1 As Double, PW As Double, P1abs As Double, Pamb As Double, T1 As Double
Dim my As Double, Tr As Double, c As Double, my0 As Double, Var1 As Double, Var2 As Double
Dim Calt As Double, Re As Double, A As Double, qm As Double, qv As Double, sg As Double
Dim C_neu As Double, Ai As Double, Re_i As Double, pi As Double, Var3 As Double, Var4 As Double
Dim counter As Integer
DB = TextBox2.Value / 1000
DR = TextBox1.Value / 1000
P1 = TextBox4.Value * 10 ^ 5
P2 = TextBox5.Value * 10 ^ 5
Pamb = TextBox7.Value * 10 ^ 5
PW = TextBox6.Value * (10 ^ 5)
T1 = TextBox3.Value + 273.15
P1abs = P1 + Pamb
Rho1 = P1abs / ((287 * T1))
ß = DB / DR
E = 1 / Sqr(1 - ß ^ 4)
Var1 = ß ^ 4
Var2 = ß ^ 8
Ep = 1 - (0.351 + 0.256 * Var1 + (0.93 * Var2)) * (1 - ((P2 / P1) ^ (1 / 1.4)))
Re = 10 ^ 5
counter = 0
pi = 3.14159265358979
Var3 = ß ^ 2
Var4 = ß ^ 3.5
Do
Calt = 0.5961 + 0.0261 * Var3 - 0.216 * Var2 + 0.000521 * ((10 ^ 6 * ß) / Re) ^ 0.7 + (0.0188 + 0.0063 * A) * (ß ^ 3.5) * (10 ^ 6 / Re) ^ 0.3
A = (19000 * ß / Re) ^ 0.8
qm = Ep * Calt * (pi / 4 * DB ^ 2) * E * ((2 * PW * Rho1) ^ 0.5) * 3600
qv = qm / Rho1
sg = qv / (3.14 / 4 * DR ^ 2) / 3600
Tr = 291.15
my0 = 18.27 * (10 ^ -6)
c = 120
my = my0 * (Tr + c) / (T1 + c) * ((T1 / Tr) ^ (3 / 2))
Re_i = sg * DR * Rho1 / my
C_neu = 0.5961 + 0.0261 * (ß ^ 2) - 0.216 * (ß ^ 8) + 0.000521 * (10 ^ 6 * ß / Re_i) ^ 0.7 + (0.0188 + 0.0063 * A) * (ß ^ 3.5) * (10 ^ 6 / Re_i) ^ 0.3
A = (19000 * ß / Re_i) ^ 0.8
counter = counter + 1
If Abs(C_neu - Calt) / Calt < 10 ^ -6 Then
Exit Do
Else
Re = Re_i
End If
Loop