No title |POLVER05_1 d(T)/d(t) = -1000 # d(nm)/d(t) = -rN*nstar - pi*dp^2*N*rC # d(N)/d(t) = rN - 0.5*rF # d(V)/d(t) = rN*nstar*v1 + pi*(dp^2)*v1*N*rC # d(x)/d(t) = 4*Q/(pi*d^2) # P = 1.013e5 # Mw = 39.948/1000 # R = 8.314 # pi = 3.141 # rhop = 2.7*1000 # kB = R/6.022E23 # m1 = 4.48E-26 # v1 = 1.23E-29 # s1 = 2.575E-19 # Qrt = 2/(1000*60) # Trt = 273 # d = 18/1000 # gama = 4.6325e-19 # Q = Qrt*T/Trt # rhog = P*Mw/(R*T) # meu = (26.69*sqrt(Mw*1000*T)/((dg*1.0e+10)^2))*1.0E-7 # ps = exp(13.07 - 36373/T)*P # lamda = 1/(sqrt(2)*gama)*kB*T/P # sigma = (948 - 0.202*T)/1000 # nms = ps/(kB*T) # theta = sigma*s1/(kB*T) # S = nm*rhog/nms # rN = if (S>1) then ((v1/rhog)*sqrt(2*sigma/(pi*m1))*(nms^2)*S*exp(theta-(4*theta^3/(27*ln(S)^2)))) else (0) # nstar = if (S>1) then ((2*theta/(3*ln(S)))^3) else (0) # dp = if ((V<0) or (N<0)) then ((6*v1/pi)^(1/3)) else ((6*V/((N+1)*pi))^(1/3)+(6*v1/pi)^(1/3)) # Kn = 2*lamda/dp # fKn = (1.333*Kn+1.333*Kn^2)/(1+1.71*Kn+1.333*Kn^2) # vp = pi*dp^3/6 # c = sqrt((8*kB*T)/(pi*rhop*vp)) # D = kB*T/(3*pi*meu*dp)*(5+4*Kn+6*Kn^2+18*Kn^3)/(5-Kn+(8+pi)*Kn^2) # Ia = 8*D/(pi*c) # g1 = 1/(3*dp*Ia)*((dp+Ia)^3-(dp^2+Ia^2)^(3/2))-dp # beta = 8*pi*D*dp*((dp/(dp+sqrt(2)*g1))+(4*sqrt(2)*D/(c*dp)))^(-1) # dg = 0.384e-9 # rC = if (S<1) then (0) else (nms*(S-1)*sqrt(kB*T/(2*pi*m1))*fKn) # rF = beta*N^2*rhog # t(0)=0 T(0)=1973 nm(0)=2.634e+022 N(0)=0 V(0)=0 x(0)=0 t(f)=1