← Back to branch summary

~ubuntu-branches/ubuntu/precise/code-saturne/precise

« back to all changes in this revision

Viewing changes to include/base/cstphy.h

  • Committer: Package Import Robot
  • Author(s): Sylvestre Ledru
  • Date: 2011-11-24 00:00:08 UTC
  • mfrom: (6.1.9 sid)
  • Revision ID: package-import@ubuntu.com-20111124000008-2vo99e38267942q5
Tags: 2.1.0-3
Install a missing file

Show diffs side-by-side

added added

removed removed

Lines of Context:
include/base/cstphy.h
1
 
!-------------------------------------------------------------------------------
2
 
 
3
 
!     This file is part of the Code_Saturne Kernel, element of the
4
 
!     Code_Saturne CFD tool.
5
 
 
6
 
!     Copyright (C) 1998-2009 EDF S.A., France
7
 
 
8
 
!     contact: saturne-support@edf.fr
9
 
 
10
 
!     The Code_Saturne Kernel is free software; you can redistribute it
11
 
!     and/or modify it under the terms of the GNU General Public License
12
 
!     as published by the Free Software Foundation; either version 2 of
13
 
!     the License, or (at your option) any later version.
14
 
 
15
 
!     The Code_Saturne Kernel is distributed in the hope that it will be
16
 
!     useful, but WITHOUT ANY WARRANTY; without even the implied warranty
17
 
!     of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
18
 
!     GNU General Public License for more details.
19
 
 
20
 
!     You should have received a copy of the GNU General Public License
21
 
!     along with the Code_Saturne Kernel; if not, write to the
22
 
!     Free Software Foundation, Inc.,
23
 
!     51 Franklin St, Fifth Floor,
24
 
!     Boston, MA  02110-1301  USA
25
 
 
26
 
!-------------------------------------------------------------------------------
27
 
 
28
 
!                              cstphy.h
29
 
!===============================================================================
30
 
 
31
 
! KELVIN
32
 
 
33
 
!       TKELVI       --> =  273,15
34
 
!       TKELVN       --> = -273,15
35
 
 
36
 
double precision   tkelvi            , tkelvn
37
 
parameter (        tkelvi =  273.15d0, tkelvn = -273.15d0)
38
 
 
39
 
! CALORIES
40
 
 
41
 
!       1 cal = XCAL2J J
42
 
 
43
 
double precision   xcal2j
44
 
parameter (        xcal2j = 4.1855d0)
45
 
 
46
 
! STEPHAN BOLTZMANN
47
 
 
48
 
double precision   stephn
49
 
parameter (        stephn = 5.6703d-8)
50
 
 
51
 
! GRAVITE
52
 
 
53
 
double precision   gx,gy,gz
54
 
common / rgravi /  gx,gy,gz
55
 
 
56
 
! Vecteur rotation
57
 
 
58
 
integer            icorio
59
 
common / irotat /  icorio
60
 
 
61
 
double precision   omegax, omegay, omegaz
62
 
double precision   irot(3,3), prot(3,3), qrot(3,3), rrot(3,3)
63
 
common / rrotat /  omegax, omegay, omegaz, &
64
 
                   irot, prot  , qrot  , rrot
65
 
 
66
 
! CONSTANTES PHYSIQUES DU FLUIDE
67
 
!   IXYZP0 : INDICATEUR DE REMPLISSAGE DE XYZP0
68
 
!   RO0    : MASSE VOLUMIQUE    DE REFERENCE
69
 
!   VISCL0 : VISCOSITE          DE REFERENCE
70
 
!   P0     : PRESSION TOTALE    DE REFERENCE
71
 
!   PRED0  : PRESSION REDUITE   DE REFERENCE
72
 
!   XYZP0  : POSITION PRESSION  DE REFERENCE
73
 
!   T0     : TEMPERATURE        DE REFERENCE
74
 
!   CP0    : CHALEUR SPECIFIQUE DE REFERENCE
75
 
 
76
 
integer           ixyzp0(nphsmx)
77
 
common / icstfl / ixyzp0
78
 
double precision  ro0(nphsmx)    , viscl0(nphsmx),                &
79
 
                  p0 (nphsmx)    , pred0 (nphsmx),                &
80
 
                  xyzp0(3,nphsmx), t0    (nphsmx),                &
81
 
                  cp0(nphsmx)
82
 
common / rcstfl / ro0        , viscl0        ,                    &
83
 
                  p0         , pred0         ,                    &
84
 
                  xyzp0      , t0            ,                    &
85
 
                  cp0
86
 
 
87
 
! TURBULENCE
88
 
!   IVISLS = 0 : VISCOSITE LAMINAIRE CONSTANTE = VISLS0
89
 
!   XKAPPA : CST DE KARMAN (~0.42)
90
 
!   CSTLOG : CST DE LA LOI LOG: 1/XKAPPA*LOG(YPLUS) + CSTLOG (~5.2)
91
 
!   YPLULI : YPLUS LIMITE 1./XKAPPA OU 10.88 SI IDEUCH=2
92
 
!   *POW   : COEFF WERNER AND WENGLE
93
 
!   CMU025 = CMU**0.25
94
 
!   CE1, CE2, SIGMAK, SIGMAE :
95
 
!            CONSTANTES DU K-EPSILON
96
 
!   C*RIJ* : CONSTANTES DU Rij-EPSILON STANDARD (LRR)
97
 
!   CSSG*  : CONSTANTES SPECIFIQUES DU RIJ-EPSILON SSG
98
 
!   CV2F*  : CONSTANTES SPECIFIQUES DU V2F PHI-MODEL
99
 
!   CKW*   : CONSTANTES SPECIFIQUES DU K-OMEGA SST
100
 
!            (SK=SIGMA_K, SW=SIGMA_W, BT=BETA, GM=GAMMA)
101
 
!   ALMAX  : ECHELLE DE LONGUEUR TURBULENTE
102
 
!   UREF   : VITESSE DE REFERENCE
103
 
!   XLOMLG : LONGUEUR POUR LONGUEUR DE MELANGE
104
 
!   XLESFL, ALES, BLES
105
 
!       DELTA = XLESFL * (ALES*VOLUME)^BLES (largeur du filtre utilise
106
 
!       en fonction du volume de la cellule)
107
 
!   CSMAGO
108
 
!       La constante de Smagorinsky theorique vaut 0.18
109
 
!       pour un canal plan, on prendra cependant plutot 0.065
110
 
!   XLESFD
111
 
!       Dans le cas d un modele dynamique, XLESFD est le rapport entre la
112
 
!       largeur du filtre explicite et celle du filtre implicite
113
 
!   SMAGMX
114
 
!       Constante de Smagorinsky maximale souhaitee (on peut prendre 10*CSMAGO)
115
 
!   IDRIES
116
 
!       Amortissement Van Driest active (=1) ou non (=0)
117
 
!   CDRIES
118
 
!       Constante de Van Driest dans (1-exp(-y+/CDRIES))
119
 
!   CE4    : Coefficient du terme interfacial dans k-eps
120
 
!            (Ce coefficient sert en Lagrangien)
121
 
!   VOLMIN : VOLUME DE CONTROLE MINIMAL
122
 
!   VOLMAX : VOLUME DE CONTROLE MAXIMAL
123
 
!   VOLTOT : VOLUME TOTAL DU DOMAINE
124
 
 
125
 
double precision  xkappa , cstlog , ypluli(nphsmx)  ,             &
126
 
                  apow   , bpow   , cpow   , dpow   ,             &
127
 
                  cmu    , cmu025 , ce1    , ce2    , ce4    ,    &
128
 
                  sigmak , sigmae ,                               &
129
 
                  crij1  , crij2  , crij3  , crijep , csrij  ,    &
130
 
                  crijp1 , crijp2 ,                               &
131
 
                  cssge2 , cssgs1 , cssgs2 ,                      &
132
 
                  cssgr1 , cssgr2 , cssgr3 , cssgr4 , cssgr5 ,    &
133
 
                  cv2fa1 , cv2fe2 , cv2fmu , cv2fc1 , cv2fc2 ,    &
134
 
                  cv2fct , cv2fcl , cv2fet ,                      &
135
 
                  ckwsk1 , ckwsk2 , ckwsw1 , ckwsw2 , ckwbt1 ,    &
136
 
                  ckwbt2 , ckwgm1 , ckwgm2 , ckwa1  , ckwc1  ,    &
137
 
                  volmin , volmax , voltot ,                      &
138
 
                  almax (nphsmx)  , uref  (nphsmx),               &
139
 
                  xlomlg(nphsmx)  ,                               &
140
 
                  xlesfl(nphsmx)  , ales  (nphsmx), bles(nphsmx), &
141
 
                  csmago(nphsmx)  , cdries(nphsmx),               &
142
 
                  xlesfd(nphsmx)  , smagmx(nphsmx),               &
143
 
                  cwale(nphsmx)
144
 
common / rturbu / xkappa , cstlog , ypluli ,                      &
145
 
                  apow   , bpow   , cpow   , dpow   ,             &
146
 
                  cmu    , cmu025 , ce1    , ce2    , ce4    ,    &
147
 
                  sigmak , sigmae ,                               &
148
 
                  crij1  , crij2  , crij3  , crijep , csrij  ,    &
149
 
                  crijp1 , crijp2 ,                               &
150
 
                  cssge2 , cssgs1 , cssgs2 ,                      &
151
 
                  cssgr1 , cssgr2 , cssgr3 , cssgr4 , cssgr5 ,    &
152
 
                  cv2fa1 , cv2fe2 , cv2fmu , cv2fc1 , cv2fc2 ,    &
153
 
                  cv2fct , cv2fcl , cv2fet ,                      &
154
 
                  ckwsk1 , ckwsk2 , ckwsw1 , ckwsw2 , ckwbt1 ,    &
155
 
                  ckwbt2 , ckwgm1 , ckwgm2 , ckwa1  , ckwc1  ,    &
156
 
                  volmin , volmax , voltot ,                      &
157
 
                  almax           , uref            ,             &
158
 
                  xlomlg          ,                               &
159
 
                  xlesfl          , ales            , bles      , &
160
 
                  csmago          , cdries          ,             &
161
 
                  xlesfd          , smagmx          ,             &
162
 
                  cwale
163
 
 
164
 
! CONSTANTES POUR LES SCALAIRES
165
 
 
166
 
! ISCSTH :
167
 
!   -1 : DE TYPE TEMPERATURE EN C (      CP POUR LA LOI DE PAROI)
168
 
!    0 : SCALAIRE PASSIF      (IE PAS DE CP POUR LA LOI DE PAROI)
169
 
!    1 : DE TYPE TEMPERATURE EN K (      CP POUR LA LOI DE PAROI)
170
 
!    2 : ENTHALPIE            (IE PAS DE CP POUR LA LOI DE PAROI)
171
 
!      LA DISTINCTION C/K SERT EN RAYONNEMENT
172
 
! IVISLS : SI POSITIF STRICTEMENT, INDIQUE QUE LA VISCOSITE ASSOCIEE
173
 
!            AU SCALAIRE EST VARIABLE, ET LA VALEUR EST LE NUMERO
174
 
!            D'ORDRE DE LA VISCOSITE DANS LE TABLEAU DES VISCOSITES
175
 
!            VARIABLES
176
 
! IVISSA : COMME IVISLS SAUF QUE SERT AU STOCKAGE DE LA VISCOSITE AU
177
 
!          PAS DE TEMPS PRECEDENT
178
 
! ICLVFL : 0 : CLIPPING DES VARIANCES A ZERO
179
 
!          1 : CLIPPING DES VARIANCES A ZERO ET A f(1-f)
180
 
!          2 : CLIPPING DES VARIANCES A MAX(ZERO,SCAMIN) ET SCAMAX
181
 
! ISCAVR : NUMERO DU SCALAIRE ASSOCIE A LA VARIANCE OU ZERO
182
 
!          SI LE SCALAIRE N'EST PAS UNE VARIANCE
183
 
! IPHSCA : NUMERO DE LA PHASE PORTEUSE
184
 
! SCAMIN, SCAMAX : MIN ET MAX POUR CLIPPING DES SCALAIRES
185
 
!                  ON NE CLIPPE QUE SI SCAMIN < SCAMAX
186
 
! VISLS0 : VISCOSITE DES SCALAIRES SI CONSTANTE
187
 
! SIGMAS : PRANDTL DES SCALAIRES
188
 
! RVARFL : COEFF DE DISSIPATION DES VARIANCES
189
 
 
190
 
integer           iscsth(nscamx),ivisls(nscamx),ivissa(nscamx),   &
191
 
                  iclvfl(nscamx),                                 &
192
 
                  iscavr(nscamx),iphsca(nscamx)
193
 
double precision  scamin(nscamx),scamax(nscamx),                  &
194
 
                  visls0(nscamx),sigmas(nscamx),                  &
195
 
                  rvarfl(nscamx)
196
 
common / iscala / iscsth        ,ivisls        ,ivissa        ,   &
197
 
                  iclvfl        ,                                 &
198
 
                  iscavr        ,iphsca
199
 
common / rscala / scamin        ,scamax        ,                  &
200
 
                  visls0        ,sigmas        ,                  &
201
 
                  rvarfl
202
 
 
203
 
! FIN