solid_drag.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Subroutine: SOLID_DRAG                                              C
!  Purpose: Accounting for the equal and opposite drag force on a      C
!  continuous solids phase due to discrete solid particles by          C
!  introducing the drag as a source term. Face centered.               C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
      SUBROUTINE solid_drag_u(A_M, B_M)

!-----------------------------------------------
! Modules
!-----------------------------------------------
      USE param
      USE param1
      USE geometry
      USE physprop
      USE indices
      USE compar
      USE discretelement
      USE fun_avg
      USE functions

      IMPLICIT NONE
!-----------------------------------------------
! Dummy Arguments
!-----------------------------------------------
! Septadiagonal matrix A_m
      DOUBLE PRECISION :: A_M(dimension_3, -3:3, 0:dimension_m)
! Vector b_m
      DOUBLE PRECISION :: B_M(dimension_3, 0:dimension_m)
!-----------------------------------------------
! Local Variables
!-----------------------------------------------
! Indices
      INTEGER :: IJK, IJKE, I
! Solids phase indices
      INTEGER :: M
!-----------------------------------------------

! currently no difference between interpolated and non-interpolated
! implementation of solid-solid drag

      DO m = 1, mmax
         DO ijk = ijkstart3, ijkend3
            IF(fluid_at(ijk)) THEN

               i = i_of(ijk)
               ijke = east_of(ijk)

               a_m(ijk,0,m) = a_m(ijk,0,m) - vol_u(ijk) *              &
                  avg_x(sdrag_am(ijk,m), sdrag_am(ijke,m), i)

               b_m(ijk,m) = b_m(ijk,m) - vol_u(ijk) *                  &
                  avg_x(sdrag_bm(ijk,1,m), sdrag_bm(ijke,1,m), i)

               ENDIF   ! end if (fluid_at(ijk))
         ENDDO   ! end do (ijk=ijkstart3,ijkend3)
      ENDDO   ! end do (cm=1,mmax)


      RETURN
      END SUBROUTINE solid_drag_u

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Subroutine: SOLID_DRAG                                              C
!  Purpose: Accounting for the equal and opposite drag force on a      C
!  continuous solids phase due to discrete solid particles by          C
!  introducing the drag as a source term. Face centered.               C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C

      SUBROUTINE solid_drag_v(A_M, B_M)

!-----------------------------------------------
! Modules
!-----------------------------------------------
      USE param
      USE param1
      USE geometry
      USE physprop
      USE indices
      USE compar
      USE discretelement
      USE fun_avg
      USE functions

      IMPLICIT NONE
!-----------------------------------------------
! Dummy Arguments
!-----------------------------------------------
! Septadiagonal matrix A_m
      DOUBLE PRECISION :: A_M(dimension_3, -3:3, 0:dimension_m)
! Vector b_m
      DOUBLE PRECISION :: B_M(dimension_3, 0:dimension_m)
!-----------------------------------------------
! Local Variables
!-----------------------------------------------
! Indices
      INTEGER :: IJK, IJKN, J
! Solids phase indices
      INTEGER :: M

      DO m = 1, mmax
         DO ijk = ijkstart3, ijkend3
            IF(fluid_at(ijk)) THEN
                  j = j_of(ijk)
                  ijkn = north_of(ijk)

                  a_m(ijk,0,m) = a_m(ijk,0,m) - vol_v(ijk) *           &
                     avg_y(sdrag_am(ijk,m), sdrag_am(ijkn,m), j)

                  b_m(ijk,m) = b_m(ijk,m) - vol_v(ijk) *               &
                     avg_y(sdrag_bm(ijk,2,m), sdrag_bm(ijkn,2,m), j)

            ENDIF   ! end if (fluid_at(ijk))
         ENDDO   ! end do (ijk=ijkstart3,ijkend3)
      ENDDO   ! end do (cm=1,mmax)


      RETURN
      END SUBROUTINE solid_drag_v


!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Subroutine: SOLID_DRAG_W                                            C
!  Purpose: Accounting for the equal and opposite drag force on a      C
!  continuous solids phase due to discrete solid particles by          C
!  introducing the drag as a source term. Face centered.               C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C

      SUBROUTINE solid_drag_w(A_M, B_M)

!-----------------------------------------------
! Modules
!-----------------------------------------------
      USE param
      USE param1
      USE geometry
      USE physprop
      USE indices
      USE compar
      USE discretelement
      USE fun_avg
      USE functions

      IMPLICIT NONE
!-----------------------------------------------
! Dummy Arguments
!-----------------------------------------------
! Septadiagonal matrix A_m
      DOUBLE PRECISION :: A_M(dimension_3, -3:3, 0:dimension_m)
! Vector b_m
      DOUBLE PRECISION :: B_M(dimension_3, 0:dimension_m)
!-----------------------------------------------
! Local Variables
!-----------------------------------------------
! Indices
      INTEGER :: IJK, IJKT, K
! Solids phase indices
      INTEGER :: M


      DO m = 1, mmax
         DO ijk = ijkstart3, ijkend3
            IF(fluid_at(ijk)) THEN
               ijkt = top_of(ijk)
               k = k_of(ijk)

               a_m(ijk,0,m) = a_m(ijk,0,m) + vol_w(ijk) *              &
                  avg_z(sdrag_am(ijk,m), sdrag_am(ijkt,m), k)
               b_m(ijk,m) = b_m(ijk,m) + vol_w(ijk) *                  &
                  avg_z(sdrag_bm(ijk,3,m), sdrag_bm(ijkt,3,m), k)

            ENDIF   ! end if (fluid_at(ijk))
         ENDDO   ! end do (ijk=ijkstart3,ijkend3)
      ENDDO   ! end do (cm=1,mmax)


      RETURN
      END SUBROUTINE solid_drag_w