Notices | Cast3M

$$$$ FLAM     NOTICE  CHAT      11/09/12    21:16:15     7124           
                                             DATE     11/09/12
                                                                                
  Operateur FLAM                           Voir aussi :
    --------------  

 1a) ZH2 ZO2 ZN2 ZH2O Q = FLAM | 'EBU'      | RHO CV MOT1 R YH2 YO2 YN2 
                              | 'LAMINAIR' | 


                YH2O T Dt | K E YH2u C_EBU (C0)    | ;
                          |                        |


 1b) RCHPO1 RCHPO2 RCHPO3 RCHPO4  = 'FLAM' 'HEAVYSID' CHPO1 FLOT1 FLOT2 
                                    FLOT3 CHPO2 CHPO3 CHPO4 CHPO5 ;
                  
 1c) RCHPO1 RCHPO2 RCHPO3 RCHPO4  = 'FLAM' 'ARRHENIU' CHPO1 FLOT1 FLOT2
                                    FLOT3 FLOT4 FLOT5 FLOT6 FLOT7 FLOT8
                                    CHPO2 CHPO3 CHPO4 CHPO5 ;

 
 2a) RCHPO1 = 'FLAM' 'CREBCOM' MOD1 FLOT1 FLOT2 CHPO1 ;

 2b) RCHPO1 RCHPO2 = 'FLAM' 'CREBCOM2' MOD1 TAB2 LMOT1 LREE1 CHPO1 
     CHPO2 CHPO3 CHPO4 CHPO5 FLOT1 FLOT2 FLOT3 ;
 


    OBJET :
    _______

 1) Cet operateur integre du temps tn au temps tn+1=tn+Dt des equations
 differentielles ordinaires qui modelisent l'evolution temporelle des
 fractions massiques d'hydrogene, d'oxygene, d'azote et de vapeur d'eau
 selon la reaction globale,

        2 H_2 + O_2 ---> 2 H_2O

 Ainsi, on resoud entre tn et tn+1,

        dYH2/dt  = C_H2  omega(...)
        dYO2/dt  = C_O2  omega(...)
        dYH2O/dt = C_H2O omega(...)

 ou Yi est la fraction massique associee a l'espece i et les C_i 
 sont des constantes dependantes de la stoechiometrie, et des masses 
 molaires des differents constituants, et omega la vitesse de reaction. 
 L'azote intervient uniquement comme espece neutre.
 Suivant l'option, la vitesse de reaction est donnee par une cinetique
 adaptee a la combustion en regime laminaire, a la combustion en regime
 turbulent de type Eddy Break-Up (EBU), ou une loi d'Arrhenius adaptee
 a la detonation (ARRHENIU), ou un modele de combustion infiniment
 rapide si la temperature est plus grande qu'une temperature seuil
 (HEAVYSID).


 2) Suivi du front de flamme selon le critere "CREBCOM"

   Commentaires :
   ______________

1a) Cas 'LAMINAR', 'EBU'

   ZH2   : fraction massique de H2 au temps tn+1 (CHPOINT)
   ZO2   : fraction massique de O2 au temps tn+1 (CHPOINT)
   ZN2   : fraction massique de N2 au temps tn+1 (CHPOINT)
   ZH2O  : fraction massique de H2O au temps tn+1 (CHPOINT)
   Q     : energie liberee par la reaction (CHPOINT)
   RHO   : densite du melange (CHPOINT)
   CV    : capacite calorifique a volume constant (CHPOINT)
   MOT1  : dependance en T des CV ('LINEAIRE' ou 'QUADRATI')
   R     : constante des gaz du melange (CHPOINT)
   YH2   : fraction massique de H2 au temps tn (CHPOINT)
   YO2   : fraction massique de O2 au temps tn (CHPOINT)
   YN2   : fraction massique de N2 au temps tn (CHPOINT)
   YH2O  : fraction massique de H2O au temps tn (CHPOINT)
   T     : temperature du melange (CHPOINT)
   Dt    : pas de temps (REEL)
   K     : energie cinetique turbulente (CHPOINT)
   E     : taux de dissipation de K (CHPOINT)
   YH2u  : fraction massique de H2 initiale (CHPOINT)
   C_EBU : constante du modele Eddy Break-Up (REEL)
   C0    : valeur seuil (defaut 1.d-4) (REEL) :
           Yh2/Yh2u > 1-C0 -> Q=0 ; Yh2/Yh2u < C0 -> Yh2=0

   N.B.: Tous les CHPOINTs ont une seule composante et le meme support
         geometrique : CENTRE ou SOMMET.

1b) Cas 'HEAVYSID'

   CHPO1   : CHPOINT contenant la temperature de seuil (en K; 
             une composante, 'SCAL').

   FLOT1   : FLOTTANT contenant l'enthalpie de H2 a 0K (en J/kg)

   FLOT2   : FLOTTANT contenant l'enthalpie de O2 a 0K (en J/kg)

   FLOT3   : FLOTTANT contenant l'enthalpie de H2O a 0K (en J/kg)
      
       
   CHPO2   : CHPOINT  contenant la masse volumique de
             l'hydrogene (en Kg/m^3; une composante, 'H2  ').

   CHPO3   : CHPOINT  contenant la masse volumique de
             l'oxygene (en Kg/m^3; une composante, 'O2  ').

   CHPO4   : CHPOINT  contenant la masse volumique de
             l'eau (en Kg/m^3; une composante, 'H2O ').

   CHPO5   : CHPOINT contenant la temperature (en K; une 
             composante, 'SCAL').

      
   RCHPO1  : CHPOINT  contenant la masse volumique de 
            l'hydrogene apres la reaction (en Kg/m^3; 
            une composante, 'H2  ').

   RCHPO2  : CHPOINT  contenant la masse volumique de
             l'oxygene  apres la reaction (en Kg/m^3;
              une composante, 'O2  ').

   RCHPO3  : CHPOINT  contenant la masse volumique de
             l'eau  apres la reaction (en Kg/m^3;  
             une composante, 'H2O ').

   RCHPO4  : CHPOINT contenant la chaleur libere' (en J/m^3; 
             une composante, 'SCAL').

   N.B.: Tous les CHPOINTs ont une seule composante et le meme support
         geometrique. 

1c) CAS 'ARRHENIU'

   CHPO1   : CHPOINT contenant la temperature de seuil (en K; 
             une composante, 'SCAL').

   FLOT1   : FLOTTANT contenant la constant A (en unites S.I.)

   FLOT2   : FLOTTANT contenant la constant b 

   FLOT3   : FLOTTANT contenant la constant c 

   FLOT4   : FLOTTANT contenant la "temperature" d'activation
             Ta (en K)

   FLOT5   : FLOTTANT contenant l'enthalpie de H2 a 0K (en J/kg)

   FLOT6   : FLOTTANT contenant l'enthalpie de O2 a 0K (en J/kg)

   FLOT7   : FLOTTANT contenant l'enthalpie de H2O a 0K (en J/kg)
      
   FLOT8   : FLOTTANT contenant le pas de temps (en s)

   CHPO2   : CHPOINT  contenant la masse volumique de
             l'hydrogene (en Kg/m^3; une composante, 'H2  ').

   CHPO3   : CHPOINT  contenant la masse volumique de
             l'oxygene (en Kg/m^3; une composante, 'O2  ').

   CHPO4   : CHPOINT  contenant la masse volumique de
            l'eau (en Kg/m^3; une composante, 'H2O ').

   CHPO5   : CHPOINT contenant la temperature (en K; une 
             composante, 'SCAL').
     
   RCHPO1  : CHPOINT  contenant la masse volumique de 
             l'hydrogene apres la reaction (en Kg/m^3; 
             une composante, 'H2  ').

   RCHPO2  : CHPOINT  contenant la masse volumique de
            l'oxygene apres la reaction (en Kg/m^3; 
            une composante, 'O2  ').

   RCHPO3  : CHPOINT  contenant la masse volumique de
             l'eau  apres la reaction (en Kg/m^3; 
             une composante, 'H2O ').

   RCHPO4  : CHPOINT contenant la chaleur libere Q (en J/m^3; 
             une composante, 'SCAL').

   N.B.:-tous les CHPOINTs ont une seule composante et le meme support 
         geometrique. 


2a) Cas 'CREBCOM'

    MOD1   : Objet MODELE.

    FLOT1  : parametre "epsilon" du critere CREBCOM 

    FLOT2  : valeur maximum de la variable de progres de la reaction
 
    CHPO1  : variable de progres de la reaction (objet de type 
             CHAMPOINT)
    
    RCHPO1 : objet de type CHAMPOINT qui vaut 1 dans la region de 
             combustion, 0 ailleurs (une composante, 'SCAL')


2b) Cas 'CREBCOM2'
    
    MOD1   : objet MODELE.
   
    TAB2   : table des proprietes du gaz (voir operateur 'PRIM')

    LMOT1  : LISTMOTS qui contient les noms des composantes qui
             interviennent dans la reaction chimique consideree

    LREE1  : LISTREEL qui contient le coefficients stoechiometriques de 
             la reaction chimique consideree (positifs pour les 
             reactifs, negatifs pour les produits; meme ordre que LMOT1)

    CHPO1  : CHPOINT  contenant la masse volumique (une composante,
             'SCAL', support geometrique SPG = 'DOMA' MOD1 'CENTRE')

    CHPO2  : CHPOINT contenant les fractions massiques des especes en 
             LMOT1 (SPG = 'DOMA' MOD1 'CENTRE')

    CHPO3  : CHPOINT contenant la fraction massique initiale de la 
             premiere espece en LMOT1 (meme nom que la premiere espece 
             en LMOT1, SPG = 'DOMA' MOD1 'CENTRE')

    CHPO4  : CHPOINT contenant la fraction massique finale de la premiere 
             espece en LMOT1 (meme nom que la premiere espece en LMOT1, 
             SPG = 'DOMA' MOD1 'CENTRE')

    CHPO5  : CHPOINT contenant la vitesse caracteristique (une composante, 
             'SCAL', SPG = 'DOMA' MOD1 'CENTRE')

    CHPO6  : CHPOINT contenant la dimension de la maille (une composante, 
             'SCAL', SPG = 'DOMA' MOD1 'CENTRE')

    FLOT1  : FLOTTANT, parametre "epsilon" du critere CREBCOM

    FLOT2  : FLOTTANT, pas de temps

    FLOT3  : FLOTTANT, parametre pour controler la variable de progres. 
             Si la variable de progres est plus petite que (-FLOT3) ou 
             plus grande que (1+FLOT3), un message d'erreur est donne

    RCHPO1 : CHPOINT, increment de l'energie totale par unite de volume
             suite a la reaction chimique (une composante, 'SCAL',
             SPG = 'DOMA' MOD1 'CENTRE')
 
    RCHPO2 : CHPOINT, increment des densite massiques des especes suite 
             a la reaction chimique (noms de composantes en LMOT1, 
             SPG = 'DOMA' MOD1 'CENTRE')

    
    Remarques :
    ___________

    'LAMINAIR' 
        Combustion hydrogene en regime laminaire, avec la vitesse
        de reaction donnee par une cinetique globale,

        omega = C1 RHO^2 YH2 YO2 exp(-C2/T)

        ou C1 et C2 sont des constantes predefinies.

    'EBU'
        Combustion hydrogene en regime turbulent, avec la vitesse
        de reaction donnee par un modele Eddy Break-Up (EBU),
        
        omega = C_EBU (E/K) YH2 ( 1 - YH2/YH2u )/(2MH2)

        ou MH2 est la masse molaire de l'hydrogene.

     'HEAVYSID'
        Combustion hydrogene avec la vitesse de reaction donnee par
        une loi de la forme:
       
        omega = +oo         si (T>Ts)
        omega = 0           sinon

        i.e. si T>Ts on a une combustion infiniment rapide.

    'ARRHENIU'
        Combustion hydrogene avec la vitesse de reaction donnee par
        une loi d'Arrhenius generale, de la forme:

        omega = A T^-b (RHO_H2)^c (RHO_O2)  exp (-Ta/T)  si (T>Ts)
        omega = 0                                        sinon

    'CREBCOM2'
        Dans le cas 'CREBCOM2' on controle la positivite de la densite,
        des fractions massiques, de la vitesse caracteristique, que la 
        variable de progres est comprise entre 0 et 1.