m6law.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"
#include "com08_c.inc"
#include "param_c.inc"
#include "scr14_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	m6law (pm, off, sig, eint, rho, rk, re, vorti, vnew, rho0, vis, psh, d1, d2, d3, d4, d5, d6, c1, c2, c3, c4, c5, c6, pc, amu, amu2, espe, ssp, df, tmu, mat, wxx, wyy, wzz, nel, dpdm, jtur, jpor)

Function/Subroutine Documentation

m6law()

subroutine m6law	(		pm,
			off,
			sig,
			eint,
			rho,
			rk,
			re,
			vorti,
			vnew,
			rho0,
			vis,
			psh,
			d1,
			d2,
			d3,
			d4,
			d5,
			d6,
			c1,
			c2,
			c3,
			c4,
			c5,
			c6,
			pc,
			amu,
			amu2,
			espe,
			ssp,
			df,
			tmu,
		integer, dimension(*)	mat,
			wxx,
			wyy,
			wzz,
		integer	nel,
			dpdm,
		integer, intent(in)	jtur,
		integer, intent(in)	jpor )

Definition at line 28 of file m6law.F.

C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "com08_c.inc"
#include      "param_c.inc"
#include      "scr14_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: JTUR
      INTEGER, INTENT(IN) :: JPOR
      INTEGER MAT(*),NEL
      my_real
     .   pm(npropm,*), off(*), sig(nel,6), eint(*), rho(*), rk(*), re(*),
     .   vorti(*)
      my_real
     .   vnew(*), rho0(*), d1(*), d2(*), d3(*), d4(*), d5(*), d6(*),
     .   pc(*), vis(*), c1(*), c2(*), c3(*), c4(*), c5(*),
     .   c6(*), amu(*), amu2(*), espe(*), ssp(*), df(*),
     .   tmu(*), psh(*), wxx(*), wyy(*), wzz(*), rho0_1,
     .   pc_1, c1_1, c2_1, c3_1,
     .   c4_1, c5_1, c6_1, tmu_1, psh_1
      my_real
     .     dpdm(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, MX
      my_real
     .   dav, rk2t, vis2, fac
C-----------------------------------------------
 
      !SAUVEGARDE VORTICITE SI ANIM
      IF((anim_e(10)==1 .OR. anim_se(10)==1).AND. dt1/=zero)THEN
         fac=four/dt1
        IF(n2d==0)THEN
          DO 5 i=1,nel
 5        vorti(i)=fac*sqrt(wxx(i)**2+wyy(i)**2+wzz(i)**2)
        ELSE
          DO 6 i=1,nel
 6        vorti(i)=fac*wzz(i)
        ENDIF
      ENDIF
 
      !EOS COEFFICIENTS
      mx     =mat(1)
      rho0_1=pm( 1,mx)
      pc_1  =pm(37,mx)
      c1_1  =pm(104,mx)-pm(88,mx)
      c2_1  =pm(32,mx)
      c3_1  =pm(33,mx)
      c4_1  =pm(34,mx)
      c5_1  =pm(35,mx)
      c6_1  =pm(36,mx)
      tmu_1 =pm(81,mx)
      psh_1=pm(88,mx)
      DO 10 i=1,nel
      rho0(i)=rho0_1
      vis(i) =pm(24,mx)*rho(i)
      pc(i)  =pc_1
      c1(i)  =c1_1
      c2(i)  =c2_1
      c3(i)  =c3_1
      c4(i)  =c4_1
      c5(i)  =c5_1
      c6(i)  =c6_1
      tmu(i) =tmu_1
      psh(i)=psh_1
   10 CONTINUE
 
      !IF TURBULENT
      IF(jtur/=0)THEN
        IF(jpor/=2)THEN
          DO i=1,nel
            tmu(i)=tmu(i)*rk(i)*rk(i) / max(em15,(re(i)*vnew(i)))
            vis(i)=vis(i)+tmu(i)
          ENDDO
        ELSE
          DO i=1,nel
            tmu(i)=zero
          ENDDO
        ENDIF
      ELSE
        DO i=1,nel
          rk(i)=zero
        ENDDO
      ENDIF
 
      !OUTPUT
      DO i=1,nel
       vis2 = two*vis(i)
       dav  = -(d1(i)+d2(i)+d3(i))*third
       !------------------------------!
       !   VISCOUS STRESS TENSOR      !
       !------------------------------!
       sig(i,1)=vis2*(d1(i)+dav)
       sig(i,2)=vis2*(d2(i)+dav)
       sig(i,3)=vis2*(d3(i)+dav)
       sig(i,4)=vis(i) *d4(i)
       sig(i,5)=vis(i) *d5(i)
       sig(i,6)=vis(i) *d6(i)
 
C         DF(I)  =RHO0(I)/RHO(I)
C         AMU(I) =ONE/DF(I)-ONE
       amu2(i)= max(zero,amu(i))**2
       espe(i)=df(i)*eint(i) / max(em15,vnew(i))
       rk2t = two*rk(i) / max(em15,three*vnew(i))
       
       !------------------------------!
       !         SOUND SPEED          !
       !------------------------------!
       
       dpdm(i) = dpdm(i)
     .       +(c5(i)+c6(i)*amu(i))*df(i)*df(i)*rk2t + rk2t*df(i)
       ssp(i)=sqrt(abs(dpdm(i))/rho0(i))
      ENDDO
 
      RETURN