fldvar_mod.f

Go to the documentation of this file.

! -*- f90 -*-
!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module: fldvar                                                      C
!  Purpose: Common block containing field variables data               C
!                                                                      C
!  Author: M. Syamlal                                 Date: dd-mmm-yy  C
!  Reviewer:                                          Date: dd-mmm-yy  C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
      MODULE fldvar


! Void fraction
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  ep_g
! Store the void/volume fraction or value of 1 depending on setting of
! governing equations
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  epg_jfac
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  epg_ifac
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  eps_ifac

! Previous-time-step value of Void fraction
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  ep_go

! Gas pressure
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  p_g

! Previous-time-step value of Gas pressure
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  p_go

! Gas density
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  ro_g

! Previous-time-step value of Gas density
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  ro_go

! Macroscopic gas density
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  rop_g

! Previous-time-step value of macroscopic gas density
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  rop_go


! Solids phase densities
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: ro_s

! Previous-time-step value Solids phase densities
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: ro_so

! Macroscopic density of solids phases
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  rop_s

! Previous-time-step value of macroscopic density of solids phases
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  rop_so

! Macroscopic density of particle diameter
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  d_p

! Previous-time-step value of particle diameter
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  d_po

! Gas phase temperature
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  t_g

! Solids phase temperature
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  t_s

! Previous-time-step value of Gas phase temperature
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  t_go

! Previous-time-step value of Solids phase temperature
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  t_so

! Gas species mass fraction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  x_g

! Solids species mass fraction
      DOUBLE PRECISION, DIMENSION(:, :, :), ALLOCATABLE ::  x_s

! Previous-time-step value of Gas species mass fraction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  x_go

! Previous-time-step value of Solids species mass fraction
      DOUBLE PRECISION, DIMENSION(:, :, :), ALLOCATABLE ::  x_so

! x-component of gas velocity
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  u_g

! Previous-time-step value of x-component of gas velocity
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  u_go

! x-component of solids phase velocity
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  u_s

! Previous-time-step value of x-component of solids phase velocity
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  u_so

! y-component of gas velocity
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  v_g

! Previous time-step value of y-component of gas velocity
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  v_go

! y-component of solids phase velocity
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  v_s

! Previous time-step value of y-component of solids phase velocity
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  v_so

! z-component of gas velocity
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  w_g

! Previous time-step value of z-component of gas velocity
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  w_go

! z-component of solids phase velocity
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  w_s

! Previous time-step value of z-component of solids phase velocity
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  w_so

! Solids pressure as a result of granular motion
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  p_s

! Solids pressure as a result of granular motion
! Collisional Contribution
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  p_s_c

! Solids pressure as a result of granular motion
! Viscous Contribution
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  p_s_v

! Solids pressure as a result of granular motion
! Frictional Contribution
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  p_s_f

! Solids pressure as a result of granular motion
! Plastic Contribution
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  p_s_p

! Solids pressure that maintains EP_g >= EP_star
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  p_star

! Previous-time-step value of Solids pressure that maintains
! EP_g >= EP_star
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  p_staro

! Granular temperature of mth phase
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  theta_m

! Previous-time-step value of Granular temperature of mth phase
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  theta_mo

! User-defined Scalars
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  scalar

! Previous time-step value of User-defined Scalars
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  scalaro

! K and Epsilon for Gas turbulence
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  k_turb_g
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  e_turb_g

! Previous time-step value of K and Epsilon
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  k_turb_go
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  e_turb_go


!--------------------------------------------------------------------------> J.Musser TBR
      DOUBLE PRECISION, DIMENSION(:, :, :), ALLOCATABLE :: ep_ss
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: err_array


      CONTAINS

! Function to calculate the solids phase M volume fraction.
      DOUBLE PRECISION FUNCTION ep_s(IJK, xxM)
      IMPLICIT NONE
      INTEGER, INTENT(IN) :: IJK   ! spatial index
      INTEGER, INTENT(IN) :: xxM   ! solids phase index
      ep_s = rop_s(ijk, xxm)/ro_s(ijk,xxm)
      END FUNCTION


!!!HPF$ align EP_g(:) with TT(:)
!!!HPF$ align EP_go(:) with TT(:)
!!!HPF$ align P_g(:) with TT(:)
!!!HPF$ align P_go(:) with TT(:)
!!!HPF$ align RO_g(:) with TT(:)
!!!HPF$ align RO_go(:) with TT(:)
!!!HPF$ align ROP_g(:) with TT(:)
!!!HPF$ align ROP_go(:) with TT(:)
!!!HPF$ align D_p(:, *) with TT(:)
!!!HPF$ align D_po(:, *) with TT(:)
!!!HPF$ align ROP_s(:, *) with TT(:)
!!!HPF$ align ROP_so(:, *) with TT(:)
!!!HPF$ align T_g(:) with TT(:)
!!!HPF$ align T_s(:, *) with TT(:)
!!!HPF$ align T_go(:) with TT(:)
!!!HPF$ align T_so(:, *) with TT(:)
!!!HPF$ align X_g(:, *) with TT(:)
!!!HPF$ align X_s(:, *, *) with TT(:)
!!!HPF$ align X_go(:, *) with TT(:)
!!!HPF$ align X_so(:, *, *) with TT(:)
!!!HPF$ align U_g(:) with TT(:)
!!!HPF$ align U_go(:) with TT(:)
!!!HPF$ align U_s(:, *) with TT(:)
!!!HPF$ align U_so(:, *) with TT(:)
!!!HPF$ align V_g(:) with TT(:)
!!!HPF$ align V_go(:) with TT(:)
!!!HPF$ align V_s(:, *) with TT(:)
!!!HPF$ align V_so(:, *) with TT(:)
!!!HPF$ align W_g(:) with TT(:)
!!!HPF$ align W_go(:) with TT(:)
!!!HPF$ align W_s(:, *) with TT(:)
!!!HPF$ align W_so(:, *) with TT(:)
!!!HPF$ align P_s(:, *) with TT(:)
!!!HPF$ align P_s_c(:, *) with TT(:)
!!!HPF$ align P_star(:) with TT(:)
!!!HPF$ align P_staro(:) with TT(:)
!!!HPF$ align THETA_m(:, *) with TT(:)
!!!HPF$ align THETA_mo(:, *) with TT(:)


      END MODULE fldvar