calcguesstincf.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	calcguesstincf (nface, dmin, vfa, umfa, cvfa, hcfa, ithermal, tincfguess, compressible)

Function/Subroutine Documentation

calcguesstincf()

subroutine calcguesstincf	(	integer	nface,
		real*8	dmin,
		real*8, dimension(0:7,*)	vfa,
		real*8, dimension(*)	umfa,
		real*8, dimension(*)	cvfa,
		real*8, dimension(*)	hcfa,
		integer	ithermal,
		real*8	tincfguess,
		integer	compressible
	)

!
!     calculates a guess for tincf based on the minimum of:
!       dmin/(local velocity)
!       density*dmin**2/(2*dynamic_viscosity)
!
!       dmin is the smallest edge length of the mesh
!
      implicit none
!
      integer nface,i,ithermal,compressible
!
      real*8 vfa(0:7,*),umax,dmin,umfa(*),tincfguess,cvfa(*),hcfa(*)
!
      tincfguess=1.d30
      do i=1,nface
!
!        convection
!
         umax=dsqrt(vfa(1,i)*vfa(1,i)+
     &              vfa(2,i)*vfa(2,i)+
     &              vfa(3,i)*vfa(3,i))
         if(umax.gt.1.d-30) tincfguess=min(tincfguess,dmin/umax)
!
!        viscous diffusion
!
         if((umfa(i).gt.0.d0).and.(vfa(5,i).gt.0.d0)) then
            tincfguess=min(tincfguess,vfa(5,i)*dmin*dmin/
     &           (2.d0*umfa(i)))
         endif
!
!        thermal diffusion
!
         if(ithermal.gt.0) then
            if((hcfa(i).gt.0.d0).and.(cvfa(i).gt.0.d0).and.
     &         (vfa(5,i).gt.0.d0)) then
               tincfguess=min(tincfguess,vfa(5,i)*cvfa(i)*dmin*dmin/
     &              (2.d0*hcfa(i)))
            endif
         endif
      enddo
!     
      return