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, pc, amu, amu2, espe, ssp, df, tmu, mat, wxx, wyy, wzz, nel, dpdm, jtur, jpor, eos_struct)

Function/Subroutine Documentation

m6law()

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

Definition at line 30 of file m6law.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
        USE eos_param_mod , ONLY : eos_param_
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(*),
     .   amu(*), amu2(*), espe(*), ssp(*), df(*),
     .   tmu(*), psh(*), wxx(*), wyy(*), wzz(*), rho0_1,
     .   pc_1, tmu_1, psh_1
      my_real, INTENT(IN) :: dpdm(nel)
      TYPE(EOS_PARAM_),INTENT(INOUT) :: EOS_STRUCT
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)
      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
      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)
       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          !
       !------------------------------!
       ssp(i)=sqrt(abs(dpdm(i))/rho0(i))
      ENDDO
 
      RETURN