calcstressheatflux.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	calcstressheatflux (sti, umel, gradvel, qfx, hcel, gradtel, nef, isti, iqfx, mi)

Function/Subroutine Documentation

calcstressheatflux()

subroutine calcstressheatflux	(	real*8, dimension(6,mi(1),*)	sti,
		real*8, dimension(*)	umel,
		real*8, dimension(3,3,*)	gradvel,
		real*8, dimension(3,mi(1),*)	qfx,
		real*8, dimension(*)	hcel,
		real*8, dimension(3,*)	gradtel,
		integer	nef,
		integer	isti,
		integer	iqfx,
		integer, dimension(*)	mi
	)

!
!     calculation of the stress and the heat flux for output
!     in CFD calculations
!
      implicit none
!
      integer nef,isti,iqfx,mi(*),i
!
      real*8 sti(6,mi(1),*),qfx(3,mi(1),*),umel(*),hcel(*),
     &  gradvel(3,3,*),gradtel(3,*)
!
      do i=1,nef
         if(isti.gt.0) then
            sti(1,1,i)=2.d0*umel(i)*gradvel(1,1,i)
            sti(2,1,i)=2.d0*umel(i)*gradvel(2,2,i)
            sti(3,1,i)=2.d0*umel(i)*gradvel(3,3,i)
            sti(4,1,i)=umel(i)*(gradvel(1,2,i)+gradvel(2,1,i))
            sti(5,1,i)=umel(i)*(gradvel(1,3,i)+gradvel(3,1,i))
            sti(6,1,i)=umel(i)*(gradvel(2,3,i)+gradvel(3,2,i))
         endif
         if(iqfx.gt.0) then
            qfx(1,1,i)=-hcel(i)*gradtel(1,i)
            qfx(2,1,i)=-hcel(i)*gradtel(2,i)
            qfx(3,1,i)=-hcel(i)*gradtel(3,i)
         endif
      enddo
!
      return