b_m_p_star.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: B_m_P_star_e(IJK, IER)                                 C
!  Purpose: Determine p_star when there is an interface at east        C
!                                                                      C
!                                                                      C
!  Author: M. Syamlal                                 Date: 6-AUG-96   C
!  Reviewer:                                          Date:            C
!                                                                      C
!                                                                      C
!  Literature/Document References:                                     C
!                                                                      C
!  Variables referenced:                                               C
!  Variables modified:                                                 C
!                                                                      C
!  Local variables:                                                    C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
!
      SUBROUTINE b_m_p_star_e(B_M, IJK)
!...Translated by Pacific-Sierra Research VAST-90 2.06G5  12:17:31  12/09/98
!...Switches: -xf
!
!  Include param.inc file to specify parameter values
!
!-----------------------------------------------
!   M o d u l e s
!-----------------------------------------------
      USE param
      USE param1
      USE scales
      USE constant
      USE physprop
      USE fldvar
      USE run
      USE rxns
      USE toleranc
      USE geometry
      USE indices
      USE compar
      USE bodyforce
      USE fun_avg
      USE functions
      IMPLICIT NONE
!-----------------------------------------------
!   G l o b a l   P a r a m e t e r s
!-----------------------------------------------
!-----------------------------------------------
!   D u m m y   A r g u m e n t s
!-----------------------------------------------
!
!                      Indices
      INTEGER          I, J, K, IJK, IJKE
!
!                      Phase index
      INTEGER          M
!
!                      Average volume fraction
      DOUBLE PRECISION EPGA
!
!                      Average density
      DOUBLE PRECISION ROPGA
!
!                      LHS (similar to A_m)
      DOUBLE PRECISION A
!
!                      RHS (similar to b_m)
      DOUBLE PRECISION BB
!
!                      b_m
      DOUBLE PRECISION b_m
!
!                      sum of ep_s
      DOUBLE PRECISION Eps
!
!                      Source terms (Surface)
      DOUBLE PRECISION Sdp, Sdps
!
!                      Source terms (Volumetric)
      DOUBLE PRECISION V0, Vmt, Vbf
!-----------------------------------------------

      i = i_of(ijk)
      j = j_of(ijk)
      k = k_of(ijk)
      ijke = east_of(ijk)
!
      a = zero
      bb = zero
      eps = zero
!
      DO m = 1, mmax
         IF (close_packed(m)) THEN
            epga = avg_x(ep_s(ijk,m),ep_s(ijke,m),i)
!
!         Surface forces
!
!         Pressure term
            sdp = -p_scale*epga*(p_g(ijke)-p_g(ijk))*ayz(ijk)
            sdps = -epga*(p_s(ijke,m)-p_s(ijk,m))*ayz(ijk)
!
!           Shear stress terms
!
!         Volumetric forces
            ropga = avg_x(rop_s(ijk,m),rop_s(ijke,m),i)
!
!         Previous time step
            v0 = avg_x(rop_so(ijk,m),rop_so(ijke,m),i)*odt
!
!         Interphase mass transfer
            vmt = avg_x(sum_r_s(ijk,m),sum_r_s(ijke,m),i)
!
!         Body force
            vbf = ropga*bfx_s(ijk,m)
!
!         Collect the terms
            a = a - ((v0 + zmax(vmt))*vol_u(ijk))
            bb=bb-(sdp+sdps+((v0+zmax((-vmt)))*u_so(ijk,m)+vbf)*vol_u(ijk))
            eps = eps + epga
         ENDIF
      END DO
      b_m = -((bb - a)/(eps*ayz(ijk))-p_star(ijk))
!
      RETURN
      END SUBROUTINE b_m_p_star_e
!
!
!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: B_m_P_star_n(IJK, IER)                                 C
!  Purpose: Determine p_star when there is an interface at north       C
!                                                                      C
!                                                                      C
!  Author: M. Syamlal                                 Date: 6-AUG-96   C
!  Reviewer:                                          Date:            C
!                                                                      C
!                                                                      C
!  Literature/Document References:                                     C
!                                                                      C
!  Variables referenced:                                               C
!  Variables modified:                                                 C
!                                                                      C
!  Local variables:                                                    C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
!
      SUBROUTINE b_m_p_star_n(B_M, IJK)
!...Translated by Pacific-Sierra Research VAST-90 2.06G5  12:17:31  12/09/98
!...Switches: -xf
!
!  Include param.inc file to specify parameter values
!
!-----------------------------------------------
!   M o d u l e s
!-----------------------------------------------
      USE param
      USE param1
      USE scales
      USE constant
      USE physprop
      USE fldvar
      USE run
      USE rxns
      USE toleranc
      USE geometry
      USE indices
      USE compar
      USE bodyforce
      USE fun_avg
      USE functions
      IMPLICIT NONE
!-----------------------------------------------
!   G l o b a l   P a r a m e t e r s
!-----------------------------------------------
!-----------------------------------------------
!   D u m m y   A r g u m e n t s
!-----------------------------------------------
!
!                      Indices
      INTEGER          I, J, K, IJK, IJKN
!
!                      Phase index
      INTEGER          M
!
!                      Average volume fraction
      DOUBLE PRECISION EPGA
!
!                      Average density
      DOUBLE PRECISION ROPGA
!
!                      LHS (similar to A_m)
      DOUBLE PRECISION A
!
!                      RHS (similar to b_m)
      DOUBLE PRECISION BB
!
!                      b_m
      DOUBLE PRECISION b_m
!
!                      sum of ep_s
      DOUBLE PRECISION Eps
!
!                      Source terms (Surface)
      DOUBLE PRECISION Sdp, Sdps
!
!                      Source terms (Volumetric)
      DOUBLE PRECISION V0, Vmt, Vbf
!
!-----------------------------------------------

      i = i_of(ijk)
      j = j_of(ijk)
      k = k_of(ijk)
      ijkn = north_of(ijk)
!
      a = zero
      bb = zero
      eps = zero
!
      DO m = 1, mmax
         IF (close_packed(m)) THEN
            epga = avg_y(ep_s(ijk,m),ep_s(ijkn,m),j)
!
!         Surface forces
!
!         Pressure term
            sdp = -p_scale*epga*(p_g(ijkn)-p_g(ijk))*axz(ijk)
            sdps = -epga*(p_s(ijkn,m)-p_s(ijk,m))*axz(ijk)
!
!           Shear stress terms
!
!         Volumetric forces
            ropga = avg_y(rop_s(ijk,m),rop_s(ijkn,m),j)
!
!         Previous time step
            v0 = avg_y(rop_so(ijk,m),rop_so(ijkn,m),j)*odt
!
!         Interphase mass transfer
            vmt = avg_y(sum_r_s(ijk,m),sum_r_s(ijkn,m),j)
!
!         Body force
            vbf = ropga*bfy_s(ijk,m)
!
!         Collect the terms
            a = a - ((v0 + zmax(vmt))*vol_v(ijk))
            bb=bb-(sdp+sdps+((v0+zmax((-vmt)))*v_so(ijk,m)+vbf)*vol_v(ijk))
            eps = eps + epga
         ENDIF
      END DO
      b_m = -((bb - a)/(eps*axz(ijk))-p_star(ijk))
!
      RETURN
      END SUBROUTINE b_m_p_star_n
!
!
!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: B_m_P_star_t(IJK, IER)                                 C
!  Purpose: Determine p_star when there is an interface at top         C
!                                                                      C
!                                                                      C
!  Author: M. Syamlal                                 Date: 6-AUG-96   C
!  Reviewer:                                          Date:            C
!                                                                      C
!                                                                      C
!  Literature/Document References:                                     C
!                                                                      C
!  Variables referenced:                                               C
!  Variables modified:                                                 C
!                                                                      C
!  Local variables:                                                    C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
!
      SUBROUTINE b_m_p_star_t(B_M, IJK)
!...Translated by Pacific-Sierra Research VAST-90 2.06G5  12:17:31  12/09/98
!...Switches: -xf
!
!  Include param.inc file to specify parameter values
!
!-----------------------------------------------
!   M o d u l e s
!-----------------------------------------------
      USE param
      USE param1
      USE scales
      USE constant
      USE physprop
      USE fldvar
      USE run
      USE rxns
      USE toleranc
      USE geometry
      USE indices
      USE compar
      USE bodyforce
      USE fun_avg
      USE functions
      IMPLICIT NONE
!-----------------------------------------------
!   G l o b a l   P a r a m e t e r s
!-----------------------------------------------
!-----------------------------------------------
!   D u m m y   A r g u m e n t s
!-----------------------------------------------
!
!                      Indices
      INTEGER          I, J, K, IJK, IJKT
!
!                      Phase index
      INTEGER          M
!
!                      Average volume fraction
      DOUBLE PRECISION EPGA
!
!                      Average density
      DOUBLE PRECISION ROPGA
!
!                      LHS (similar to A_m)
      DOUBLE PRECISION A
!
!                      RHS (similar to b_m)
      DOUBLE PRECISION BB
!
!                      b_m
      DOUBLE PRECISION b_m
!
!                      sum of ep_s
      DOUBLE PRECISION Eps
!
!                      Source terms (Surface)
      DOUBLE PRECISION Sdp, Sdps
!
!                      Source terms (Volumetric)
      DOUBLE PRECISION V0, Vmt, Vbf
!
!-----------------------------------------------

      i = i_of(ijk)
      j = j_of(ijk)
      k = k_of(ijk)
      ijkt = top_of(ijk)
!
      a = zero
      bb = zero
      eps = zero
!
      DO m = 1, mmax
         IF (close_packed(m)) THEN
            epga = avg_z(ep_s(ijk,m),ep_s(ijkt,m),k)
!
!         Surface forces
!
!         Pressure term
            sdp = -p_scale*epga*(p_g(ijkt)-p_g(ijk))*axy(ijk)
            sdps = -epga*(p_s(ijkt,m)-p_s(ijk,m))*axy(ijk)
!
!           Shear stress terms
!
!         Volumetric forces
            ropga = avg_z(rop_s(ijk,m),rop_s(ijkt,m),k)
!
!         Previous time step
            v0 = avg_z(rop_so(ijk,m),rop_so(ijkt,m),k)*odt
!
!         Interphase mass transfer
            vmt = avg_z(sum_r_s(ijk,m),sum_r_s(ijkt,m),k)
!
!         Body force
            vbf = ropga*bfz_s(ijk,m)
!
!         Collect the terms
            a = a - ((v0 + zmax(vmt))*vol_w(ijk))
            bb=bb-(sdp+sdps+((v0+zmax((-vmt)))*w_so(ijk,m)+vbf)*vol_w(ijk))
            eps = eps + epga
         ENDIF
      END DO
      b_m = -((bb - a)/(eps*axy(ijk))-p_star(ijk))
!
      RETURN
      END SUBROUTINE b_m_p_star_t
!
!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: B_m_P_star_w(b_m, IJK, IER)                            C
!  Purpose: Determine p_star when there is an interface at west        C
!                                                                      C
!                                                                      C
!  Author: M. Syamlal                                 Date: 6-AUG-96   C
!  Reviewer:                                          Date:            C
!                                                                      C
!                                                                      C
!  Literature/Document References:                                     C
!                                                                      C
!  Variables referenced:                                               C
!  Variables modified:                                                 C
!                                                                      C
!  Local variables:                                                    C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
!
      SUBROUTINE b_m_p_star_w(B_M, IJK)
!...Translated by Pacific-Sierra Research VAST-90 2.06G5  12:17:31  12/09/98
!...Switches: -xf
!
!  Include param.inc file to specify parameter values
!
!-----------------------------------------------
!   M o d u l e s
!-----------------------------------------------
      USE param
      USE param1
      USE scales
      USE constant
      USE physprop
      USE fldvar
      USE run
      USE rxns
      USE toleranc
      USE geometry
      USE indices
      USE compar
      USE bodyforce
      USE fun_avg
      USE functions
      IMPLICIT NONE
!-----------------------------------------------
!   G l o b a l   P a r a m e t e r s
!-----------------------------------------------
!-----------------------------------------------
!   D u m m y   A r g u m e n t s
!-----------------------------------------------
!
!                      Indices
      INTEGER          IM, J, K, IJK, IJKW, IMJK
!
!                      Phase index
      INTEGER          M
!
!                      Average volume fraction
      DOUBLE PRECISION EPGA
!
!                      Average density
      DOUBLE PRECISION ROPGA
!
!                      LHS (similar to A_m)
      DOUBLE PRECISION A
!
!                      RHS (similar to b_m)
      DOUBLE PRECISION BB
!
!                      b_m
      DOUBLE PRECISION b_m
!
!                      sum of ep_s
      DOUBLE PRECISION Eps
!
!                      Source terms (Surface)
      DOUBLE PRECISION Sdp, Sdps
!
!                      Source terms (Volumetric)
      DOUBLE PRECISION V0, Vmt, Vbf
!
!-----------------------------------------------

      im = im1(i_of(ijk))
      j = j_of(ijk)
      k = k_of(ijk)
      ijkw = west_of(ijk)
      imjk = im_of(ijk)
!
      a = zero
      bb = zero
      eps = zero
!
      DO m = 1, mmax
         IF (close_packed(m)) THEN
            epga = avg_x(ep_s(ijkw,m),ep_s(ijk,m),im)
!
!         Surface forces
!
!         Pressure term
            sdp = -p_scale*epga*(p_g(ijk)-p_g(ijkw))*ayz(imjk)
            sdps = -epga*(p_s(ijk,m)-p_s(ijkw,m))*ayz(imjk)
!
!           Shear stress terms
!
!         Volumetric forces
            ropga = avg_x(rop_s(ijkw,m),rop_s(ijk,m),im)
!
!         Previous time step
            v0 = avg_x(rop_so(ijkw,m),rop_so(ijk,m),im)*odt
!
!         Interphase mass transfer
            vmt = avg_x(sum_r_s(ijkw,m),sum_r_s(ijk,m),im)
!
!         Body force
            vbf = ropga*bfx_s(ijkw,m)
!
!         Collect the terms
            a = a - ((v0 + zmax(vmt))*vol_u(ijkw))
            bb=bb-(sdp+sdps+((v0+zmax((-vmt)))*u_so(imjk,m)+vbf)*vol_u(ijkw))
            eps = eps + epga
         ENDIF
      END DO
      b_m = -(((-(bb - a)/(eps*ayz(imjk))))-p_star(ijk))
!
      RETURN
      END SUBROUTINE b_m_p_star_w
!
!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: B_m_P_star_s(b_m, IJK, IER)                            C
!  Purpose: Determine p_star when there is an interface at south       C
!                                                                      C
!                                                                      C
!  Author: M. Syamlal                                 Date: 6-AUG-96   C
!  Reviewer:                                          Date:            C
!                                                                      C
!                                                                      C
!  Literature/Document References:                                     C
!                                                                      C
!  Variables referenced:                                               C
!  Variables modified:                                                 C
!                                                                      C
!  Local variables:                                                    C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
!
      SUBROUTINE b_m_p_star_s(B_M, IJK)
!...Translated by Pacific-Sierra Research VAST-90 2.06G5  12:17:31  12/09/98
!...Switches: -xf
!
!  Include param.inc file to specify parameter values
!
!-----------------------------------------------
!   M o d u l e s
!-----------------------------------------------
      USE param
      USE param1
      USE scales
      USE constant
      USE physprop
      USE fldvar
      USE run
      USE rxns
      USE toleranc
      USE geometry
      USE indices
      USE compar
      USE bodyforce
      USE fun_avg
      USE functions
      IMPLICIT NONE
!-----------------------------------------------
!   G l o b a l   P a r a m e t e r s
!-----------------------------------------------
!-----------------------------------------------
!   D u m m y   A r g u m e n t s
!-----------------------------------------------
!
!
!                      Indices
      INTEGER          I, JM, K, IJK, IJKS, IJMK
!
!                      Phase index
      INTEGER          M
!
!                      Average volume fraction
      DOUBLE PRECISION EPGA
!
!                      Average density
      DOUBLE PRECISION ROPGA
!
!                      LHS (similar to A_m)
      DOUBLE PRECISION A
!
!                      RHS (similar to b_m)
      DOUBLE PRECISION BB
!
!                      b_m
      DOUBLE PRECISION b_m
!
!                      sum of ep_s
      DOUBLE PRECISION Eps
!
!                      Source terms (Surface)
      DOUBLE PRECISION Sdp, Sdps
!
!                      Source terms (Volumetric)
      DOUBLE PRECISION V0, Vmt, Vbf
!-----------------------------------------------

      i = i_of(ijk)
      jm = jm1(j_of(ijk))
      k = k_of(ijk)
      ijks = south_of(ijk)
      ijmk = jm_of(ijk)
!
      a = zero
      bb = zero
      eps = zero
!
      DO m = 1, mmax
         IF (close_packed(m)) THEN
            epga = avg_y(ep_s(ijks,m),ep_s(ijk,m),jm)
!
!         Surface forces
!
!         Pressure term
            sdp = -p_scale*epga*(p_g(ijk)-p_g(ijks))*axz(ijk)
            sdps = -epga*(p_s(ijk,m)-p_s(ijks,m))*axz(ijk)
!
!           Shear stress terms
!
!         Volumetric forces
            ropga = avg_y(rop_s(ijks,m),rop_s(ijk,m),jm)
!
!         Previous time step
            v0 = avg_y(rop_so(ijks,m),rop_so(ijk,m),jm)*odt
!
!         Interphase mass transfer
            vmt = avg_y(sum_r_s(ijks,m),sum_r_s(ijk,m),jm)
!
!         Body force
            vbf = ropga*bfy_s(ijks,m)
!
!         Collect the terms
            a = a - ((v0 + zmax(vmt))*vol_v(ijks))
            bb=bb-(sdp+sdps+((v0+zmax((-vmt)))*v_so(ijmk,m)+vbf)*vol_v(ijks))
            eps = eps + epga
         ENDIF
      END DO
      b_m = -(((-(bb - a)/(eps*axz(ijk))))-p_star(ijk))
!
      RETURN
      END SUBROUTINE b_m_p_star_s
!
!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: B_m_P_star_b(b_m, IJK, IER)                            C
!  Purpose: Determine p_star when there is an interface at bottom      C
!                                                                      C
!                                                                      C
!  Author: M. Syamlal                                 Date: 6-AUG-96   C
!  Reviewer:                                          Date:            C
!                                                                      C
!                                                                      C
!  Literature/Document References:                                     C
!                                                                      C
!  Variables referenced:                                               C
!  Variables modified:                                                 C
!                                                                      C
!  Local variables:                                                    C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
!
      SUBROUTINE b_m_p_star_b(B_M, IJK)
!...Translated by Pacific-Sierra Research VAST-90 2.06G5  12:17:31  12/09/98
!...Switches: -xf
!
!  Include param.inc file to specify parameter values
!
!-----------------------------------------------
!   M o d u l e s
!-----------------------------------------------
      USE param
      USE param1
      USE scales
      USE constant
      USE physprop
      USE fldvar
      USE run
      USE rxns
      USE toleranc
      USE geometry
      USE indices
      USE compar
      USE bodyforce
      USE fun_avg
      USE functions
      IMPLICIT NONE
!-----------------------------------------------
!   G l o b a l   P a r a m e t e r s
!-----------------------------------------------
!-----------------------------------------------
!   D u m m y   A r g u m e n t s
!-----------------------------------------------
!
!                      Indices
      INTEGER          I, J, KM, IJK, IJKB, IJKM
!
!                      Phase index
      INTEGER          M
!
!                      Average volume fraction
      DOUBLE PRECISION EPGA
!
!                      Average density
      DOUBLE PRECISION ROPGA
!
!                      LHS (similar to A_m)
      DOUBLE PRECISION A
!
!                      RHS (similar to b_m)
      DOUBLE PRECISION BB
!
!                      b_m
      DOUBLE PRECISION b_m
!
!                      sum of ep_s
      DOUBLE PRECISION Eps
!
!                      Source terms (Surface)
      DOUBLE PRECISION Sdp, Sdps
!
!                      Source terms (Volumetric)
      DOUBLE PRECISION V0, Vmt, Vbf
!-----------------------------------------------

      i = i_of(ijk)
      j = j_of(ijk)
      km = km1(k_of(ijk))
      ijkb = bottom_of(ijk)
      ijkm = km_of(ijk)
!
      a = zero
      bb = zero
      eps = zero
!
      DO m = 1, mmax
         IF (close_packed(m)) THEN
            epga = avg_z(ep_s(ijkb,m),ep_s(ijk,m),km)
!
!         Surface forces
!
!         Pressure term
            sdp = -p_scale*epga*(p_g(ijk)-p_g(ijkb))*axy(ijk)
            sdps = -epga*(p_s(ijk,m)-p_s(ijkb,m))*axy(ijk)
!
!           Shear stress terms
!
!         Volumetric forces
            ropga = avg_z(rop_s(ijkb,m),rop_s(ijk,m),km)
!
!         Previous time step
            v0 = avg_z(rop_so(ijkb,m),rop_so(ijk,m),km)*odt
!
!         Interphase mass transfer
            vmt = avg_z(sum_r_s(ijkb,m),sum_r_s(ijk,m),km)
!
!         Body force
            vbf = ropga*bfz_s(ijkb,m)
!
!         Collect the terms
            a = a - ((v0 + zmax(vmt))*vol_w(ijkb))
            bb=bb-(sdp+sdps+((v0+zmax((-vmt)))*w_so(ijkm,m)+vbf)*vol_w(ijkb))
            eps = eps + epga
         ENDIF
      END DO
      b_m = -(((-(bb - a)/(eps*axy(ijk))))-p_star(ijk))
!
      RETURN
      END SUBROUTINE b_m_p_star_b