paramTurb.F90

Go to the documentation of this file.

module paramturb
!
!       Module that contains the constants for the turbulence models
!       as well as some global variables/parameters for the turbulent
!       routines.
!
    use constants, only: realtype, inttype
    use inputphysics, only: rsak => sakappa, rsacb1 => sacb1, rsacb2 => sacb2, &
                            rsacb3 => sasigma, rsacv1 => sacv1, rsacw2 => sacw2, &
                            rsacw3 => sacw3, rsact1 => sact1, rsact2 => sact2, &
                            rsact3 => sact3, rsact4 => sact4, rsacrot => sacrot
    implicit none
    save
!
!       Spalart-Allmaras constants. rsaK, rsaCb1, rsaCb2, rsaCb3, rsaCv1,
!       rsaCw2, rsaCw3, rsaCt1, rsaCt2, rsaCt3, rsaCt4, rsaCrot are imported
!       directly from inputPhysics via USE renaming above.
!       rsaCw1 is derived: cb1/kappa^2 + (1+cb2)/sigma, computed at init.
!
    real(kind=realtype) :: rsacw1
!
!       K-omega constants.
!
    real(kind=realtype), parameter :: rkwk = 0.41_realtype
    real(kind=realtype), parameter :: rkwsigk1 = 0.5_realtype
    real(kind=realtype), parameter :: rkwsigw1 = 0.5_realtype
    real(kind=realtype), parameter :: rkwsigd1 = 0.5_realtype
    real(kind=realtype), parameter :: rkwbeta1 = 0.0750_realtype
    real(kind=realtype), parameter :: rkwbetas = 0.09_realtype
!
!       K-omega SST constants.
!
    real(kind=realtype), parameter :: rsstk = 0.41_realtype
    real(kind=realtype), parameter :: rssta1 = 0.31_realtype
    real(kind=realtype), parameter :: rsstbetas = 0.09_realtype
 
    real(kind=realtype), parameter :: rsstsigk1 = 0.85_realtype
    real(kind=realtype), parameter :: rsstsigw1 = 0.5_realtype
    real(kind=realtype), parameter :: rsstbeta1 = 0.0750_realtype
 
    real(kind=realtype), parameter :: rsstsigk2 = 1.0_realtype
    real(kind=realtype), parameter :: rsstsigw2 = 0.856_realtype
    real(kind=realtype), parameter :: rsstbeta2 = 0.0828_realtype
!
!       K-tau constants.
!
    real(kind=realtype), parameter :: rktk = 0.41_realtype
    real(kind=realtype), parameter :: rktsigk1 = 0.5_realtype
    real(kind=realtype), parameter :: rktsigt1 = 0.5_realtype
    real(kind=realtype), parameter :: rktsigd1 = 0.5_realtype
    real(kind=realtype), parameter :: rktbeta1 = 0.0750_realtype
    real(kind=realtype), parameter :: rktbetas = 0.09_realtype
!
!       V2-f constants.
!
    real(kind=realtype), parameter :: rvfc1 = 1.4_realtype
    real(kind=realtype), parameter :: rvfc2 = 0.3_realtype
    real(kind=realtype), parameter :: rvfbeta = 1.9_realtype
    real(kind=realtype), parameter :: rvfsigk1 = 1.0_realtype
    real(kind=realtype), parameter :: rvfsige1 = 0.7692307692_realtype
    real(kind=realtype), parameter :: rvfsigv1 = 1.00_realtype
    real(kind=realtype), parameter :: rvfcn = 70.0_realtype
 
    real(kind=realtype), parameter :: rvfn1cmu = 0.190_realtype
    real(kind=realtype), parameter :: rvfn1a = 1.300_realtype
    real(kind=realtype), parameter :: rvfn1b = 0.250_realtype
    real(kind=realtype), parameter :: rvfn1cl = 0.300_realtype
    real(kind=realtype), parameter :: rvfn6cmu = 0.220_realtype
    real(kind=realtype), parameter :: rvfn6a = 1.400_realtype
    real(kind=realtype), parameter :: rvfn6b = 0.045_realtype
    real(kind=realtype), parameter :: rvfn6cl = 0.230_realtype
 
    real(kind=realtype) :: rvflimitk, rvflimite, rvfcl
    real(kind=realtype) :: rvfcmu
!
!       Variables to store the parameters for the wall functions fits.
!       As these variables depend on the turbulence model they are set
!       during runtime. Allocatables are used, because the number of
!       fits could be different for the different models.
!       The curve is divided in a number of intervals and is
!       constructed such that both the function and the derivatives
!       are continuous. Consequently cubic polynomials are used.
!
    ! nFit:               Number of intervals of the curve.
    ! ypT(0:nFit):        y+ values at the interval boundaries.
    ! reT(0:nFit):        Reynolds number at the interval
    !                     boundaries, where the Reynolds number is
    !                     defined with the local velocity and the
    !                     wall distance.
    ! up0(nFit):          Coefficient 0 in the fit for the
    !                     nondimensional tangential velocity as a
    !                     function of the Reynolds number.
    ! up1(nFit):          Idem for coefficient 1.
    ! up2(nFit):          Idem for coefficient 2.
    ! up3(nFit):          Idem for coefficient 3.
    ! tup0(nFit,nt1:nt2): Coefficient 0 in the fit for the
    !                     nondimensional turbulence variables as a
    !                     function of y+.
    ! tup1(nFit,nt1:nt2): Idem for coefficient 1.
    ! tup2(nFit,nt1:nt2): Idem for coefficient 2.
    ! tup3(nFit,nt1:nt2): Idem for coefficient 3.
    ! tuLogFit(nt1:nt2):  Whether or not the logarithm of the variable
    !                     has been fitted.
 
    integer(kind=intType) :: nfit
 
    real(kind=realtype), dimension(:), allocatable :: ypt, ret
    real(kind=realtype), dimension(:), allocatable :: up0, up1
    real(kind=realtype), dimension(:), allocatable :: up2, up3
 
    real(kind=realtype), dimension(:, :), allocatable :: tup0, tup1
    real(kind=realtype), dimension(:, :), allocatable :: tup2, tup3
#ifndef USE_TAPENADE
    real(kind=realtype), dimension(:), allocatable :: yptb, retb
    real(kind=realtype), dimension(:), allocatable :: up0b, up1b
    real(kind=realtype), dimension(:), allocatable :: up2b, up3b
#endif
 
    logical, dimension(:), allocatable :: tulogfit
 
end module paramturb