drialg.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac3d\drialg.f
00002 !
00080                      SUBROUTINE DRIALG
00081 !                    *****************
00082 !
00083      & (FRV, FRT, RI, NPOIN3)
00084 !
00085 !***********************************************************************
00086 ! TELEMAC3D   V6P1                                   21/08/2010
00087 !***********************************************************************
00088 !
00089 !
00090 !
00091 !reference  "TURBULENCE MODELS FROM STATISTICAL PARAMETERS AT
00092 !+              ONE POINT : K-EPSILON MODEL"
00093 !+              E.D.F. BULLETIN DE LA DIRECTION DES ETUDES ET RECHERCHES -
00094 !+              SERIE A NUCLEAIRE, HYDRAULIQUE, THERMIQUE. NO1 1987
00095 !+              PP 35-47. (P.L. VIOLLET 'S PAPER)
00096 !reference "MODELE TRIDIMENSIONNEL D'ECOULEMENTS MESO-METEOROLOGIQUES.
00097 !+              ETUDE DE LA CONVECTION-DIFFUSION D'UN POLLUANT PASSIF A
00098 !+              CETTE ECHELLE"
00099 !+             (PP 37 TO 44 AND APPENDIX 1). J.C. SOLIVA 'S THESIS.
00100 !
00101 !
00102 !
00103 !
00104 !
00105 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00106 !| FRT            |<->| DAMPING FUNCTION FOR TRACERS
00107 !| FRV            |<->| DAMPING FUNCTION FOR VELOCITIES
00108 !| NPOIN3         |-->| NUMBER OF 3D POINTS
00109 !| RI             |<->| RICHARDSON NUMBER
00110 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00111 !
00112       IMPLICIT NONE
00113       INTEGER LNG,LU
00114       COMMON/INFO/LNG,LU
00115 !
00116 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00117 !
00118       INTEGER, INTENT(IN) :: NPOIN3
00119       DOUBLE PRECISION, INTENT(INOUT) :: FRV(NPOIN3), FRT(NPOIN3)
00120       DOUBLE PRECISION, INTENT(INOUT) :: RI(NPOIN3)
00121 !
00122 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00123 !
00124       INTEGER I
00125       DOUBLE PRECISION Q2, UW, TW, Q
00126 !
00127       INTRINSIC SQRT
00128 !
00129 !***********************************************************************
00130 !                            ____      ____      ____
00131 ! ADIMENSIONAL CORRELATIONS  T'W' AND  U'W' ( =  V'W' )
00132 ! (MAKES SURE THAT RI=0., FRI=1. FOR VELOCITIES)
00133 !
00134       DO I=1,NPOIN3
00135 !
00136 ! TURBULENT KINETIC ENERGY
00137 !
00138         IF(RI(I).LT.1.16D0) THEN
00139 !
00140       Q2=SQRT(4.507462D0*RI(I)**2+3.740367D-1*RI(I)+6.028062D-2)
00141      &   -2.412442D0*RI(I)+2.455211D-1
00142 !
00143       Q = SQRT(Q2)
00144 !                          ____
00145 ! ADIMENSIONAL CORRELATION T'W'
00146 !
00147       TW=(Q**3*(7.421047D-1*Q**2+
00148      &    2.147382D-1*RI(I)))/(9.999995D-1*Q**2*RI(I)-1.033492D0*Q**2+
00149      &    2.893635D-1*RI(I)**2-3.280470D0*RI(I))
00150 !                          ____
00151 ! ADIMENSIONAL CORRELATION U'W'
00152 !
00153       UW=(Q**3*(4.857410D-1*Q**2+1.541820D0*RI(I)))/(9.999995D-1*Q
00154      & **2*RI(I)-1.033492D0*Q**2+2.893635D-1*RI(I)**2-3.280470D0*RI(I))
00155 !
00156 ! DIVIDES BY -UW FOR RI=0. TO GET 1.
00157 ! IT IS INDEED SOUGHT THAT -UW (FOR RI=0.) BE EQUAL TO THE SQUARE OF
00158 ! THE KARMAN CONSTANT. IN VISCLM, THE MIXING LENGTH HAS ALREADY BEEN
00159 ! MULTIPLIED BY THE SQUARE OF THE KARMAN CONSTANT.
00160 !
00161           FRV(I) = -UW/0.161724D0
00162           FRT(I) = -TW/0.161724D0
00163 !
00164         ELSE
00165 !
00166           FRV(I) = 0.D0
00167           FRT(I) = 0.D0
00168 !
00169         ENDIF
00170 !
00171       ENDDO
00172 !
00173 !-----------------------------------------------------------------------
00174 !
00175       RETURN
00176       END SUBROUTINE DRIALG

Generated on Fri Aug 31 2013 18:12:58 by S.E.Bourban (HRW) using doxygen 1.7.0