>	restart:#"m12_p09"

Determinar la velocidad de evaporación de una gota de 1 mm de diámetro moviéndose a 10 m/s en aire a 25 ºC, 50%RH y 100 kPa.

Datos:

>	read`../therm_eq.m`:read`../therm_const.m`:read`../therm_proc.m`:with(therm_proc):

>	su1:="Aire":su2:="H2O":dat:=[D=0.001*m_,v=10*m_/s_,T0=(25+273)*K_,phi0=0.5,Di=22e-6*m_^2/s_,Sh=2+0.6*Re_^(1/2)*Sc^(1/3),Nu=2+0.6*Re_^(1/2)*Pr^(1/3)];

Esquema:

Eqs. const.:

>	Adat:=get_gas_data(su1):Adat:=subs(c[p]=c[pa],R=R[a],M=M[a],T[b]=nada,[Adat]):Wgdat:=get_gas_data(su2):Wgdat:=subs(c[p]=c[pv],R=R[v],M=M[v],[Wgdat]):Wldat:=get_liq_data(su2):Wdat:=op(Wgdat),Wldat:get_pv_data(su2):dat:=op(dat),Const,SI2,SI1:

>

mu_:=subs(Adat,mu);rho_:=subs(Adat,dat,p0/(R[a]*T0));nu_:=subs(Adat,dat,mu/rho_);Sc_:=subs(dat,nu_/Di);Re_:=subs(dat,v*D/nu_);Sh__:=subs(dat,Sc=Sc_,dat,Sh);hm:=Sh*Di/D;hm_:=subs(dat,Sh__*Di/D);eqTwet:=h1=hsa;h=h(T,w);eq8_8;w0_:=subs(dat,evalf(subs(dat,w(phi0,T0,p0))));h0_:=subs(Adat,Wdat,T=T0,dat,h(T,w0_));Tsa_:=fsolve(subs(dat,SI0,h0_=subs(Adat,Wdat,h(T,w(1,T,p0)))),T=200..400)*K_;Tsa=TKC(%);wsa_:=subs(dat,evalf(subs(dat,w(1,Tsa_,p0))));A_:=evalf(subs(dat,Pi*D^2));mdot_:=subs(dat,hm_*A_*rho_*(wsa_-w0_));m0_:=evalf(subs(Wdat,dat,rho*Pi*D^3/6));tevap_:=subs(dat,m0_/mdot_);

>	Pr_:=subs(Adat,dat,mu_*c[pa]/k);Nu__:=subs(dat,Pr=Pr_,dat,Nu);h_:=subs(Adat,dat,Nu__*k/D);Qdot_:=subs(dat,h_*A_*(T0-Tsa_));mdot_:=subs(Wdat,dat,Qdot_/h[lv0]);tevap_:=subs(dat,m0_/mdot_);

>