qmomk_initial_conditions.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: QMOMK_INITIAL_CONDITIONS                               C
!  Author: Alberto Passalacqua                        Date: 30-Jul-08  C
!  Reviewer:                                          Date:            C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C

      SUBROUTINE qmomk_initial_conditions

!-----------------------------------------------
!     M o d u l e s
!-----------------------------------------------
      USE param
      USE param1
      USE constant
      USE physprop
      USE fldvar
      USE geometry
      USE compar
      USE indices
      USE qmom_kinetic_equation
      USE qmomk_quadrature
      USE qmomk_parameters
      USE ic
      USE functions

      IMPLICIT NONE

      DOUBLE PRECISION InitVal
      INTEGER :: I, M, IJK

      DO ijk = ijkstart3, ijkend3

        IF(.NOT.fluid_at(ijk)) cycle
        DO m = 1, mmax
          ! A.P.
          ! Initializing weights as volume_fraction/number_of_nodes
          ! Note that MFIX doesn't initialize the volume fraction in
          ! boundaries, and QMOM need a positive value there.
          ! Here the boundaries are not considered, and they will
          ! be initialized separetely in the correspoding module

          DO i = 1, qmomk_nn
             qmomk_n0(i, ijk, m) = rop_s(ijk, m)/(qmomk_nn * ro_s(ijk,m))
          END DO

          ! A.P. Granular temperature minimum value is bounded
          !      This is required by QMOM
          initval = max(sqrt(theta_m(ijk,m)), minimum_theta)

          qmomk_u0(1, ijk, m) = -initval + u_s(ijk, m)
          qmomk_u0(2, ijk, m) = +initval + u_s(ijk, m)
          qmomk_u0(3, ijk, m) = -initval + u_s(ijk, m)
          qmomk_u0(4, ijk, m) = +initval + u_s(ijk, m)
          qmomk_u0(5, ijk, m) = -initval + u_s(ijk, m)
          qmomk_u0(6, ijk, m) = +initval + u_s(ijk, m)
          qmomk_u0(7, ijk, m) = -initval + u_s(ijk, m)
          qmomk_u0(8, ijk, m) = +initval + u_s(ijk, m)

          qmomk_v0(1, ijk, m) = -initval + v_s(ijk, m)
          qmomk_v0(2, ijk, m) = -initval + v_s(ijk, m)
          qmomk_v0(3, ijk, m) = +initval + v_s(ijk, m)
          qmomk_v0(4, ijk, m) = +initval + v_s(ijk, m)
          qmomk_v0(5, ijk, m) = -initval + v_s(ijk, m)
          qmomk_v0(6, ijk, m) = -initval + v_s(ijk, m)
          qmomk_v0(7, ijk, m) = +initval + v_s(ijk, m)
          qmomk_v0(8, ijk, m) = +initval + v_s(ijk, m)

          qmomk_w0(1, ijk, m) = -initval + w_s(ijk, m)
          qmomk_w0(2, ijk, m) = -initval + w_s(ijk, m)
          qmomk_w0(3, ijk, m) = -initval + w_s(ijk, m)
          qmomk_w0(4, ijk, m) = -initval + w_s(ijk, m)
          qmomk_w0(5, ijk, m) = +initval + w_s(ijk, m)
          qmomk_w0(6, ijk, m) = +initval + w_s(ijk, m)
          qmomk_w0(7, ijk, m) = +initval + w_s(ijk, m)
          qmomk_w0(8, ijk, m) = +initval + w_s(ijk, m)
        END DO
      END DO

      DO ijk = ijkstart3, ijkend3
       DO m = 1, mmax
         IF (fluid_at(ijk)) THEN
           CALL moments_twenty_eight_nodes (qmomk_n0(:,ijk,m), &
                qmomk_u0(:,ijk,m), qmomk_v0(:,ijk,m), qmomk_w0(:,ijk,m), qmomk_m0(:,ijk,m))

           CALL eight_node_3d (qmomk_m0(:,ijk,m), qmomk_n0(:,ijk,m), &
                qmomk_u0(:,ijk,m), qmomk_v0(:,ijk,m), qmomk_w0(:,ijk,m))

           CALL moments_twenty_eight_nodes (qmomk_n0(:,ijk,m), &
                qmomk_u0(:,ijk,m), qmomk_v0(:,ijk,m), qmomk_w0(:,ijk,m), qmomk_m0(:,ijk,m))
         END IF
       END DO
      END DO

      END SUBROUTINE qmomk_initial_conditions