ebcs4_vel_mod Module Reference

Functions/Subroutines
subroutine	ebcs4_vel (nseg, iseg, segvar, a, v, x, liste, nod, irect, la, fv, ms, stifn, ebcs, output, dt1, time)

Function/Subroutine Documentation

ebcs4_vel()

subroutine ebcs4_vel_mod::ebcs4_vel	(	integer	nseg,
		integer, dimension(nseg)	iseg,
		type(t_segvar)	segvar,
			a,
			v,
			x,
		integer, dimension(nod)	liste,
		integer	nod,
		integer, dimension(4,nseg)	irect,
			la,
			fv,
			ms,
			stifn,
		type(t_ebcs_vel), intent(in)	ebcs,
		type(output_), intent(inout)	output,
		intent(in)	dt1,
		intent(in)	time )

Parameters

[in,out]

output

output structure

Definition at line 40 of file ebcs4_vel.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE ebcs_mod
      USE segvar_mod
      USE output_mod , ONLY : output_
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NSEG,NOD,ISEG(NSEG),LISTE(NOD),IRECT(4,NSEG)
      my_real a(3,*),v(3,*),x(3,*),la(3,nod),ms(*),stifn(*),fv(*)
      TYPE(t_ebcs_vel), INTENT(IN) :: EBCS
      TYPE(t_segvar) :: SEGVAR
      my_real,INTENT(IN) :: dt1 !< time step
      my_real,INTENT(IN) :: time !< current time
      TYPE(OUTPUT_), INTENT(INOUT) :: OUTPUT !< output structure
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: I,IS,KSEG,N1,N2,N3,N4,NG1,NG2,NG3,NG4,N,IVX,IVY,IVZ,IRHO,IENER
      my_real :: orient,rho,c,roc,fac,
     .           x13,y13,z13,x24,y24,z24,nx,ny,nz,s,
     .           roou,enou,vmx,vmy,vmz,fluxi,fluxo,p,dvx,dvy,dvz,ener,
     .           vx,vy,vz
      my_real :: de_in, de_out, dm_in, dm_out
C-----------------------------------------------
      ivx=ebcs%ivx
      ivy=ebcs%ivy
      ivz=ebcs%ivz
      irho=ebcs%irho
      iener=ebcs%iener
      de_in = zero
      de_out = zero
      dm_in = zero
      dm_out = zero
 
      IF(ivx > 0)THEN
        vx=ebcs%vx*fv(ivx)
      ELSE
        vx=ebcs%vx
      ENDIF
      IF(ivy > 0)THEN
        vy=ebcs%vy*fv(ivy)
      ELSE
        vy=ebcs%vy
      ENDIF
      IF(ivz > 0)THEN
        vz=ebcs%vz*fv(ivz)
      ELSE
        vz=ebcs%vz
      ENDIF
       IF(irho > 0)THEN
        rho=ebcs%rho*fv(irho)
      ELSE
        rho=ebcs%rho
      ENDIF
      IF(iener > 0)THEN
        ener=ebcs%ener*fv(iener)
      ELSE
        ener=ebcs%ener
      ENDIF
      c=ebcs%c
      roc=rho*c
 
!INIT. MASS DENSITY AND ENERGY
! NORMAL AT NODES
      DO i=1,nod
         la(1,i)=zero
         la(2,i)=zero
         la(3,i)=zero
      ENDDO
 
      DO is=1,nseg
        kseg=abs(iseg(is))
        orient=float(iseg(is)/kseg)
        n1=irect(1,is)
        n2=irect(2,is)
        n3=irect(3,is)
        n4=irect(4,is)
        IF(n4==0 .OR. n4==n3) THEN
          fac=one_over_6*orient
          n4=n3
        ELSE
          fac=one_over_8*orient
        ENDIF
        ng1=liste(n1)
        ng2=liste(n2)
        ng3=liste(n3)
        ng4=liste(n4)
        x13=x(1,ng3)-x(1,ng1)
        y13=x(2,ng3)-x(2,ng1)
        z13=x(3,ng3)-x(3,ng1)
        x24=x(1,ng4)-x(1,ng2)
        y24=x(2,ng4)-x(2,ng2)
        z24=x(3,ng4)-x(3,ng2)
        nx=(y13*z24-z13*y24)*fac
        ny=(z13*x24-x13*z24)*fac
        nz=(x13*y24-y13*x24)*fac
        la(1,n1)=la(1,n1)+nx
        la(2,n1)=la(2,n1)+ny
        la(3,n1)=la(3,n1)+nz
        la(1,n2)=la(1,n2)+nx
        la(2,n2)=la(2,n2)+ny
        la(3,n2)=la(3,n2)+nz
        la(1,n3)=la(1,n3)+nx
        la(2,n3)=la(2,n3)+ny
        la(3,n3)=la(3,n3)+nz
        vmx=v(1,ng1)+v(1,ng2)+v(1,ng3)
        vmy=v(2,ng1)+v(2,ng2)+v(2,ng3)
        vmz=v(3,ng1)+v(3,ng2)+v(3,ng3)
        IF(n4/=n3) THEN
           la(1,n4)=la(1,n4)+nx
           la(2,n4)=la(2,n4)+ny
           la(3,n4)=la(3,n4)+nz
           vmx=vmx+v(1,ng4)
           vmy=vmy+v(2,ng4)
           vmz=vmz+v(3,ng4)
        ENDIF
! mass and energy balance
        roou = segvar%RHO(kseg)
        enou = segvar%EINT(kseg)
        fluxo=(vmx*nx+vmy*ny+vmz*nz)*dt1
        fluxi=min(fluxo,zero)
        fluxo=max(fluxo,zero)
        dm_out=dm_out-fluxo*roou
        dm_in=dm_in-fluxi*rho
        de_out=de_out-fluxo*enou
        de_in=de_in-fluxi*ener
! masse density and incoming energy stored in face buffer
        segvar%RHO(kseg)=rho
        segvar%EINT(kseg)=ener
      ENDDO
 
!$OMP CRITICAL
      output%DATA%INOUT%DM_IN  = output%DATA%INOUT%DM_IN + dm_in
      output%DATA%INOUT%DM_OUT = output%DATA%INOUT%DM_OUT + dm_out
      output%DATA%INOUT%DE_IN = output%DATA%INOUT%DE_IN + de_in
      output%DATA%INOUT%DE_OUT = output%DATA%INOUT%DE_OUT + de_out
!$OMP END CRITICAL
 
      IF(time == zero)THEN
        DO i=1,nod
          n=liste(i)
          v(1,n)=vx
          v(2,n)=vy
          v(3,n)=vz
        ENDDO
      ENDIF
 
      DO i=1,nod
        n=liste(i)
        s=sqrt(la(1,i)**2+la(2,i)**2+la(3,i)**2)
        dvx=v(1,n)-vx
        dvy=v(2,n)-vy
        dvz=v(3,n)-vz
        p=roc*(dvx*la(1,i)+dvy*la(2,i)+dvz*la(3,i))/s
        a(1,n)=a(1,n)-p*la(1,i)
        a(2,n)=a(2,n)-p*la(2,i)
        a(3,n)=a(3,n)-p*la(3,i)
        stifn(n)=stifn(n)+(two*(s*roc)**2)/ms(n)
      ENDDO
 
      RETURN