TITLE      'Profili termici in un condotto, calcolo dei valori medi '         { CondottiB }
SELECT

COORDINATES
 ycylinder('r','z')

VARIABLES
ITb	Por

DEFINITIONS
transfermesh('transfer.dat',T,vz,vr,p,dTsudr,dvsudr)  {Importa i dati da "CondottiA"}

r1=5.0e-3
L=5*r1
vz0=1e-1
visc=2.5e-5
dens=0.785	
cond=0.0363	
Cp=1020	
rcp=dens*Cp
T0=eval(T,0,0)	{ Temperatura di ingresso }
Tw=eval(T,r1,z)  { Temperatura di parete }

Tb=ITb/Por
vav=Por/(PI*r1^2)

Re=dens* vz0*2*r1/ visc   { Numero di Reynolds }
Pr=Cp*visc/cond	             { Numero di Prandtl }

Lcin=2*r1*0.05*Re        { Lunghezza di imbocco cinematica}
Lter=2*r1*0.05*Re*Pr     { Lunghezza di imbocco termica}

q=Eval(-cond*dTsudr,r1,z)	 	tau=Eval(-visc*dvsudr,r1,z)
h=abs(q/(Tb-Tw))			f=tau/(1/2*dens*vz0^2)
hfinal=eval(h,r1,L)		ffinal=eval(f,r1,L)	
GrL=Re*Pr*2*r1/L
NuL=hfinal*2*r1/cond


EQUATIONS
ITb:	dr(ITb)-2*PI*r*vz*T= 1e-16*dzz(ITb)
Por:	 dr(Por)-2*PI*r*vz= 1e-16*dzz(Por)

BOUNDARIES
region 'domain' start 'outer' (0,0)
line to (r1,0) { In }
natural(ITb)=0
natural(Por)=0	
line to (r1,L) { Wall }
natural(ITb)=0
natural(Por)=0	
 line to (0,L) { Out }
value(ITb)=0  	
value(Por)=0
 line to close { Axis }

MONITORS
elevation(Tb) as "Temperatura di bulk" from(r1,-L) to (r1,L)
elevation(vav,vz0) as "velocità media" from(r1,-L) to (r1,L)
elevation(q) as "Flusso termico" from(r1,-L) to (r1,L)
elevation(f) as "Fattore di attrito" from(r1,-L) to (r1,L)
elevation(h) as "Coefficiente di scambio termico"  from(r1,-L) to (r1,L)


PLOTS
elevation(Tb) as "Temperatura di bulk"  from(r1,-L) to (r1,L) export format "#y#b#1" file="Tbtable.txt"
elevation(q) as "velocità media" from(r1,-L) to (r1,L)  export format "#y#b#1" file="qtable.txt"
elevation(f) as "Fattore di attrito" from(r1,-L) to (r1,L) export format "#y#b#1" file="ftable.txt"
elevation(h) as "Coefficiente di scambio termico" from(r1,-L) to (r1,L) export format "#y#b#1" file="htable.txt"
surface(T)  export format "#x#y#b#1" file="T.txt"


summary
report(Re)
report(Pr)
report(GrL)
report(NuL)
report(hfinal)
report(ffinal)
report(Lcin)
report(Lter)

END