calcdifflux_kdif.f90

国外大名顶顶的“台风”并行计算流体力学CFD软件的早期版本的源代码

!------------------------------------------------------------------------------!
! Procedure : calcdifflux_kdif            Auteur : E. Radenac
!                                         Date   : Juillet 2003
! Fonction                                Modif  : 
!  Calcul de la difference de flux a appliquer lorsque l'echange entre deux
!  zones n'est pas a chaque pas de temps (correction de flux). Solver thermique
!
! Defauts/Limitations/Divers :
!------------------------------------------------------------------------------!
 
subroutine calcdifflux_kdif(etatcons1, etatcons2, nfacelim, corcoef, connface2)

use OUTPUT
use DEFZONE
use DEFFIELD
use GEO3D
use TYPHMAKE

implicit none

! -- Declaration des entrees --
integer                    :: nfacelim            ! nombre de faces limites
real(krp), dimension(nfacelim) &
                           :: corcoef             ! coeff correction de flux
integer, dimension(nfacelim) &
                           :: connface2

! -- Declaration des entrees/sorties --
type(st_scafield), dimension(3) &
                           :: etatcons1, etatcons2 ! stockage des flux cumules
                                                   ! et des differences de flux
                                                   ! pour les deux zones

! -- Declaration des variables internes --
integer                        :: i              ! indice de face
real(krp)                      :: dif_enflux     ! difference des energies d'interface dans les deux zones

! -- Debut de la procedure --

! Supplement de flux pour les echanges espaces : calcul de la difference a appliquer

do i=1, nfacelim

! La difference est la "somme" des flux cumules car ce sont les valeurs algebriques dont on dispose
! (les flux sortant de part et d'autre)

dif_enflux = etatcons2(1)%scal(connface2(i)) + etatcons1(1)%scal(i)

  ! difference de flux pour la zone1 : tient compte du reste du cycle precedent
  ! (etatcons(3))
  etatcons1(2)%scal(i) = -corcoef(i)*dif_enflux + etatcons1(3)%scal(i)
  ! initialisation du reste a la difference initiale
  etatcons1(3)%scal(i) = etatcons1(2)%scal(i) 
  etatcons2(2)%scal(connface2(i)) = (-1._krp + corcoef(i))*dif_enflux + &
                                    etatcons2(3)%scal(connface2(i))
  ! initialisation du reste a la difference initiale
  etatcons2(3)%scal(connface2(i)) = etatcons2(2)%scal(connface2(i))

enddo

endsubroutine calcdifflux_kdif

!------------------------------------------------------------------------------!
! Historique des modifications
!
! juillet 2003 (v0.0.1b): creation de la procedure
! oct 2003              : ajout coefficient correction de flux
!------------------------------------------------------------------------------!