Télécharger konv.notice

Retour à la liste thématique

Afficher cette notice en

Numérotation des lignes :
   1 : $$$$ KONV     NOTICE  CHAT      11/09/12    21:16:46     7124           
   2 :                                              DATE     11/09/12
   3 :          
   4 :   Operateur KONV                          Voir aussi : NAVI
   5 :      --------------  
   6 : 
   7 : 
   8 :      OBJET   :
   9 : 
  10 :  I Formulation Volume Finis (OPTI VF) :
  11 :  ______________________________________
  12 : 
  13 :       Discretise l'operateur de convection par des schemas volumes finis 
  14 :     Il calcule un champ point FACE  qui represente
  15 :     le flux de chaque arete du maillage selon 3 schemas.
  16 :     On ecrit le flux dans l'increment KIZG.
  17 : 
  18 :      SYNTAXE 1 :  KONV VNF VNC option;
  19 :                   VNF = U calcule aux sommets (par ksof)
  20 :                   VNC = U calcule aux centres (par knol)
  21 :                   Option = 'upwind' | 'quick' | 'muscl'
  22 : 
  23 :      SYNTAXE 2 :  en utilisant EQEX
  24 :                   RV = EQEX TABDOM ALFA .. ITMA ..
  25 :                        'ZONE' MOD1 'OPER' 'KONV' VNF VNC 'INCO' 'TN'
  26 :                         CLIM 'TN' TIMP ENTREE ... ;
  27 : 
  28 :            on met l'option dans klop : klop ind rv 'option' ;
  29 : 
  30 :    Erreur possible dans le calcul de DT si VN est petit changer EPSILON
  31 : 
  32 : 
  33 :       COMMENTAIRES :
  34 : 
  35 : 
  36 : Le flux est evalue sur chaque maille par : FI = VNF TK LGR
  37 : ou TK est la temperature sur la face K calculee de facons differentes selon
  38 : les schemas a partir des temperatures sur les mailles connexes: TIM ,TI ET TIP.
  39 : (Voir le rapport pour le choix de ces trois temperatures)
  40 : 
  41 : 
  42 : 
  43 : -Schema UPWIND :
  44 : 
  45 :         TK = TI
  46 : 
  47 : 
  48 : -Schema QUICK : 
  49 :   
  50 :         TK = (TI+TIP)/2 + (TIM+TIP-2*TI)/8
  51 : 
  52 : 
  53 : -Schema MUSCL : 
  54 : 
  55 : 
  56 :         TK = TI + DELTA_I/2
  57 : 
  58 : 
  59 : ou DELTA_I est la pente entre deux mailles :
  60 : 
  61 :         GAMMA_I = SIGNE (TIP-TI)
  62 :         DELTA_I = GAMMA_I * MAX (0,MIN(GAMMA_I*(TIP-TI),GAMMA_I*(TI-TIM)))
  63 : 
  64 : 
  65 : 
  66 :       
  67 : 
  68 : Le bilan des flux precedents est fait dans l'operateur AVCT
  69 : 
  70 : 
  71 :  II Formulation Element Finis (OPTI EF ou EFM1) :
  72 :  ________________________________________________
  73 : 
  74 :       Discretise l'operateur de convection en Element Finis.
  75 :    Suivant l'option l'operateur est traites sous forme
  76 :   conservative ou non conservative.
  77 : 
  78 :     SYNTAXE - EQEX  Cf operateur EQEX
  79 :     ______________
  80 : 
  81 :     'OPER' 'KONV' ROC UN  LAM           'INCO' 'TN' :
  82 : 
  83 : 
  84 :      1/ Formulation non conservative
  85 : 
  86 :          roc ( u Grad T )
  87 : 
  88 :      2/ Formulation conservative
  89 : 
  90 :           Div ( roc u T )
  91 : 
  92 : 
  93 :     Commentaires :
  94 :     ______________
  95 : 
  96 :      roc  capacite calorique (J/M**3/°C)
  97 :           FLOTTANT ou CHPOINT SCAL CENTRE ou CHPOINT SCAL SOMMET ou MOT
  98 :      lam  conductivite thermique (W/M/°C)
  99 :           cette donnee est necessaire pour l'evaluation
 100 :           du Peclet de maille. En general c'est le coefficent
 101 :           de l'operateur LAPN. Si c'est operateur est absent
 102 :           mettre lam=0.
 103 :           FLOTTANT ou CHPOINT SCAL CENTRE ou CHPOINT SCAL SOMMET ou MOT
 104 :      un   Champ de vitesse transportant
 105 :           CHPOINT VECT SOMMET ou MOT
 106 :      tn   Champ de temperature
 107 :           CHPOINT SCAL SOMMET ou MOT
 108 : 
 109 :  Un coefficient de type MOT indique que l'operateur va chercher le
 110 :  coefficient dans la table INCO a l'indice MOT.
 111 : 
 112 : 
 113 :     Options : (EQEX)
 114 :     _________
 115 : 
 116 :  La discretisation des termes de convection peut etre :
 117 : 
 118 :  centree                              OPTION CENTREE
 119 :  decentree                            OPTION SUPG
 120 :  decentree avec capture de choc       OPTION SUPGCC   Option par defaut
 121 :  Crank Nicholson generalise           OPTION CNG
 122 :  (ordre 4 en temps)
 123 : 
 124 : 
 125 :  Formulation non conservative         OPTION NOCONS   Option par defaut
 126 :  Formulation conservative             OPTION CONS
 127 : 
 128 :  Formulation EF                       OPTION EF
 129 : 
 130 : 
 131 : 
 132 :  III Discretisation des Equations d'Euler 
 133 :  _________________________________________
 134 : 
 135 :  IIIa : gaz parfait mono-constituent polytropique
 136 :  ________________________________________________
 137 : 
 138 : 
 139 :  Discretisation en VF "cell-centered" des equations d'Euler pour un gaz 
 140 :  parfait mono-constituent polytropique
 141 : 
 142 :  Inconnues:
 143 : 
 144 :  densite, quantite de mouvement (qdm), energie totale par unite de volume 
 145 :  (variables conservatives)
 146 : 
 147 :  ou
 148 : 
 149 :  densite, vitesse, pression (variables primitives)
 150 : 
 151 :  On peut calculer:
 152 : 
 153 :  IIIa.1. Le flux numerique
 154 :  IIIa.2  Le residu
 155 :  IIIa.3  La matrice jacobienne du residu par rapport aux variables  
 156 :          conservatives
 157 :  IIIa.4  La matrice jacobienne du residu par rapport aux variables  
 158 :          primitives
 159 :  IIIa.5  La matrice de preconditionnent bas Mach par rapport aux 
 160 :          variables conservatives
 161 :  IIIa.6  La matrice de preconditionnent bas Mach par rapport aux
 162 :          variables primitives
 163 :  IIIa.7  La contribution de quelque condition limite au residu et
 164 :          a la matrice jacobienne
 165 : 
 166 : 
 167 :  IIIa.1 et IIIa.2  Le flux numerique et le residu
 168 :  ________________
 169 : 
 170 : 
 171 :  RCHPO1 RFLOT1 = 'KONV' 'VF' 'PERFMONO' MOT1 MOT2  MOD1 LMOT1 MCHAM1 
 172 :                  MCHAM2 MCHAM3 MCHAM4 (CHPO5 CHPO6) (MAIL1) ;
 173 : 
 174 :  ENTRÉES 
 175 : 
 176 :  MOT1   : objet de type MOT
 177 :            Il vaut 'RESI' si on veut calculer le residu
 178 :            Il vaut 'FLUX' si on veut calculer le flux
 179 : 
 180 :  MOT2   : objet de type MOT
 181 :            Il indique la methode:
 182 :            'GODUNOV'  = solveur exacte
 183 :            'VANLEER'  = solveur de van Leer
 184 :            'VLH'      = solveur de van Leer Hanel
 185 :            'HUSVL'    = HUS (van Leer + Osher) 
 186 :            'HUSVLH'   = HUS (van Leer Hanel + Osher) 
 187 :            'AUSMPLUS' = AUSM+ 
 188 :            'ROE'      = solveur de Roe
 189 :            'SS'       = solveur choc-choc
 190 :            'AUSMPLM'  = AUSM+ low Mach
 191 :            'RUSANOV'  = solveur de Rusanov 
 192 :            'RUSANOLM' = solveur de Rusanov pour le bas-Mach
 193 :            'CENTERED' = solveur centre
 194 :            'ROELM'    = solveur de Roe-Turkel pour le bas-Mach
 195 :            'HLLC'     = solveur HLLC
 196 :            'HLLCLM'   = solveur HLLC-Turkel pour le bas-Mach
 197 :            'AUSMPUP'  = solveur AUSM+up low Mach
 198 : 
 199 :   
 200 :  LMOT1   : objet de type LISTMOTS
 201 :            Noms de composantes du resultat (RCHPO1)
 202 :            Il contient dans l'ordre suivant: le noms de la densite,
 203 :            de la qdm, de l'energie totale par unite de volume
 204 : 
 205 :  MOD1    : objet MODELE.
 206 : 
 207 :  MCHAM1  : MCHAML contenant la masse volumique, qui a comme 
 208 :            SPG (support geometrique) 'DOMA' MOD1 'FACEL'
 209 :            (une composante, 'SCAL')
 210 :            (voir operateur PRET)
 211 : 
 212 :  MCHAM2  : MCHAML  contenant la vitesse  et les cosinus 
 213 :            directeurs du repere locale (n,t) dans le repere 
 214 :            global (x,y) (dans le cas 2D 6 composantes:
 215 :            * 'UN' = vitesse normale  (SPG = 'DOMA' MOD1 'FACEL')
 216 :            * 'UT' = vitesse tangentielle  (SPG = 'DOMA' MOD1 FACEL')
 217 :            * 'NX' = n.x (SPG = 'FACE')
 218 :            * 'NY' = n.y (SPG = 'FACE')
 219 :            * 'TX' = t.x (SPG = 'FACE')
 220 :            * 'TY' = t.y (SPG = 'FACE')).
 221 :            (voir operateur PRET)
 222 : 
 223 :  MCHAM3  : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la pression du
 224 :            gaz (une  seule composante, 'SCAL').
 225 :            (voir operateur PRET)
 226 : 
 227 :  MCHAM4  : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant le "gamma" du 
 228 :            gaz (une seule composante, 'SCAL').
 229 :            (voir operateur PRET)
 230 : 
 231 :  CHPO5   : CHPOINT contenant la vitesse de cut-off 
 232 :            (SPG = 'DOMA' MOD1 'CENTRE', une seule composante, 
 233 :             'SCAL'). 
 234 :            A donner dans le cas bas Mach
 235 : 
 236 :  CHPO6   : CHPOINT contenant la deuxieme vitesse de cut-off 
 237 :            (SPG = 'DOMA' MOD1 'CENTRE', une seule composante, 
 238 :             'SCAL'). 
 239 :            A donner dans le cas bas Mach
 240 : 
 241 : (MAIL1)  : MAILLAGE de POI1 ou, si appartenant a  ('DOMA' MOD1 'FACE'),
 242 :            on ne calcule pas la contribution au flux ou au residu.
 243 : 
 244 :  SORTIES
 245 : 
 246 :  RCHPO1  : objet de type CHPOINT (composantes =  LMOT1)
 247 :            Residu  si MOT2 = 'RESI' (SPG = 'DOMA' MOD1 'CENTRE')
 248 :            Flux si MOT2 = 'FLUX'    (SPG = 'DOMA' MOD1 'FACE') 
 249 : 
 250 :  RFLOT1  : objet de type FLOTTANT
 251 :            Il est le temps caracteristique associe a l'onde la plus
 252 :            rapide (meme dans le cas bas Mach, on considere le systeme 
 253 :            non preconditionne)
 254 : 
 255 :  Remarque
 256 :  --------
 257 : 
 258 :  RCHPO1 est egal a:
 259 :  * la derive temporelle des inconnues si l'option 'RESI' est utilisee
 260 :  * la projection du flux  sur ('DOMA' MOD1 'XXNORMAF') si l'option 
 261 :    'FLUX' est utilisee
 262 : 
 263 : 
 264 :  IIIa.3   La matrice jacobienne du residu par rapport aux variables  
 265 :  _____    conservatives
 266 : 
 267 : 
 268 :  RMAT1  = 'KONV' 'VF' 'PERFMONO' 'JACOCONS' MOD1 LMOT1 (MAIL1) MOT1 
 269 :           CHPO1 CHPO2 CHPO3 CHPO4 (CHPO5 CHPO6) ;
 270 :  
 271 :  ENTRÉES 
 272 : 
 273 :  MOT1   : objet de type MOT
 274 :           'VLH'       : jacobien du residu pour la methode VLH 
 275 :           'AUSMPLUS'  : jacobien du residu pour la methode AUSM+ 
 276 :           'AUSMPLM'   : jacobien du residu pour la methode AUSM+ low Mach
 277 :           'RUSANOLM'  : jacobien du residu pour la methode de Rusanov low Mach
 278 :           'CENTERED'  : jacobien du residu pour la methode centree
 279 : 
 280 :  LMOT1   : objet de type LISTMOTS
 281 :            Noms des variables conservatives
 282 :            Il contient dans l'ordre suivant: le noms de la densite,
 283 :            de la qdm, de l'energie totale par unite de volume
 284 : 
 285 :  MOD1    : objet MODELE.
 286 : 
 287 :  CHPO1   : CHPOINT contenant la masse volumique
 288 :            (SPG = 'DOMA' DOM1 'CENTRE', une seule composante, 
 289 :             'SCAL').
 290 : 
 291 :  CHPO2   : CHPOINT contenant la vitesse
 292 :            (SPG = 'DOMA' DOM1 'CENTRE', deux/trois composantes
 293 :             'UX', 'UY', 'UZ')
 294 : 
 295 :  CHPO3   : CHPOINT contenant la pression du gaz
 296 :            (SPG = 'DOMA' DOM1 'CENTRE', une seule composante, 
 297 :             'SCAL').
 298 : 
 299 :  CHPO4   : CHPOINT contenant le "gamma" du gaz
 300 :            (SPG = 'DOMA' DOM1 'CENTRE', une seule composante, 
 301 :             'SCAL').
 302 : 
 303 :  CHPO5   : CHPOINT contenant la vitesse de cut-off
 304 :            (SPG = 'DOMA' DOM1 'CENTRE', une seule composante, 
 305 :             'SCAL'). 
 306 :            A donner dans le cas bas Mach
 307 : 
 308 :  CHPO6   : CHPOINT contenant la deuxieme vitesse de cut-off 
 309 :            (SPG = 'DOMA' DOM1 'CENTRE', une seule composante, 
 310 :             'SCAL'). 
 311 :            A donner dans le cas bas Mach
 312 : 
 313 :  (MAIL1) : MAILLAGE de POI1 ou, si appartenant a  ('DOMA' DOM1 'FACE'),
 314 :            on ne calcule pas la contribution au jacobien.
 315 : 
 316 :  SORTIES
 317 : 
 318 :  RMAT1   : objet de type MATRIK
 319 :            (SPG =  'DOMA' DOM1 'CENTRE')
 320 :            (inconnues primales = inconnues duales = LMOT1)
 321 :            Il contient le jacobien du residu par rapport aux variables 
 322 :            conservatives.
 323 : 
 324 : 
 325 :  IIIa.4   La matrice jacobienne du residu par rapport aux variables  
 326 :  ______   primitives
 327 : 
 328 :  RMAT1  = 'KONV' 'VF' 'PERFMONO' 'JACOPRIM' MOD1 LMOT1 LMOT2 (MAIL1) MOT1 
 329 :           CHPO1 CHPO2 CHPO3 CHPO4 (CHPO5 CHPO6) ;
 330 :  
 331 :  ENTREES 
 332 : 
 333 :  MOT1   : objet de type MOT
 334 :           'VLH'       : jacobien du residu pour la methode VLH 
 335 :           'AUSMPLUS'  : jacobien du residu pour la methode AUSM+ 
 336 :           'AUSMPLM'   : jacobien du residu pour la methode AUSM+ low Mach
 337 : 
 338 :  LMOT1   : objet de type LISTMOTS
 339 :            Noms des variables conservatives
 340 :            Il contient dans l'ordre suivant: le noms de la densite,
 341 :            de la qdm, de l'energie totale par unite de volume
 342 : 
 343 :  LMOT2   : objet de type LISTMOTS
 344 :            Noms des variables primitives
 345 :            Il contient dans l'ordre suivant: le noms de la densite,
 346 :            de la vitesse, de la pression
 347 : 
 348 :  MOD1    : objet MODELE.
 349 : 
 350 :  CHPO1   : CHPOINT contenant la masse volumique
 351 :            (SPG = 'DOMA' MOD1 'CENTRE', une seule composante, 
 352 :             'SCAL').
 353 : 
 354 :  CHPO2   : CHPOINT contenant la vitesse
 355 :            (SPG = 'DOMA' MOD1 'CENTRE', deux/trois composantes
 356 :             'UX', 'UY', 'UZ')
 357 : 
 358 :  CHPO3   : CHPOINT contenant la pression du gaz
 359 :            (SPG = 'DOMA' MOD1 'CENTRE', une seule composante, 
 360 :             'SCAL').
 361 : 
 362 :  CHPO4   : CHPOINT contenant le "gamma" du gaz
 363 :            (SPG = 'DOMA' MOD1 'CENTRE', une seule composante, 
 364 :             'SCAL').
 365 : 
 366 :  CHPO5   : CHPOINT contenant la vitesse de cut-off 
 367 :            (SPG = 'DOMA' MOD1 'CENTRE', une seule composante, 
 368 :             'SCAL'). 
 369 :            A donner dans le cas bas Mach
 370 : 
 371 :  CHPO6   : CHPOINT contenant la deuxieme vitesse de cut-off 
 372 :            (SPG = 'DOMA' MOD1 'CENTRE', une seule composante, 
 373 :             'SCAL'). 
 374 :            A donner dans le cas bas Mach
 375 : 
 376 :  (MAIL1) : MAILLAGE de POI1 ou, si appartenant a  ('DOMA' MOD1 'FACE'),
 377 :            on ne calcule pas la contribution au jacobien.
 378 : 
 379 :  SORTIES
 380 : 
 381 :  RMAT1   : objet de type MATRIK
 382 :            (SPG = 'DOMA' MOD1 'CENTRE')
 383 :            (inconnues primales = variables primitives  = LMOT2)
 384 :            (inconnues duales = variables conservatives = LMOT1)
 385 :            Il contient le jacobien du residu par rapport aux variables 
 386 :            primitives.
 387 : 
 388 : 
 389 :  IIIa.5   La matrice de preconditionnement bas Mach par rapport aux 
 390 :  ______   variables conservatives (divise par le pas de temps locale;
 391 :           le pas de temps est calcule pour le systeme preconditionnne
 392 :           (intervalle de temps pendant lequel le signale le plus rapide
 393 :            traverse une cellule reguliere))
 394 : 
 395 :  RMAT1 =  'KONV' 'VF' 'PERFMONO' 'GAMMCONS' MAIL1 LMOT1
 396 :            CHPO1 CHPO2 CHPO3 CHPO4 CHPO5 CHPO6 ;
 397 : 
 398 :  LMOT1   : objet de type LISTMOTS
 399 :            Noms des variables conservatives
 400 :            Il contient dans l'ordre suivant: le noms de la densite,
 401 :            de la qdm, de l'energie totale par unite de volume
 402 : 
 403 :  MAIL1   : SPG des CHPOINTs
 404 : 
 405 :  CHPO0   : CHPOINT contenant le diametre des elts
 406 :            (une seule composante, 'SCAL').
 407 : 
 408 :  CHPO1   : CHPOINT contenant la masse volumique
 409 :            (une seule composante, 'SCAL').
 410 : 
 411 :  CHPO2   : CHPOINT contenant la vitesse
 412 :            (deux/trois composantes
 413 :             'UX', 'UY', 'UZ')
 414 : 
 415 :  CHPO3   : CHPOINT contenant la pression du gaz
 416 :            (une seule composante, 
 417 :             'SCAL').
 418 : 
 419 :  CHPO4   : CHPOINT contenant le "gamma" du gaz
 420 :            (une seule composante, 
 421 :             'SCAL').
 422 : 
 423 :  CHPO5   : CHPOINT contenant la vitesse de cut-off 
 424 :            (SPG = 'CENTRE', une seule composante, 
 425 :             'SCAL'). 
 426 :             A donner dans le cas bas Mach
 427 : 
 428 :  CHPO6   : CHPOINT contenant la deuxieme vitesse de cut-off 
 429 :            (SPG = 'CENTRE', une seule composante, 
 430 :             'SCAL'). 
 431 :              A donner dans le cas bas Mach
 432 : 
 433 :  RMAT1   : objet de type MATRIK
 434 :            (SPG =  MAIL1)
 435 :            (inconnues primales = variables duales  = LMOT1)
 436 : 
 437 :  IIIa.6   La matrice de preconditionnement bas Mach par rapport aux 
 438 :  ______   variables primitives (divise par le pas de temps locale;
 439 :           le pas de temps est calcule pour le systeme preconditionnne
 440 :           (intervalle de temps pendant lequel le signale le plus rapide
 441 :            traverse une cellule reguliere))
 442 : 
 443 :  RMAT1 =  'KONV' 'VF' 'PERFMONO' 'GAMMPRIM' MAIL1 LMOT1 LMOT2
 444 :            CHPO1 CHPO2 CHPO3 CHPO4 CHPO5 CHPO6 ;
 445 : 
 446 :  LMOT1   : objet de type LISTMOTS
 447 :            Noms des variables conservatives
 448 :            Il contient dans l'ordre suivant: le noms de la densite,
 449 :            de la qdm, de l'energie totale par unite de volume
 450 : 
 451 :  LMOT2   : objet de type LISTMOTS
 452 :            Noms des variables primitives
 453 :            Il contient dans l'ordre suivant: le noms de la densite,
 454 :            de la vitesse, de la pression
 455 : 
 456 :  MAIL1   : SPG des CHPOINTs
 457 : 
 458 :  CHPO0   : CHPOINT contenant le diametre des elts
 459 :            (une seule composante, 'SCAL').
 460 : 
 461 :  CHPO1   : CHPOINT contenant la masse volumique
 462 :            (une seule composante, 'SCAL').
 463 : 
 464 :  CHPO2   : CHPOINT contenant la vitesse
 465 :            (deux/trois composantes
 466 :             'UX', 'UY', 'UZ')
 467 : 
 468 :  CHPO3   : CHPOINT contenant la pression du gaz
 469 :            (une seule composante, 
 470 :             'SCAL').
 471 : 
 472 :  CHPO4   : CHPOINT contenant le "gamma" du gaz
 473 :            (une seule composante, 
 474 :             'SCAL').
 475 : 
 476 :  CHPO5   : CHPOINT contenant la vitesse de cut-off 
 477 :            (SPG = 'CENTRE', une seule composante, 
 478 :             'SCAL'). 
 479 :            A donner dans le cas bas Mach
 480 : 
 481 :  CHPO6   : CHPOINT contenant la deuxieme vitesse de cut-off 
 482 :            (SPG = 'CENTRE', une seule composante, 
 483 :             'SCAL'). 
 484 :            A donner dans le cas bas Mach
 485 : 
 486 :  RMAT1   : objet de type MATRIK
 487 :            (SPG = 'CENTRE')
 488 :            (inconnues primales = variables primitives  = LMOT2)
 489 :            (inconnues duales = variables conservatives = LMOT1)
 490 : 
 491 : 
 492 :  IIIa.7   La contribution de quelque condition limite au residu et
 493 :  ------  a la matrice jacobienne
 494 : 
 495 :  RCHPLI RCHPRE = 'KONV' 'VF' 'PERFMONO' 'CLIM'
 496 :                  'RESI' $MOD1 $MOD2 LMOTC LMOTP CHPRN CHPVN CHPPN 
 497 :                   CHPGN  MOT1 CHPLI ;
 498 : 
 499 :  ou 
 500 : 
 501 :  RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'JACOCONS' 
 502 :          $MOD1 $MOD2 LMOTC LMOTP CHPRN CHPVN CHPPN 
 503 :                   CHPGN  MOT1 CHPLI ;
 504 : 
 505 :  RJACO = 'KONV' 'VF' 'PERFMONO' 'CLIM'  'JACOPRIM' 
 506 :          $MOD1 $MOD2 LMOTC LMOTP CHPRN CHPVN CHPPN 
 507 :                   CHPGN  MOT1 CHPLI ;
 508 : 
 509 :  $MOD1   : l'objet modele du domaine total
 510 : 
 511 :  $MOD2   : l'objet modele du domaine du bord 
 512 : 
 513 :  LMOTC   : LISTMOTS, noms des variables conservatives
 514 : 
 515 :  LMOTP   : LISTMOTS, noms des variables primitives
 516 : 
 517 :  CHPRN   : densite  (SPG = 'DOMA' $MOD1 'CENTRE', une seule
 518 :            composante, 'SCAL')
 519 : 
 520 :  CHPVN   : vitesse (SPG = 'DOMA' $MOD1 'CENTRE',
 521 :            composantes: 'UX', 'UY', ('UZ'))
 522 : 
 523 :  CHPPN   : pression  (SPG = 'DOMA' $MOD1 'CENTRE', une seule 
 524 :            composante, 'SCAL')
 525 : 
 526 :  CHPGN   : gamma  (SPG = 'DOMA' $MOD1 'CENTRE', une seule 
 527 :            composante, 'SCAL')
 528 : 
 529 :  MOT1    : MOT, type de condition limite:
 530 :            'INRI' : entree subsonique. On utilise les invariants de
 531 :                     Riemann. On donne densite 'RN', vitesse 'UX' 'UY' 
 532 :                     ('UZ'), pression 'PN'.
 533 :            'INSU' : entree subsonique. On donne l'enthalpie totale 'HT'
 534 :                     par unite de masse, l'entropie 'S' (pression divisee 
 535 :                     par la densite a la puissance gamma). On impose que 
 536 :                     la vitesse tangentielle est nulle. On recupere la 
 537 :                     pression a l'interieur. La contribution au residu et 
 538 :                     a la matrice jacobienne sont calculee via la methode
 539 :                     'AUSMPLUS'.
 540 :            'OUTRI': sortie subsonique. On utilise les invariants de
 541 :                     Riemann. On donne densite 'RN', vitesse 'UX' 'UY' 
 542 :                     ('UZ'), pression 'PN'.
 543 :            'OUTP' : sortie subsonique. On donne la pression 'PN', 
 544 :                     on recupere la vitesse et la densite a l'interieur.
 545 :                     La contribution au residu et a la matrice jacobienne 
 546 :                     sont calculee via la methode 'AUSMPLUS'.
 547 :            'INSS' : entree supersonique. On donne densite 'RN', vitesse 
 548 :                     'UX' 'UY' ('UX'), pression 'PN'.
 549 :            'OUTSS': sortie supersonique. On donne un CHPOINT vide.
 550 :            'INJE' : condition limite d'injection compressible. 
 551 :                     On donne le flux de masse 'MOME' et 'RT' (temperature fois
 552 :                     la constante de gaz); on impose que la vitesse
 553 :                     tangentielle est nulle, on recupere la pression.
 554 :           'INJELM': condition limite d'injection faiblement compressible. 
 555 :                     On donne le flux de masse et RT (temperature fois
 556 :                     la constante de gaz); on impose que la vitesse
 557 :                     tangentielle est nulle, on recupere la pression.
 558 : 
 559 :  CHPLI   : condition limite impose (SPG ='DOMA' $MOD2 'CENTRE')
 560 :            Les composantes dependent de MOT1
 561 : 
 562 :  Resultats
 563 : 
 564 :  RCHPLI : CHPOINT qui contient les valeur de la densite, de la vitesse et 
 565 :           de pression associe a la condition limite.
 566 :           (SPG ='DOMA' $MOD2 'CENTRE', composantes = LMOTP)
 567 : 
 568 :  RCHPRE : CHPOINT qui contient la contribution de la condition limite au 
 569 :           residu (SPG en 'DOMA' $MOD1 'CENTRE', composantes = LMOTC)
 570 : 
 571 :  RJACO  : MATRIK qui contient la contribution de la condition limite a la
 572 :           matrice jacobienne du residu.
 573 :  
 574 : 
 575 :  IIIb :gaz thermiquement parfait
 576 :  ____
 577 : 
 578 :  Discretisation en VF "cell-centered" des equations d'Euler pour
 579 :  un melange de gaz parfaits.
 580 :  Inconnues: densite (totale) du melange, quantite de mouvement, 
 581 :  energie totale par unite de volume, densites des constituants du gaz, 
 582 :  scalaires passifs (multiplies par la densite du melange)
 583 :  
 584 :  RCHPO1 RFLOT1 = 'KONV' 'VF' 'PERFTEMP' MOT1 MOT2
 585 :                         MOD1 TAB2 LMOT1  MCHAM1 MCHAM2 MCHAM3
 586 :                         (MCHAM4) (MCHAM5) ;
 587 :  
 588 :  ENTREES 
 589 :  
 590 :  
 591 :  MOT1   : objet de type MOT
 592 :           Il vaut 'RESI' si on veut calculer le residu
 593 :           Il vaut 'FLUX' si on veut calculer le flux
 594 :  
 595 :  MOT2   : objet de type MOT
 596 :           Il indique la methode de decentrement:
 597 :           'VLH'     = solveur de van Leer Hanel
 598 :           'SS'      = solveur choc-choc
 599 :  
 600 :  MOD1  : objet MODELE
 601 : 
 602 :  TAB2  : table qui contient les proprietes du gaz. Plus precisement:
 603 :          * le nom de l'espece qui n'est pas dans les Equations d'Euler 
 604 :            en TAB2 . 'ESPNEULE' (MOT);
 605 :          * les noms des especes qui apparaissent explicitement dans les
 606 :            equations d'Euler en  TAB2 . 'ESPEULE' (LISTMOTS); dans le 
 607 :            cas mono-espece cet indice n'existe pas.
 608 :          * le degre des polynomes cv_i=cv_i(T), en TAB2 . 'NORD'
 609 :            (ENTIER >= 0)
 610 :          * les proprietes de chaque gaz 'ESPI', dans la table
 611 :            TAB2 . 'ESPI':
 612 :            - TAB2 . 'ESPI' . 'A' (LISTREEL)
 613 :              qui contient les (TAB2.'NORD')+1 coefficients des polynomes
 614 :              cv(T), (A0,A1,...); si le cv(T) sont supposes etre en 
 615 :              J/kg/K, les Ai doivent etre en unites coherentes; 
 616 :            - TAB2 . 'ESPI' . 'R' (J/kg/K dans le SI, FLOTTANT) qui 
 617 :              contient la constante du gaz parfait
 618 :            - TAB2 . 'ESPI' . 'H0K' (J/kg, FLOTTANT) qui contient
 619 :              l'enthalpie de formation du gaz a 0K
 620 :          * TAB2 . 'SCALPASS' = si existante, noms des scalaires passifs 
 621 :            a transporter (LISTMOTS)
 622 : 
 623 :  LMOT1 : objet de type LISTMOTS
 624 :          Noms de composantes du resultat (RCHPO1)
 625 :          Il contient dans l'ordre suivant: le nom de la densite
 626 :          totale, de la quantite de mouvement, de l'energie totale 
 627 :          par unite de volume, des densites des constituants du gaz 
 628 :          (en TAB2 . 'ESPEULE'), des scalaires passifs 
 629 :          (en TAB2 . 'SCALPASS') fois la densite
 630 :  
 631 :  MCHAM1: MCHAML contenant la masse volumique, qui a comme SPG (support 
 632 :          geometrique) 'DOMA' MOD1 'FACEL' (une composante, 'SCAL')
 633 :          (voir operateur 'PRET')
 634 :  
 635 :  MCHAM2: MCHAML  contenant la vitesse  et les cosinus 
 636 :           directeurs du repere locale (n,t) dans le repere 
 637 :           global (x,y) (dans le cas 2D 6 composantes:
 638 :           * 'UN' = vitesse normale  (SPG = 'DOMA' MOD1 'FACEL')
 639 :           * 'UT' = vitesse tangentielle  (SPG = 'DOMA' MOD1 FACEL')
 640 :           * 'NX' = n.x (SPG = 'FACE')
 641 :           * 'NY' = n.y (SPG = 'FACE')
 642 :           * 'TX' = t.x (SPG = 'FACE')
 643 :           * 'TY' = t.y (SPG = 'FACE')).
 644 :           (voir operateur 'PRET')
 645 :  
 646 :   MCHAM3: MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la pression du
 647 :           gaz (une  seule composante, 'SCAL').
 648 :           (voir operateur 'PRET')
 649 :  
 650 :  (MCHAM4):MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant les fractions 
 651 :           massiques des especes qui apparaissent explicitement dans 
 652 :           les equations d'Euler (voir TAB2 . 'ESPEULE') ;
 653 :           Dans le cas mono-espece, MCHAM4 n'est pas a donner
 654 :           (voir operateur 'PRET')
 655 : 
 656 :  (MCHAM5):MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant les scalaires
 657 :           transportes (noms des composantes en TAB2 . 'SCALPASS'') ;
 658 :           (voir operateur 'PRET')
 659 :  
 660 :   SORTIES
 661 :  
 662 :  
 663 :   RCHPO1: objet de type CHPOINT (composantes =  LMOT1)
 664 :           Residu  si MOT1 = 'RESI' (SPG = 'DOMA' MOD1 'CENTRE')
 665 :           Flux si MOT1 = 'FLUX'    (SPG = 'DOMA' MOD1 'FACE') 
 666 :  
 667 :   RFLOT1: objet de type FLOTTANT
 668 :           Il est le temps caracteristique associe a l'onde la plus
 669 :           rapide.
 670 :  
 671 :   Remarque
 672 :   --------
 673 :  
 674 :   RCHPO1 est egal a:
 675 :   * la derive temporelle des inconnues si l'option 'RESI' est utilisee
 676 :   * la composante normale a la face du flux convectif (les normales 
 677 :     etant donnees par 'DOMA' MOD1 'XXNORMAF') si l'option 'FLUX' 
 678 :     est utilisee
 679 : 
 680 : 
 681 : 
 682 :  IIIc : gaz  parfait, "Free matrix method"
 683 :  ____
 684 : 
 685 :  Discretisation en VF "cell-centered" des equations d'Euler pour un gaz 
 686 :  parfait.
 687 :  Inconnues: U (densite, quantite de mouvement, energie totale par unite 
 688 :  de volume).
 689 :  Implicitation via une methode "sans matrice" (a stoker), ou la 
 690 :  convergence a l'etat stationnaire est effectue en utilisant la methode 
 691 :  de Rusanov. Dans la cellule i-eme on doit calculer U_i^{n+1} t.q.
 692 : 
 693 :     (U_i^{n+1} - U_i^{n})  * AN_i(U^{n}) = 
 694 :         RES_i(U^{n}) + (BN_i(U^{n}) - BN_i(U_i^{n+1}))
 695 : 
 696 :  ou 
 697 :  AN_i contient les contributions liees au reciproque du pas de temps 
 698 :  (local) et  aux diffusivites numeriques (Rusanov) aux interfaces;
 699 :  RES_i contient le residu calcule avec n'import quelle methode 
 700 :  numerique (voir operateur KONV, cas IIIa).
 701 :  BN_i contient les contributions aux interfaces liees au flux 
 702 :  numerique centres et la diffusivite numerique (Rusanov) multiplie par
 703 :  l'etat voisin.
 704 :  On remarque que a l'etat stationnaire on a
 705 : 
 706 :     RES_i(U)=0
 707 : 
 708 :  Pour calculer le CHPOINT RES on utilise l'operateur KONV, cas IIIa.
 709 :  Pour calculer les CHPOINTs  AN et BN:
 710 : 
 711 :  AN = 'KONV' 'VF' 'PMONOFMM'  'AN'  LMOT1 MOD1  
 712 :       CHPO1 CHPO2 CHPO3 CHPO4 FLOT1 ('CLIM' LMOT2 CHPO5) ;
 713 : 
 714 :  BN = 'KONV' 'VF' 'PMONOFMM'  'BN'  LMOT1 MOD1  
 715 :       CHPO1 CHPO2 CHPO3 CHPO4 ('CLIM' LMOT2 CHPO5) ;
 716 :  
 717 :  ENTREES 
 718 :  
 719 :  MOD1  : objet MODELE 'EULER'
 720 : 
 721 :  LMOT1 : objet de type LISTMOTS
 722 :          Noms de composantes du resultat (RCHPO1)
 723 :          Il contient dans l'ordre suivant: le nom de la densite
 724 :          totale, de la quantite de mouvement, de l'energie totale 
 725 :          par unite de volume.
 726 :  
 727 :  CHPO1 : CHPOINT  contenant la masse volumique (en kg/m3; une 
 728 :          composante, 'SCAL').
 729 : 
 730 :  CHPO2 : CHPOINT  contenant les debits (en  kg/s/m2; deux 
 731 :          composantes en 2D, 'UX  ','UY  ', trois composantes
 732 :          en 3D, 'UX  ','UY  ', 'UZ  ').
 733 :      
 734 :  CHPO3 : CHPOINT  contenant l'energie totale par unite de volume 
 735 :          (en J/m3; une composante, 'SCAL').
 736 : 
 737 :  CHPO4 : CHPOINT  contenant le "gamma" du gaz (une composante, 
 738 :          'SCAL').
 739 : 
 740 :  FLOT1 : flottant, le double de la CFL
 741 : 
 742 :  LMOT2 : objet de type LISTMOTS
 743 :          Noms de composantes de conditions aux bords  (CHPO5)
 744 :          Il contient dans l'ordre suivant: le nom de la densite
 745 :          totale, de la vitesse, de la pression.
 746 : 
 747 :  CHPO5 : CHPOINT contenant les conditions aux bords (densite, 
 748 :          vitesse et pression sur le bord).
 749 :  
 750 :  SORTIES
 751 : 
 752 :  AN    : objet de type CHPOINT (SPG = 'DOMA' MOD1 'CENTRE', 
 753 :          une composante, 'SCAL').
 754 :  
 755 :  BN    : objet de type CHPOINT (SPG = 'DOMA' MOD1 'CENTRE', 
 756 :           composantes =  LMOT1)
 757 : 
 758 :  IIId : gaz  parfait, "Free matrix method", Euler/NS bas Mach
 759 :  ____
 760 : 
 761 :  Discretisation en VF "cell-centered" des equations d'Euler/NS pour un gaz 
 762 :  parfait.
 763 :  Inconnues: U (densite, quantite de mouvement, energie totale par unite 
 764 :  de volume).
 765 :  Implicitation via une methode "sans matrice" (a stoker), ou la 
 766 :  convergence a l'etat stationnaire est effectue en utilisant la methode 
 767 :  de Rusanov (voir IIIc). Dans la cellule i-eme on doit calculer 
 768 : 
 769 :     
 770 :     (U_i^{n+1} - U_i^{n}) = DUN 
 771 : 
 772 :  Syntaxe:
 773 : 
 774 :      DUN IPRO = 'KONV' 'VF' 'PMON1FMM' MOT1 LMOT1 MOD1
 775 :                  CHPORE  CHPO1 CHPO2 CHPO3 CHPO4 CHPO5 
 776 :                  FLOT1 FLOT2 
 777 :                  NJAC 'CLIM' LMOT2 CHPO6 CHPO7 ;
 778 : 
 779 :  ENTREES 
 780 :  
 781 :  MOT1  : objet de type MOT, methode d'inversion  ('PJACO', point Jacobi, 
 782 :          'LJACOF', 'LJACOB', 'LJACOFB' Gauss-Seidel Forward, Backward et
 783 :          Symmetric)
 784 : 
 785 :  LMOT1 : objet de type LISTMOTS
 786 :          Noms de composantes du resultat (RCHPO1)
 787 :          Il contient dans l'ordre suivant: le nom de la densite
 788 :          totale, de la quantite de mouvement, de l'energie totale 
 789 :          par unite de volume.
 790 :  
 791 :  MOD1  : objet MODELE 'EULER'
 792 : 
 793 :  CHPORE: CHPOINT  contenant le residu explicite (Euler/NS)
 794 : 
 795 :  CHPO1 : CHPOINT  contenant la masse volumique (en kg/m3; une 
 796 :          composante, 'SCAL').
 797 : 
 798 :  CHPO2 : CHPOINT  contenant les debits (en  kg/s/m2; deux 
 799 :          composantes en 2D, 'UX  ','UY  ', trois composantes
 800 :          en 3D, 'UX  ','UY  ', 'UZ  ').
 801 :      
 802 :  CHPO3 : CHPOINT  contenant l'energie totale par unite de volume 
 803 :          (en J/m3 en SI, une composante, 'SCAL').
 804 : 
 805 :  CHPO4 : CHPOINT contenant le "gamma" du gaz (une composante, 
 806 :          'SCAL').
 807 : 
 808 :  CHPO5 : CHPOINT contenat la vitesse de cut-off (une composante, 
 809 :          'SCAL').
 810 : 
 811 :  FLOT1 : flottant, le pas de temps physique 
 812 : 
 813 :  FLOT2 : flottant, le double de la CFL "duale"
 814 : 
 815 :  LMOT2 : objet de type LISTMOTS
 816 :          Noms de composantes de conditions aux bords  (CHPO6)
 817 :          Il contient dans l'ordre suivant: le nom de la densite
 818 :          totale, de la vitesse, de la pression.
 819 : 
 820 :  CHPO6 : CHPOINT contenant les conditions aux bords (densite, 
 821 :          vitesse et pression sur le bord).
 822 :  
 823 :  CHPO7 : CHPOINT contenant le rayon spectral visquex (une composante, 
 824 :          'SCAL').
 825 :  
 826 :  SORTIES
 827 : 
 828 :  DUN   : objet de type CHPOINT (SPG = 'DOMA' MOD1 'CENTRE', 
 829 :          composantes = LMOT1).
 830 :  
 831 :  IPRO  : entier, 0 si tout s'est bien passe
 832 : 
 833 : 
 834 :  IIIg. Gas parfait multi-constituant
 835 :  ___________________________________
 836 : 
 837 :  IIIg.1  Le residu 
 838 :  __________________
 839 : 
 840 : 
 841 :  RESIDU DELTAT = 'KONV' 'VF' 'PERFMULT' 'RESI' MOT1 LMOT1 MOD1 
 842 :                   MCHAM1 MCHAM2 MCHAM3 MCHAM4 MCHAM5
 843 :                   TABGAS (CHPO6) (CHPO7) (MAIL1) ;
 844 : 
 845 :  ENTRÉES 
 846 : 
 847 : 
 848 :  MOT1   : objet de type MOT
 849 :            Il indique la methode:
 850 :            'GODUNOV'  = solveur exacte
 851 :            'VANLEER'  = solveur de van Leer
 852 :            'VLH'      = solveur de van Leer Hanel
 853 :            'HUSVL'    = HUS (van Leer + Osher) 
 854 :            'HUSVLH'   = HUS (van Leer Hanel + Osher) 
 855 :            'AUSMPLUS' = AUSM+ 
 856 :            'ROE'      = solveur de Roe
 857 :            'SS'       = solveur choc-choc
 858 :            'AUSMPLM'  = AUSM+ low Mach
 859 :            'RUSANOV'  = solveur de Rusanov
 860 :            'RUSANOLM' = solveur de Rusanov pour le bas-Mach
 861 :            'CENTERED' = schema centre
 862 :            'ROELM'    = solveur de Roe-Turkel pour le bas-Mach
 863 :            'HLLC'     = solveur HLLC
 864 :            'HLLCLM'   = solveur HLLC-Turkel pour le bas-Mach
 865 :            'AUSMPUP'  = solveur AUSM+up low Mach
 866 : 
 867 :   
 868 :  LMOT1   : objet de type LISTMOTS
 869 :            Noms de composantes du resultat (RCHPO1)
 870 :            Il contient dans l'ordre suivant: le noms de la densite,
 871 :            de la qdm, de l'energie totale par unite de
 872 :            volume, de l'especes qui sont dans TABGAS.'ESPEULE'.
 873 : 
 874 :  MOD1    : objet MODELE.
 875 : 
 876 :  MCHAM1  : MCHAML contenant la masse volumique, qui a comme 
 877 :            SPG (support geometrique) 'DOMA' MOD1 'FACEL'
 878 :            (une composante, 'SCAL')
 879 :            (voir operateur PRET)
 880 : 
 881 :  MCHAM2  : MCHAML contenant la vitesse  et les cosinus 
 882 :            directeurs du repere locale (n,t) dans le repere 
 883 :            global (x,y) (dans le cas 2D 6 composantes:
 884 :            * 'UN' = vitesse normale  (SPG = 'DOMA' MOD1 'FACEL')
 885 :            * 'UT' = vitesse tangentielle  (SPG = 'DOMA' MOD1 FACEL')
 886 :            * 'NX' = n.x (SPG = 'FACE')
 887 :            * 'NY' = n.y (SPG = 'FACE')
 888 :            * 'TX' = t.x (SPG = 'FACE')
 889 :            * 'TY' = t.y (SPG = 'FACE')).
 890 :            (voir operateur PRET)
 891 : 
 892 :  MCHAM3  : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la pression du
 893 :            gaz (une  seule composante, 'SCAL').
 894 :            (voir operateur PRET)
 895 : 
 896 :  MCHAM4  : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant le "gamma" du 
 897 :            gaz (une seule composante, 'SCAL').
 898 :            (voir operateur PRET)
 899 : 
 900 :  MCHAM5 : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant les
 901 :           fractions massiques (nombre des composants egal 
 902 :           celui dans TABGAS.'ESPEULE') 
 903 : 
 904 :  TABGAS : la table contenant les properties de gas
 905 :           (voir PRET ou PRIM)
 906 : 
 907 :  CHPO6   : CHPOINT contenant la vitesse de cut-off 
 908 :            (SPG = 'DOMA' MOD1 'CENTRE', une seule composante, 
 909 :             'SCAL'). 
 910 :            A donner dans le cas bas Mach
 911 : 
 912 :  CHPO7   : CHPOINT contenant la deuxieme vitesse de cut-off 
 913 :            (SPG = 'DOMA' MOD1 'CENTRE', une seule composante, 
 914 :             'SCAL'). 
 915 :            A donner dans le cas bas Mach
 916 : 
 917 : (MAIL1)  : MAILLAGE de POI1 ou, si appartenant a  ('DOMA' MOD1 'FACE'),
 918 :            on ne calcule pas la contribution au flux ou au residu.
 919 : 
 920 :  SORTIES
 921 : 
 922 :  RESIDU  : objet de type CHPOINT (composantes =  LMOT1)
 923 :            Residu  si MOT2 = 'RESI' (SPG = 'DOMA' MOD1 'CENTRE')
 924 : 
 925 :  DELTAT  : objet de type FLOTTANT
 926 :            Il est le temps caracteristique associe a l'onde la plus
 927 :            rapide (meme dans le cas bas Mach, on considere le systeme 
 928 :            non preconditionne)
 929 : 
 930 :  Remarque
 931 :  --------
 932 : 
 933 :  RESIDU est egal a:
 934 :  * la derive temporelle des inconnues si l'option 'RESI' est utilisee
 935 : 
 936 : 
 937 : 
 938 :  IIIg.2   La matrice jacobienne du residu par rapport aux variables  
 939 :  _____    conservatives
 940 : 
 941 : 
 942 :  RMAT1  = 'KONV' 'VF' 'PERFMULT' 'JACOCONS' MOT1 MOD1 
 943 :           TABGAS LMOT1 (MAIL1) 
 944 :           CHPO1 CHPO2 CHPO3 CHPO4 (CHPO5 CHPO6) ;
 945 :  
 946 : 
 947 :  ENTRÉES 
 948 : 
 949 :  MOT1   : objet de type MOT
 950 :           'VLH'       : jacobien du residu pour la methode VLH 
 951 :           'AUSMPLUS'  : jacobien du residu pour la methode AUSM+ 
 952 :           'AUSMPLM'   : jacobien du residu pour la methode AUSM+ low Mach
 953 : 
 954 :  TABGAS : la table contenant les properties de gas
 955 :           (voir PRET ou PRIM)
 956 : 
 957 :  LMOT1   : objet de type LISTMOTS
 958 :            Noms des variables conservatives
 959 :            Il contient dans l'ordre suivant: le noms de la densite,
 960 :            de la qdm, de l'energie totale par unite de volume,
 961 :            de l'especes qui sont dans TABGAS.'ESPEULE'.
 962 : 
 963 :  MOD1    : objet MODELE.
 964 : 
 965 :  CHPO1   : CHPOINT contenant la masse volumique
 966 :            (SPG = 'DOMA' DOM1 'CENTRE', une seule composante, 
 967 :             'SCAL').
 968 : 
 969 :  CHPO2   : CHPOINT contenant la vitesse
 970 :            (SPG = 'DOMA' DOM1 'CENTRE', deux/trois composantes
 971 :             'UX', 'UY', 'UZ')
 972 : 
 973 :  CHPO3   : CHPOINT contenant la pression du gaz
 974 :            (SPG = 'DOMA' DOM1 'CENTRE', une seule composante, 
 975 :             'SCAL').
 976 : 
 977 :  CHPO4   : CHPOINT contenant les especes 
 978 :            (SPG = 'DOMA' DOM1 'CENTRE', nombre des composantes
 979 :             egal a celui dans TABGAS.'ESPEULE').
 980 : 
 981 :  CHPO5   : CHPOINT contenant la vitesse de cut-off
 982 :            (SPG = 'DOMA' DOM1 'CENTRE', une seule composante, 
 983 :             'SCAL'). 
 984 :            A donner dans le cas bas Mach
 985 : 
 986 :  CHPO6   : CHPOINT contenant la deuxieme vitesse de cut-off 
 987 :            (SPG = 'DOMA' DOM1 'CENTRE', une seule composante, 
 988 :             'SCAL'). 
 989 :            A donner dans le cas bas Mach
 990 : 
 991 :  (MAIL1) : MAILLAGE de POI1 ou, si appartenant a  ('DOMA' DOM1 'FACE'),
 992 :            on ne calcule pas la contribution au jacobien.
 993 : 
 994 :  SORTIES
 995 : 
 996 :  RMAT1   : objet de type MATRIK
 997 :            (SPG =  'DOMA' DOM1 'CENTRE')
 998 :            (inconnues primales = inconnues duales = LMOT1)
 999 :            Il contient le jacobien du residu par rapport aux variables 
1000 :            conservatives.
1001 : 
1002 :  IIIg.3   La contribution de quelque condition limite au residu et
1003 :  ------   a la matrice jacobienne
1004 : 
1005 : 
1006 :  RCHPLI RCHPRE = 'KONV' 'VF' 'PERFMULT' 'CLIM'  'RESI'
1007 :                  $MOD1 $MOD2 TABGAS LMOTC LMOTP
1008 :                  CHPRN CHPVN CHPPN CHPYN
1009 :                  CHPLI MOT1 ;
1010 : 
1011 :  ou 
1012 : 
1013 :  RJACO = 'KONV' 'VF' 'PERFMULT' 'CLIM'  'JACOCONS'
1014 :                  $MOD1 $MOD2 TABGAS LMOTC LMOTP
1015 :                  CHPRN CHPVN CHPPN CHPYN
1016 :                  CHPLI MOT1 ;
1017 : 
1018 : 
1019 :  $MOD1   : l'objet modele du domaine total
1020 : 
1021 :  $MOD2   : l'objet modele du domaine du bord 
1022 : 
1023 :  TABGAS : la table contenant les properties de gas
1024 :           (voir PRET ou PRIM)
1025 : 
1026 :  LMOTC   : LISTMOTS, noms des variables conservatives
1027 : 
1028 :  LMOTP   : LISTMOTS, noms des variables primitives
1029 : 
1030 :  CHPRN   : densite  (SPG = 'DOMA' $MOD1 'CENTRE', une seule
1031 :            composante, 'SCAL')
1032 : 
1033 :  CHPVN   : vitesse (SPG = 'DOMA' $MOD1 'CENTRE',
1034 :            composantes: 'UX', 'UY', ('UZ'))
1035 : 
1036 :  CHPPN   : pression  (SPG = 'DOMA' $MOD1 'CENTRE', une seule 
1037 :            composante, 'SCAL')
1038 : 
1039 :  CHPYN   : les especes  (SPG = 'DOMA' $MOD1 'CENTRE', le nombre 
1040 :            des composantes egal a celui dans TABGAS.'ESPEULE')
1041 : 
1042 :  MOT1    : MOT, type de condition limite:
1043 :            'INRI' : entree subsonique. On utilise les invariants de
1044 :                     Riemann. On donne densite 'RN', vitesse 'UX' 'UY' 
1045 :                     ('UZ'), pression 'PN', les fractions massiques
1046 :            'INSU' : entree subsonique. On donne l'enthalpie totale 'HT'
1047 :                     par unite de masse, l'entropie 'S' (pression divisee 
1048 :                     par la densite a la puissance gamma),les fractions massiques
1049 :                     On impose que 
1050 :                     la vitesse tangentielle est nulle. On recupere la 
1051 :                     pression a l'interieur. La contribution au residu et 
1052 :                     a la matrice jacobienne sont calculee via la methode
1053 :                     'AUSMPLUS'.
1054 :            'OUTP' : sortie subsonique. On donne la pression 'PN', 
1055 :                     on recupere la vitesse, la densite 
1056 :                     et les fractions massiques a l'interieur.
1057 :                     La contribution au residu et a la matrice jacobienne 
1058 :                     sont calculee via la methode 'AUSMPLUS'.
1059 :            'INJE' : condition limite d'injection compressible. 
1060 :                     On donne le flux de masse 'MOME', 'RT' (temperature fois
1061 :                     la constante de gaz) et les fractions massiques; 
1062 :                     on impose que la vitesse
1063 :                     tangentielle est nulle, on recupere la
1064 :                      pression.
1065 :            'INSS':  entree supersonique. On impose les
1066 :                     conditions limites dans le flux. 
1067 :                     On donne la densite, la vitesse,
1068 :                     la pression et les fractions massiques
1069 :            'OUTSS': sortie supersonique. On donn le
1070 :                     Champ par points vide.
1071 :                     Toute l'information est recoupere de 
1072 :                     l'interieur.  
1073 :            'RESE' : condition limite de type reservoir. 
1074 :                     On donne la pression 'PN' la densite 'RN' et les 
1075 :                     fractions massiques ; on impose que la vitesse
1076 :                     tangentielle est nulle, on recupere la
1077 :                     pression interieure et, selon sa valeur,
1078 :                     on impose les conditions de col sonique ou 
1079 :                     subsonique.
1080 : 
1081 :  CHPLI   : condition limite impose (SPG ='DOMA' $MOD2 'CENTRE')
1082 :            Les composantes dependent de MOT1
1083 : 
1084 :  Resultats
1085 : 
1086 :  RCHPLI : CHPOINT qui contient les valeur de la densite, de la vitesse, 
1087 :           de pression et des fractions massiques associe a la condition limite.
1088 :           (SPG ='DOMA' $MOD2 'CENTRE', composantes = LMOTP)
1089 : 
1090 :  RCHPRE : CHPOINT qui contient la contribution de la condition limite au 
1091 :           residu (SPG en 'DOMA' $MOD1 'CENTRE', composantes = LMOTC)
1092 : 
1093 :  RJACO  : MATRIK qui contient la contribution de la condition limite a la
1094 :           matrice jacobienne du residu.
1095 : 
1096 : 
1097 :  IV Transport des scalaires
1098 :  __________________________
1099 : 
1100 : 
1101 :  Discretisation en VF "cell-centered" du transport des scalaires
1102 :           --> -->
1103 : d/dt S + div . (u   S) = 0
1104 : 
1105 :  RCHPO1 RFLOT1 = 'KONV' 'VF' 'CLAUDEIS' 'FACE' MOT1 MOT2
1106 :                         MOD1 CHPO1 MCHAM1 ; 
1107 :   
1108 :  or
1109 : 
1110 :  RMAT1 = 'KONV' 'VF' 'CLAUDEIS' 'FACE' 'JACO' MOT2
1111 :                         MOD1 CHPO1 MCHAM1 ;
1112 : 
1113 :  ENTRÉES 
1114 : 
1115 :  MOT1    : objet de type MOT
1116 :            Il vaut 'RESI' si on veut calculer le residu,
1117 :            (i.e. l' increment des scalaires passives par unite
1118 :            de temps)
1119 :            Il vaut 'FLUX' si on veut calculer le flux aux
1120 :            interfaces
1121 : 
1122 :  MOT2    : objet de type MOT
1123 :            Il indique la methode de calcul du flux
1124 :            Pour l'instant
1125 :            'UPWIND'
1126 :            'CENTERED'
1127 : 
1128 :  MOD1    : objet MODELE.
1129 : 
1130 :  CHPO1   : CHPOINT contenant la vitesse aux FACEs
1131 :            (SPG = 'DOMA' MOD1 'FACE')
1132 :            2/3 composantes, 'UX', 'UY', 'UZ'
1133 : 
1134 :  MCHAM1  : MCHAML contenant les scalaires a transporter
1135 :            SPG (support geometrique) 'DOMA' MOD1 'FACEL' 
1136 : 
1137 :  SORTIES
1138 : 
1139 :  RCHPO1  : objet de type CHPOINT (memes composantes que MCHAM1)
1140 :            Residu si MOT2 = 'RESI' (SPG = 'DOMA' MOD1 'CENTRE')
1141 :            Flux si MOT2 = 'FLUX'  (SPG = 'DOMA' MOD1 'FACE')
1142 : 
1143 :  RFLOT1  : objet de type FLOTTANT
1144 :            Il est le temps caracteristique associe a l'onde la plus
1145 :            rapide.
1146 : 
1147 :  RMAT1   : objet de type MATRIK
1148 :            (SPG =  'DOMA' MOD1 'CENTRE')
1149 :            (inconnues primales = inconnues duales; memes composantes 
1150 :             que MCHAM1)
1151 : 
1152 : 
1153 :  V Flux bifluide
1154 :  ___________________________________________________
1155 : 
1156 : 
1157 :  Discretisation en VF "cell-centered" des equations bifluides pour un
1158 :  melange d'eau et air.
1159 : 
1160 :  Inconnues:
1161 : 
1162 :  mass volumique, quantite de mouvement (qdm) et energie totale par 
1163 :  unite de volume (variables conservatives) de chaque fluide
1164 : 
1165 :  ou
1166 : 
1167 :  taux de vide, vitesses du gaz et du liquide, pression,
1168 :  temperature du gaz et du liquide (variables primitives)
1169 : 
1170 :  On peut calculer:
1171 : 
1172 :  RCHPO2 RCHPO1 RFLOT1 = 'KONV' 'VF' 'TWOFLUID' MOT1 MOT2  MOD1 LMOT1 
1173 :       MCHAM1 MCHAM2 MCHAM3 MCHAM4 MCHAM5 MCHAM6 MCHAM7 MCHAM8;
1174 : 
1175 :   ENTRÉES 
1176 : 
1177 :  MOT1   : objet de type MOT
1178 :            Il vaut 'RESI' si on veut calculer le residu
1179 :            Il vaut 'FLUX' si on veut calculer le flux
1180 : 
1181 :  MOT2   : object de type MOT
1182 :           Il indique la methode:
1183 :            'AUSMP1'   = AUSM+
1184 :            'AUSMP2'   = AUSM+ preconditione
1185 :            'AUSMDV1'  = AUSMDV
1186 :            'AUSMDV2'  = AUSMDV preconditione
1187 :   
1188 :  LMOT1   : object de type LISTMOTS
1189 :            Noms de composantes du resultat (RCHPO1)
1190 :            Il contient dans l'ordre suivant: le noms de la densite,
1191 :            de la qdm, 
1192 :              de la masse volumique du gaz,
1193 :              de la masse volumique du liquide,
1194 :              de la quantite de mouvement du gaz,
1195 :              de la quantite de mouvement du liquide,
1196 :              de l'energie totale du gaz par unite de volume           
1197 :              de l'energie totale du liquide par unite de volume           
1198 : 
1199 :  MOD1    :  object MODELE
1200 : 
1201 :  MCHAM1  : MCHAML contenant la taux de vide, qui a comme 
1202 :            SPG (support geometrique) 'DOMA' MOD1 'FACEL'
1203 :            (une composante, 'SCAL')
1204 :            (voir operateur PRET)
1205 : 
1206 :  MCHAM2  : MCHAML contenant la vitesse du gaz et les cosinus 
1207 :            directeurs du repere locale (n,t) dans le repere 
1208 :            global (x,y) (dans le cas 2D 6 composantes:
1209 :            * 'UN' = vitesse normale  (SPG = 'DOMA' MOD1 'FACEL')
1210 :            * 'UT' = vitesse tangentielle  (SPG = 'DOMA' MOD1 FACEL')
1211 :            * 'NX' = n.x (SPG = 'FACE')
1212 :            * 'NY' = n.y (SPG = 'FACE')
1213 :            * 'TX' = t.x (SPG = 'FACE')
1214 :            * 'TY' = t.y (SPG = 'FACE')).
1215 :            (voir operateur PRET)
1216 : 
1217 :  MCHAM3  : MCHAML contenant la vitesse du liquide et les cosinus 
1218 :            directeurs du repere locale (n,t) dans le repere 
1219 :            global (x,y) (dans le cas 2D 6 composantes:
1220 :            * 'UN' = vitesse normale  (SPG = 'DOMA' MOD1 'FACEL')
1221 :            * 'UT' = vitesse tangentielle  (SPG = 'DOMA' MOD1 FACEL')
1222 :            * 'NX' = n.x (SPG = 'FACE')
1223 :            * 'NY' = n.y (SPG = 'FACE')
1224 :            * 'TX' = t.x (SPG = 'FACE')
1225 :            * 'TY' = t.y (SPG = 'FACE')).
1226 :            (voir operateur PRET)
1227 : 
1228 :  MCHAM4  : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la pression
1229 :            (une  seule composante, 'SCAL').
1230 :            (voir operateur PRET)
1231 : 
1232 :  MCHAM5  : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la temperature du
1233 :            gaz (une  seule composante, 'SCAL').
1234 :            (voir operateur PRET)
1235 : 
1236 :  MCHAM6  : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la temperature du
1237 :            liquide (une  seule composante, 'SCAL').
1238 :            (voir operateur PRET)
1239 : 
1240 :  MCHAM7  : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la densite du
1241 :            gaz (une  seule composante, 'SCAL').
1242 :            (voir operateur PRET)
1243 : 
1244 :  MCHAM8  : MCHAML (SPG = 'DOMA' MOD1 'FACEL') contenant la densite du
1245 :            liquide (une  seule composante, 'SCAL').
1246 :            (voir operateur PRET)
1247 : 
1248 : 
1249 :  SORTIES
1250 :  -------
1251 : 
1252 :  RCHPO1  : objet de type CHPOINT (composantes =  LMOT1)
1253 :            Residu si MOT2 = 'RESI' ou flux si MOT2 = 'FLUX'
1254 :           (SPG = 'DOMA' MOD1 'CENTRE')
1255 : 
1256 :  RFLOT1  : objet de type FLOTTANT
1257 :            Il est le temps caracteristique associe a l'onde la plus
1258 :            rapide
1259 : 
1260 : 
1261 :  Vl  Ecoulement Reactif avec resolution  "Discrete Equation Method".
1262 :  ___________________________________________________________________
1263 : 
1264 : 
1265 :  
1266 :  Discretisation en VF "cell-centered" des equations d'Euler pour
1267 :   un gaz parfait multi-constituent. Methode "DEM".
1268 : 
1269 : 
1270 :   RCHPO1 RFLOT1 = 'KONV' 'VF' 'DEM' MOT1 MOT2 MOT3 MOD1 TABG LMOT1  
1271 :                                     CHPA1 CHPA2 
1272 :                                     MCHAA1  MCHAA2 
1273 :                                     MCHAR1  MCHAR2
1274 :                                     MCHAV1  MCHAV2
1275 :                                     MCHAP1  MCHAP2 
1276 :                                     GRALP1 K0 EPS  MAILLIM 
1277 :                                    (CHPV1 CHPV2) ;
1278 : 
1279 :   ARGUMENTS :
1280 :   ----------
1281 : 
1282 :    MOT1   : objet de type MOT
1283 :             Il vaut  'RESI' si on veut calculer le
1284 :             residu
1285 : 
1286 :    MOT2   : objet de type MOT
1287 :            Il indique la methode de decentrement (partie non reactive)
1288 :            'SS'     =  solveur choc-choc 
1289 :            'VLH'    = solver de Van Leer Hanel
1290 :            'AUSMPUP'= solver AUSM+up (bas Mach)
1291 : 
1292 :    MOT3   : objet de type MOT
1293 :             'CONS'  = la vitesse fondamentale est constante (objet REEL)
1294 :             'VARI'  = la vitesse fondamentale est variable  (objet CHAMPOINT)
1295 : 
1296 :    MOD1   : objet modele, type EULER
1297 :  
1298 :    TABG   : objet de type TABLE;
1299 :            la table contenant les properties de gas
1300 :            (voir PRET ou PRIM)
1301 :      
1302 :    LMOT1   : objet de type LISTMOTS
1303 :             contenant les noms des composants du vecteur RCHPO1.
1304 :             Ils sont donnes dans l'ordre suivant : noms de alpha, 
1305 :             de la  densite, de la vitesse, de
1306 :             l'energie specifique totale.
1307 : 
1308 :    CHPA1   : CHPOINT contenant la fraction volumique \alpha de
1309 :             la phase 1, qui a comme SPG (support geometric)   
1310 :             l'index 'CENTRE' de la table MOD1 (une
1311 :             composante, 'SCAL')
1312 : 
1313 :    CHPA2   : CHPOINT contenant la fraction volumique \alpha de
1314 :             la phase 2, qui a comme SPG (support geometric)   
1315 :             l'index 'CENTRE' de la table MOD1 (une
1316 :             composante, 'SCAL')
1317 : 
1318 :    MCHAA1  : MCHAML contenant la fraction volumique \alpha de
1319 :             la phase 1, qui a comme SPG (support geometric)   
1320 :             l'index 'FACEL' de la table MOD1 (une
1321 :             composante, 'SCAL')
1322 : 
1323 :   MCHAA2  : MCHAML contenant la fraction volumique \alpha de
1324 :             la phase 2, le meme SPG que pour MCHAA1, (une
1325 :             composante, 'SCAL') 
1326 : 
1327 :   MCHAR1  : MCHAML contenant la densite de la
1328 :             phase 1, le meme SPG que pour MCHAA1, (une
1329 :             composante, 'SCAL') 
1330 :                
1331 :   MCHAR2  : MCHAML contenant la densite de la
1332 :             phase 2, le meme SPG que pour MCHAA1, (une
1333 :             composante, 'SCAL') 
1334 : 
1335 : 
1336 :   MCHAV1  : MCHAML contenant la vitesse du liquide de la
1337 :             phase 1  et les cosinus 
1338 :             directeurs du repere locale (n,t) dans le repere 
1339 :             global (x,y) (dans le cas 2D 6 composantes:
1340 :             * 'UN' = vitesse normale  (SPG = 'DOMA' MOD1 'FACEL')
1341 :             * 'UT' = vitesse tangentielle  (SPG = 'DOMA' MOD1 FACEL')
1342 :             * 'NX' = n.x (SPG = 'FACE')
1343 :             * 'NY' = n.y (SPG = 'FACE')
1344 :             * 'TX' = t.x (SPG = 'FACE')
1345 :             * 'TY' = t.y (SPG = 'FACE')).
1346 :             (voir operateur PRET)
1347 : 
1348 : 
1349 :  MCHAV2  : MCHAML contenant la vitesse du liquide de la
1350 :             phase 2  et les cosinus 
1351 :             directeurs du repere locale (n,t) dans le repere 
1352 :             global (x,y) (la meme structure que pour MCHAV1)
1353 : 
1354 : 
1355 : 
1356 :  MCHAP1  : MCHAML (SPG =('DOMA' MOD1 'FACEL')) contenant la
1357 :            la pression de la phase 1
1358 :            (une composante, 'SCAL').
1359 : 
1360 :  MCHAP2  : MCHAML (SPG =('DOMA' MOD1 'FACEL')) contenant la
1361 :            la pression de la phase 2
1362 :            (une composante, 'SCAL'). 
1363 : 
1364 :  K0      : FLOTTANT/CHAMPOINT (voir MOT3), vitesse fondamentale de la 
1365 :            flamme
1366 :  
1367 :  GRALP1  : CHPOINT, grad(alp1)/|grad(alp1)|
1368 : 
1369 :  EPSILON : FLOTTANT t.q. a < EPSILON => a = 0
1370 : 
1371 :  MAILLIM : MAILLAGE -- la partie du maillage ou le flux
1372 :                        n'est pas determine ; il est
1373 :                        calcule via les procedures pour
1374 :                        les Conditions Limites.
1375 : 
1376 :  CHPV1   : CHPOINT, cut-off speed dans la methode AUSM+up,
1377 :            phase 1
1378 :  
1379 :  CHPV2   : CHPOINT, cut-off speed dans la methode AUSM+up,
1380 :            phase 2
1381 :  
1382 :  RESULTS :
1383 : -----------
1384 : 
1385 :  RCHPO1  : objet de type CHPOINT (le noms des composantes
1386 :            sont dans LMOT1)
1387 :            Residu  MOT2 = 'RESI' (SPG =('DOMA' MOD1 'CENTRE'))
1388 : 
1389 :  RFLOT1  : objet de type FLOTTANT
1390 :            Il est le temps caracteristique associe a l'onde la plus
1391 :            rapide
1392 : 
1393 :  VII Ghost fluid method for the poor.
1394 :  ___________________________________________________
1395 : 
1396 : 
1397 :  Discretisation en VF "cell-centered" des equations bifluides pour le
1398 :  transport d'une interface.
1399 : 
1400 :  Inconnues: phi, densite, quantite de mouvement, energie totale par
1401 :  unite de volume, densites et fractions volumiques des constituants. 
1402 : 
1403 : 
1404 :  RCHPO1 RFLOT1 = 'KONV' 'VF' 'GFMP' MOT1 MOT2 MOD1 
1405 :                         TABG
1406 :                         LMOT1  MCHAPH MCHAR MCHAV MCHAP (MCHAY 
1407 :                         MCHAA) MCHPPH  LOG1 MAILLIM ;
1408 :  ENTREES 
1409 : 
1410 : 
1411 :  MOT1   : objet de type MOT
1412 :           Il vaut 'RESI' si on veut calculer le residu
1413 : 
1414 :  MOT2   : objet de type MOT
1415 :           Il indique la methode de decentrement:
1416 :           'GODUNOV' 
1417 : 
1418 :  MOD1   : objet modele de type Euler
1419 : 
1420 :  TABG   : objet de type TABLE
1421 :           il contient les proprietes du gaz
1422 :           (voir operateur PRIM)
1423 :  
1424 :  LMOT1   : objet de type LISTMOTS
1425 :           Noms de composantes du resultat (RCHPO1)
1426 :           Il contient dans l'ordre suivant les noms de
1427 :            phi, de la densites, de la vitesse, de l'energie totale 
1428 :            par unite de volume, des densites des especes, des
1429 :            fractions volumiques.
1430 : 
1431 :  MCHAPH  : MCHAML contenant la fonction phi, 
1432 :            qui a comme SPG (support geometrique) l'indice 'FACEL' de la
1433 :            table associee a MOD1 (une composante, 'SCAL')
1434 : 
1435 :  MCHAR   : MCHAML contenant la masse volumique, 
1436 :            meme SPG que MCHAA1, une composante, 'SCAL'
1437 : 
1438 :  MCHAV   : MCHAML  contenant la vitesse et les cosinus
1439 :            directeurs du repere locale (n,t) dans le repere
1440 :            global (x,y), meme SPG que MCHAA1, 
1441 :            (dans le cas 2D 6 composantes:
1442 :            * 'UN' = vitesse normale  (SPG =('DOMA' MOD1 'FACEL'))
1443 :            * 'UT' = vitesse tangentielle  (SPG =('DOMA' MOD1 'FACEL'))
1444 :            * 'NX' = n.x (SPG = 'FACE')
1445 :            * 'NY' = n.y (SPG = 'FACE')
1446 :            * 'TX' = t.x (SPG = 'FACE')
1447 :            * 'TY' = t.y (SPG = 'FACE')).
1448 : 
1449 :  MCHAP   : MCHAML (SPG =('DOMA' MOD1 'FACEL')) contenant la pression 
1450 :           (une  seule composante, 'SCAL').
1451 : 
1452 :  MCHAY   : MCHAML (SPG =('DOMA' MOD1 'FACEL')) contenant les fractions
1453 :            massiques (composantes dans TABG . 'ESPEULE').
1454 : 
1455 :  MCHAA   : MCHAML (SPG =('DOMA' MOD1 'FACEL')) contenant les fractions
1456 :            volumiques (composantes dans TABG . 'ESPEULE').
1457 : 
1458 :  MCHPPH  : CHPO contenant phi,
1459 :            qui a comme SPG (support geometrique) l'indice 'CENTRE' de la
1460 :            table associee a MOD1 (une composante, 'SCAL')
1461 : 
1462 :  LOG1    : LOGIQUE, si VRAI phi et les fractions volumiques sont
1463 :            traites avec une approche conservative
1464 : 
1465 :  MAILLIM : MAILLAGE -- points faces ou le flux n'est pas calcule.
1466 : 
1467 :  SORTIES
1468 : 
1469 :  RCHPO1  : objet de type CHPOINT (composantes =  LMOT1)
1470 :            Residu  si MOT2 = 'RESI' (SPG =('DOMA' MOD1 'CENTRE'))
1471 : 
1472 :  RFLOT1  : objet de type FLOTTANT
1473 :            Il est le temps caracteristique associe a l'onde la plus
1474 :            rapide.
1475 : 
1476 : 
1477 : 

© Cast3M 2003 - Tous droits réservés.
Mentions légales