File: N:\mfix\model\GhdTheory\pressure.f

1     !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv
2     !
3     !  subroutine name: pressure(s,alpha,ni,n,mu,sigma,chi,T,Ti,p)
4     !
5     !  author:  C. Hrenya, Jan 2009
6     !
7     !  Purpose: find pressure 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 pressure(s,alpha,ni,n,mu,sigma,chi,T,Ti,p)
15           Implicit NONE
16     
17           integer s
18     
19           double precision alpha(s,s),ni(s),n,mu(s,s),sigma(s,s), &
20                            chi(s,s),T,Ti(s),p
21     
22           integer i,j
23           double precision pkin, pcol
24     
25           double precision pi
26           parameter (pi=3.14159265458979323846d0)
27     
28           pkin = n*T                   !kinetic pressure
29           pcol = 0.d0
30           do i=1,s
31              do j=1,s
32                pcol = pcol + mu(j,i)*(1.d0+alpha(i,j))*sigma(i,j)**3* &
33                                 chi(i,j)*ni(i)*ni(j)*Ti(i)
34              enddo
35           enddo
36           pcol = 2.d0*pi/3.d0*pcol     !collisional pressure
37           p = pkin + pcol              !total pressure (p 9 CMH notes)
38     
39           return
40           end subroutine pressure
41     
42