1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
|
!-------------------------------------------------------------------------------
! This file is part of Code_Saturne, a general-purpose CFD tool.
!
! Copyright (C) 1998-2016 EDF S.A.
!
! This program is free software; you can redistribute it and/or modify it under
! the terms of the GNU General Public License as published by the Free Software
! Foundation; either version 2 of the License, or (at your option) any later
! version.
!
! This program is distributed in the hope that it will be useful, but WITHOUT
! ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
! FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
! details.
!
! You should have received a copy of the GNU General Public License along with
! this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
! Street, Fifth Floor, Boston, MA 02110-1301, USA.
!-------------------------------------------------------------------------------
subroutine cs_coal_param
!=======================
!===============================================================================
! FONCTION :
! ---------
! INIT DES OPTIONS DES VARIABLES TRANSPORTEES
! ET DES VARIABLES ALGEBRIQUES
! COMBUSTION CHARBON PULVERISE
!-------------------------------------------------------------------------------
! Arguments
!__________________.____._____.________________________________________________.
! name !type!mode ! role !
!__________________!____!_____!________________________________________________!
!__________________!____!_____!________________________________________________!
! TYPE : E (ENTIER), R (REEL), A (ALPHANUMERIQUE), T (TABLEAU)
! L (LOGIQUE) .. ET TYPES COMPOSES (EX : TR TABLEAU REEL)
! MODE : <-- donnee, --> resultat, <-> Donnee modifiee
! --- tableau de travail
!===============================================================================
!===============================================================================
! Module files
!===============================================================================
use paramx
use dimens
use numvar
use optcal
use cstphy
use entsor
use cstnum
use ppppar
use ppthch
use coincl
use cpincl
use ppincl
use ihmpre
use ppcpfu
use cs_coal_incl
use field
!===============================================================================
implicit none
integer ii , jj , iok
integer icha , isc
double precision wmolme
!===============================================================================
! 1. VARIABLES TRANSPORTEES
!===============================================================================
! --> Nature des scalaires transportes
do isc = 1, nscapp
! ---- Type de scalaire (0 passif, 1 temperature en K
! -1 temperature en C
! 2 enthalpie)
! La distinction -1/1 sert pour le rayonnement
iscacp(iscapp(isc)) = 0
enddo
! ---- On resout en enthalpie avec un CP constant (Cf. cpvarp)
itherm = 2
iscacp(iscalt) = 0
! --> Donnees physiques ou numeriques propres aux scalaires CP
do isc = 1, nscapp
jj = iscapp(isc)
if ( iscavr(jj).le.0 ) then
! En combustion on considere que la viscosite turbulente domine
! ON S'INTERDIT DONC LE CALCUL DES FLAMMES LAMINAIRES AVEC Le =/= 1
visls0(jj) = viscl0
endif
! ------ Schmidt ou Prandtl turbulent
sigmas(jj) = 0.7d0
! ------ Coeff dissipation des fluctuations
rvarfl(jj) = 0.8d0
ii = isca(iscapp(isc))
! ------ Niveau de detail des impressions pour les variables et
! donc les scalaires (valeurs 0 ou 1)
! Si = -10000 non modifie par l'utilisateur -> niveau 1
if (iwarni(ii).eq.-10000) then
iwarni(ii) = 1
endif
! - Interface Code_Saturne:
! ======================
! ihisvr is already filled in UINUM1 routine
if (iihmpr.ne.1) then
! ------ Informations relatives a la resolution des scalaires
! - Facteur multiplicatif du pas de temps
cdtvar(ii) = 1.d0
! - Schema convectif % schema 2ieme ordre
! = 0 : upwind
! = 1 : second ordre
blencv(ii) = 0.d0
! - Type de schema convetif second ordre (utile si BLENCV > 0)
! = 0 : Second Order Linear Upwind
! = 1 : Centre
ischcv(ii) = 1
! - Test de pente pour basculer d'un schema centre vers l'upwind
! = 0 : utlisation automatique du test de pente
! = 1 : calcul sans test de pente
isstpc(ii) = 0
! - Reconstruction des flux de convetion et de diffusion aux faces
! = 0 : pas de reconstruction
ircflu(ii) = 0
endif
enddo
!===============================================================================
! 2. INFORMATIONS COMPLEMENTAIRES
!===============================================================================
! ---> Initialisation
! ---- Calcul de RO0 a partir de T0 et P0
! (loi des gaz parfaits applliquee a l'air)
! On initialise RO0 avec l'oxydant 1 qui est sense etre
! l'oxydant majoritaire
wmolme = ( wmole(io2) *oxyo2(1) +wmole(in2) *oxyn2(1) &
+wmole(ih2o)*oxyh2o(1)+wmole(ico2)*oxyco2(1)) &
/(oxyo2(1)+oxyn2(1)+oxyh2o(1)+oxyco2(1))
ro0 = p0*wmolme / (cs_physical_constants_r*t0)
! ---- Initialisation pour la masse volumique du coke
do icha = 1, ncharb
rhock(icha) = rho0ch(icha)
enddo
! On met SRROM et DIFTL0 a -GRAND pour forcer l'utilisateur a les
! definir dans cs_user_combustion
! ---> Coefficient de relaxation de la masse volumique
! RHO(n+1) = SRROM * RHO(n) + (1-SRROM) * RHO(n+1)
srrom = -grand
! ---> Viscosite laminaire associee au scalaire enthalpie
! DIFTL0 (diffusivite dynamique en kg/(m s))
! C'est cette valeur par defaut qui est TOUJOURS utilisee dans les
! calculs charbon (un peu etonnant de ne pas prendre en
! compte les variations de ce parametre physique si on
! recherche des informations sur les flux thermiques aux parois)
diftl0 =-grand
! ---> Masse volumique variable et viscosite constante (pour les suites)
irovar = 1
ivivar = 0
!===============================================================================
! 3. ON REDONNE LA MAIN A L'UTLISATEUR
!===============================================================================
! - Interface Code_Saturne
! ======================
if (iihmpr.eq.1) then
call uicpi1(srrom, diftl0)
!==========
diftl0 = 4.25d-5
endif
call cs_user_combustion
!===============================================================================
! 4. VERIFICATION DES DONNERS FOURNIES PAR L'UTLISATEUR
!===============================================================================
iok = 0
call cs_coal_verify (iok)
!=====================
if (iok.gt.0) then
write(nfecra,9999)iok
call csexit (1)
!==========
else
write(nfecra,9998)
endif
!--------
! Formats
!--------
9998 format( &
' ',/,&
' Pas d erreur detectee lors de la verification des donnees ',/,&
' (cs_user_combustion).',/)
9999 format( &
'@ ',/,&
'@ ',/,&
'@ ',/,&
'@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,&
'@ ',/,&
'@ @@ ATTENTION : ARRET A L''ENTREE DES DONNEES ',/,&
'@ ========= ',/,&
'@ LES PARAMETRES DE CALCUL SONT INCOHERENTS OU INCOMPLETS ',/,&
'@ ',/,&
'@ Le calcul ne sera pas execute (',I10,' erreurs). ',/,&
'@ ',/,&
'@ Se reporter aux impressions precedentes pour plus de ',/,&
'@ renseignements. ',/,&
'@ Verifier cs_user_combustion.' ,/,&
'@ ',/,&
'@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@',/,&
'@ ',/)
!----
! End
!----
return
end subroutine
|
