File: RELATIVE:/../../../mfix.git/model/GhdTheory/bulk_viscosity.f

1     !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv
2     !
3     !  subroutine name: bulk_viscosity(s,mi,alpha,ni,v0,mu,sigma,chi,theta,kappa)
4     !
5     !  author:  C. Hrenya, Jan 2009
6     !
7     !  Purpose: find bulk viscosity according to GHD polydisperse KT
8     !
9     !  Literature/References:
10     !     C. Hrenya handwritten notes & Garzo, Hrenya, Dufty papers (PRE, 2007)
11     !
12     !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
13     
14           subroutine bulk_viscosity(s,mi,alpha,ni,v0,mu,sigma,chi, &
15                                     theta,kappa)
16           Implicit NONE
17     
18           integer s
19     
20           double precision mi(s),alpha(s,s),ni(s),v0,mu(s,s),sigma(s,s), &
21                            chi(s,s),theta(s),kappa
22     
23           integer i,j
24     
25           double precision pi
26           parameter (pi=3.14159265458979323846d0)
27     
28           kappa = 0.d0
29           do i=1,s
30              do j=1,s
31                 kappa = kappa + mu(i,j)*mi(j)*ni(i)*ni(j)*v0* &
32                           sigma(i,j)**4*chi(i,j)*(1.d0+alpha(i,j)) &
33                       *dsqrt((theta(i)+theta(j))/(theta(i)*theta(j)))
34              enddo
35           enddo
36           kappa = 2.d0*dsqrt(pi)/9.d0*kappa        !p 12 CMH notes
37     
38           return
39           end subroutine bulk_viscosity
40     
41