1 |
gezelter |
3119 |
program nanodecay |
2 |
|
|
|
3 |
|
|
implicit none |
4 |
|
|
real(kind=8) :: k, R, g, T0, pi, kappa, u, t, M |
5 |
|
|
real(kind=8) :: BigK, BigG, rho, c, cp, umin, umax, u2, du |
6 |
|
|
real(kind=8) :: rhop, sum, numer, denom, Tf0 |
7 |
|
|
integer :: i, j, kk, nu |
8 |
|
|
|
9 |
|
|
! fixed quantities: |
10 |
|
|
pi = 4.0d0 * datan(1.0d0) |
11 |
|
|
! initial particle temperature |
12 |
|
|
T0 = 950.0 |
13 |
|
|
! fluid phase temperature |
14 |
|
|
Tf0 = 300.0 |
15 |
|
|
! |
16 |
|
|
! specific heat of particle (J g^-1 K^-1) |
17 |
|
|
cp = 0.4*0.385 + 0.6*0.235 |
18 |
|
|
! density of particle (g m^-3) |
19 |
|
|
rhop = (0.4*8920 + 0.6 * 10490) * 1000.0 |
20 |
|
|
! thermal conductivity of fluid (W m^-1 K^-1) |
21 |
|
|
BigK = 0.6 |
22 |
|
|
! density of fluid |
23 |
|
|
rho = 1.0e6 |
24 |
|
|
! specific heat of fluid (J g^-1 K^-1) |
25 |
|
|
c = 4.184 |
26 |
|
|
! interfacial conductance (W m^-2 K^-1) |
27 |
|
|
BigG = 105.0 * 1e6 |
28 |
|
|
|
29 |
|
|
! derived quantities |
30 |
|
|
! |
31 |
|
|
kappa = BigK / (rho * c) |
32 |
|
|
g = BigG / BigK |
33 |
|
|
|
34 |
|
|
umin = 0.001 |
35 |
|
|
umax = 1000.0 |
36 |
|
|
|
37 |
|
|
nu = 10000 |
38 |
|
|
du = (umax-umin)/real(nu) |
39 |
|
|
|
40 |
|
|
do i = 20, 50, 10 |
41 |
|
|
! radius of particle (in m) |
42 |
|
|
R = real(i)*1e-10 |
43 |
|
|
! mass of particle |
44 |
|
|
M = 4.0*pi*R*R*R*rhop/3.0 |
45 |
|
|
|
46 |
|
|
k = 4.0 * pi * R*R*R * rho * c / (M * cp) |
47 |
|
|
|
48 |
|
|
do j = 0, 1000000, 100 |
49 |
|
|
! time in s |
50 |
|
|
t = real(j)*1e-15 |
51 |
|
|
|
52 |
|
|
sum = 0.0 |
53 |
|
|
|
54 |
|
|
do kk = 0, nu |
55 |
|
|
|
56 |
|
|
u = umin + (umax-umin)*real(kk)/nu |
57 |
|
|
u2 = u*u |
58 |
|
|
|
59 |
|
|
numer = exp(-kappa * u2 * t / (R*R)) * u2 |
60 |
|
|
denom = (u2*(1.0+ R*g) - k*R*g)**2 + (u2*u - k*R*g*u)**2 |
61 |
|
|
|
62 |
|
|
sum = sum + du*numer/denom |
63 |
|
|
|
64 |
|
|
enddo |
65 |
|
|
sum = sum*2.0*k*R*R*g*g*(T0-Tf0)/pi |
66 |
|
|
|
67 |
|
|
write(*,*) t*1e15, Tf0 + sum |
68 |
|
|
enddo |
69 |
|
|
write(*,*) '&' |
70 |
|
|
enddo |
71 |
|
|
end program nanodecay |
72 |
|
|
|
73 |
|
|
|
74 |
|
|
|
75 |
|
|
|
76 |
|
|
|
77 |
|
|
|
78 |
|
|
|
79 |
|
|
|