Codice:
Program hwup
Implicit Real*8(a-z)
Parameter (bolz=1.381d-23, pi=3.14159, c=1d0, prh2o=5.8, preg=150, a=4d4, m=2.5)
Integer*2 x
print *,"Desideri usare acqua o glicoletilenico?"; print*, "1)Acqua"; print*,"2)EG"
read*, x
print *, "Calcolo il valore con Hamilton-Crosser"
print*, "Inserisci il valore della conduttivita' termica delle nanoparticelle"
read*, kp
print*, "Inserisci la sfericita' delle nanoparticelle"
read*, n
print*, "Inserisci la frazione volumetrica delle nanoparticelle"
read*, f
if (x==1) then; ham = hamiltoncrosser(kp,0.613d0,n,f); else; ham = hamiltoncrosser(kp,0.253d0,n,f); endif
print*, "Valore calcolato con la formula di Hamilton-Crosser:", ham
print*, "Calcolo il valore con Kumar"
print*, "Inserisci il valore della temperatura"
read*, t
print*, "Inserisci il valore del diametro medio nanoparticellare"
read*, dp
print*, "Inserisci il valore della viscosita' dinamica del fluido"
read*, v
if (x==1) then
kum=kumar(t,dp,f,0d0,0.613d0,v)
else
kum=kumar(t,dp,f,0d0,0.253d0,v)
end if
print*, "Valore calcolato con la formula di Kumar:", kum
print*, "Calcolo il valore con Prasher"
print*, "Inserisci il valore della densita' della nanoparticella"
read*, p
if(x==1) then
rey=reynolds(0.613d0,t,p,dp,v)
pra=prasher(kp,0.613d0,f,prh2o,rey)
else
rey=reynolds(0.253d0,t,p,dp,v);
pra=prasher(kp,0.253d0,f,preg,rey);
endif
print*, "Valore calcolato con la formula di Prasher:", pra
CONTAINS
Function hamiltoncrosser(kp, kb, n, f)
hamiltoncrosser=((kp+(((3d0/n)-1d0)*kb)-(((3d0/n)-1)*(kb-kp)*f))/(kp+(((3d0/n)-1)*kb)+((kb-kp)*f)));
End Function
Function reynolds(kb, t, p, dp, v)
reynolds=sqrt((18d0*(kb*t))/(pi*p*dp))/v
End Function
Function kumar(t, dp, f, db, kb, v)
kumar=(((((2*bolz*t)/(pi*v*dp*dp)))*((f*db/2)/(kb*(1-f)*(dp/2))))*c)+1;
End Function
Function prasher(kp, kb, f, pr, rey)
prasher=((kp+(2*kb)+2*(kp-kb)*f)/(kp+(2*kb)-(kp-kb)*f))*(a*rey**m*pr**0.333*f)
End Function
End Program hwup