ale51_upwind2.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "vect01_c.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	ale51_upwind2 (pm, x, ixq, flux, flu1, ale_connect, itrimat, ddvol, qmv, iflg)

Function/Subroutine Documentation

ale51_upwind2()

subroutine ale51_upwind2	(		pm,
			x,
		integer, dimension(nixq,numelq)	ixq,
			flux,
			flu1,
		type(t_ale_connectivity), intent(in)	ale_connect,
		integer	itrimat,
			ddvol,
			qmv,
		integer	iflg )

Definition at line 33 of file ale51_upwind2.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE ale_connectivity_mod
      use element_mod , only : nixq
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "com04_c.inc"
#include      "vect01_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IXQ(NIXQ,NUMELQ), ITRIMAT,IFLG,IAD2
      my_real pm(npropm,nummat), flux(4,*), flu1(*),ddvol(*),qmv(8,*), x(3,numnod)
      TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECT
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J,II,MAT(MVSIZ)
      my_real reduc,qmv5,qmv6,qmv7,qmv8,
     .        flux1(mvsiz), flux2(mvsiz), flux3(mvsiz), flux4(mvsiz),
     .        upwl(4,mvsiz), r1, r2, r3, r4
C-----------------------------------------------
      DO i=lft,llt
        flux1(i) = flux(1,i)
        flux2(i) = flux(2,i)
        flux3(i) = flux(3,i)
        flux4(i) = flux(4,i)
      ENDDO
C-----------------------------------------------
      DO i=lft,llt
        ii=i+nft
        mat(i)=ixq(1,ii)
      ENDDO
C-----------------------------------------------
C     UPWIND
C-----------------------------------------------
      DO j=1,4
        DO i=lft,llt
          upwl(j,i)=pm(16,mat(i))
        ENDDO
      ENDDO
C
      !======================================================!
      !  BOUNDARY FACE : no volume flux by default           !
      !    slip wall bc                                      !
      !======================================================!
      DO i=lft,llt
       iad2 = ale_connect%ee_connect%iad_connect(i + nft) 
       reduc=pm(92,mat(i))
       ii=ale_connect%ee_connect%connected(iad2 + 1 - 1)
       IF(ii == 0)THEN
        flux1(i)=flux1(i)*reduc
       ENDIF
C
       ii=ale_connect%ee_connect%connected(iad2 + 2 - 1)
       IF(ii == 0)THEN
        flux2(i)=flux2(i)*reduc
       ENDIF
C
       ii=ale_connect%ee_connect%connected(iad2 + 3 - 1)
       IF(ii == 0)THEN
        flux3(i)=flux3(i)*reduc
       ENDIF
C
       ii=ale_connect%ee_connect%connected(iad2 + 4 - 1)
       IF(ii == 0)THEN
        flux4(i)=flux4(i)*reduc
       ENDIF
C
      ENDDO
C     
      IF(iflg == 1)THEN
         IF (n2d  ==  1) THEN
            !!! Planar integration of momentum when 2D axi
            DO i=lft,llt
               ii = i + nft
               flux(1,i) = flux1(i)-upwl(1,i)*abs(flux1(i))
               flux(2,i) = flux2(i)-upwl(2,i)*abs(flux2(i))
               flux(3,i) = flux3(i)-upwl(3,i)*abs(flux3(i))
               flux(4,i) = flux4(i)-upwl(4,i)*abs(flux4(i))
C     
               !!! Divide by mean radius
               qmv(5,i)  = flux1(i)+upwl(1,i)*abs(flux1(i))
               qmv(6,i)  = flux2(i)+upwl(2,i)*abs(flux2(i))
               qmv(7,i)  = flux3(i)+upwl(3,i)*abs(flux3(i))
               qmv(8,i)  = flux4(i)+upwl(4,i)*abs(flux4(i))
C                                   
               flu1(i) = qmv(5,i)  +  qmv(6,i)  +  qmv(7,i)  +  qmv(8,i)
 
               r1 = half * (x(2, ixq(2, ii)) + x(2, ixq(3, ii)))
               r2 = half * (x(2, ixq(3, ii)) + x(2, ixq(4, ii)))
               r3 = half * (x(2, ixq(4, ii)) + x(2, ixq(5, ii)))
               r4 = half * (x(2, ixq(5, ii)) + x(2, ixq(2, ii)))
               
               !!! Beware of r=0 axis
               IF (r1 /= zero) qmv(5,i) = qmv(5,i) / r1
               IF (r2 /= zero) qmv(6,i) = qmv(6,i) / r2
               IF (r3 /= zero) qmv(7,i) = qmv(7,i) / r3
               IF (r4 /= zero) qmv(8,i) = qmv(8,i) / r4
            ENDDO
         ELSE
            DO i=lft,llt
               flux(1,i) = flux1(i)-upwl(1,i)*abs(flux1(i))
               flux(2,i) = flux2(i)-upwl(2,i)*abs(flux2(i))
               flux(3,i) = flux3(i)-upwl(3,i)*abs(flux3(i))
               flux(4,i) = flux4(i)-upwl(4,i)*abs(flux4(i))
               
               qmv(5,i)  = flux1(i)+upwl(1,i)*abs(flux1(i))
               qmv(6,i)  = flux2(i)+upwl(2,i)*abs(flux2(i))
               qmv(7,i)  = flux3(i)+upwl(3,i)*abs(flux3(i))
               qmv(8,i)  = flux4(i)+upwl(4,i)*abs(flux4(i))
                 
               flu1(i) = qmv(5,i)  +  qmv(6,i)  +  qmv(7,i)  +  qmv(8,i)
            ENDDO
         ENDIF
      ELSE
         DO i=lft,llt
            flux(1,i) = flux1(i)-upwl(1,i)*abs(flux1(i))
            flux(2,i) = flux2(i)-upwl(2,i)*abs(flux2(i))
            flux(3,i) = flux3(i)-upwl(3,i)*abs(flux3(i))
            flux(4,i) = flux4(i)-upwl(4,i)*abs(flux4(i))
 
            qmv5  = flux1(i)+upwl(1,i)*abs(flux1(i))
            qmv6  = flux2(i)+upwl(2,i)*abs(flux2(i))
            qmv7  = flux3(i)+upwl(3,i)*abs(flux3(i))
            qmv8  = flux4(i)+upwl(4,i)*abs(flux4(i))
 
            flu1(i) = qmv5  +  qmv6  +  qmv7  +  qmv8
         ENDDO
      END IF
C
      IF(itrimat > 0)THEN
        DO i=lft,llt
          ddvol(i)=half * ( flu1(i)+flux(1,i)+flux(2,i)+flux(3,i)+flux(4,i) )
        ENDDO
      ENDIF
C----------------------------------------------- 
      RETURN